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EXECUTIVE SUMMARY 


Aircraft certification requires -he evaluation of the effects of ice accretions on aircraft 
aerodynamic components such as wings, tails, and control surfaces. Ice accumulation on 
aerodynamic surfaces can have a significant impact on aircraft performance handling qualities, 
and thus, aircraft safety. In general, ice accumulation on aerodynamic surfaces can cause flow 
separation whose extent over the aerodynamic surface is t function of ice shape and wing 
geometry. The term critical is often used to identify ice shapes responsible for large degradation 
in the aerodynamic performance of lifting surfaces. 

A number of experimental studies have been conducted over the years in an effort to assess the 
effect of various forms of ice accretions on aircraft aerodynamic performance and handling 
qualities. Most of these studies, however, have been limited to two-dimensional airfoil sections 
due to tunnel and model cost constraints. Only a small number of investigations have addressed 
the impact of ice shapes on iced three-dimensional (3D) finite wings. 

To address the lack of experimental data for iced 3D finite wing configurations, the Federal 
Aviation Administration (FAA), Wichita State University (WSU), the National Aeronautics and 
Space Administration (NASA), and general aviation aircraft manufacturers initiated a 
collaborative research program in the fall of 200C for a systematic evaluation of ice accretion 
effeCiS on finite wings. The main objective of this research program was the development of a 
3D experimental database of ice accretion effects on a swept 3D finite wing. The database can 
be used for developing certification guidance material and for improving and validating 
simulation tools for aerodynamic analysis and design. A research grant was awarded to WSU to 
design and fabricate a wind tunnel model and to conduct icing and aerodynamic performance 
tests to develop the required database. 

The specific goals of the research program were to (1) investigate wing sensitivity to various 
forms of ice accretions, (2) compare the effects of actual ice shapes from icing tunnel tests with 
equivalent simulated ice shapes, (3) assess the effects of glaze ice shape features such as hom 
angle, hom height and surface roughness on wing aerodynamic performance, and (4) develop an 
experimental database of ice shape effects on the aerodynamic performance of a swept finite 
span wing with an aileron control surface. The principal accomplishments of this research 
program, which was completed in the fall of 2002, are summarized below. 

A number of meetings were conducted with the FAA, NASA, and general aviation 
manufacturers to select a wing model, ice shapes, and test conditions, and to define experimental 
methodologies and procedures. A 5 -ft semispan swept finite wing reflection plane model 
equipped with an aileron control surface was designed, fabricated, and instrumented at WSU. 
Tests were conducted at the NASA Glenn Icing Research Tunnel (IRT) to generate five glaze ice 
shapes with complete and incomplete scallop features and one rime ice shape. Ice shape castings 
were produced by NASA personnel from the actual ice accretions for aerodynamic testing. 
Simulated 3D ice shapes were defined with the NASA Glenn LEW1CE ice accretion code for the 
same icing conditions used in the IRT icing tests. The 3D ice shapes were fabricated out of 
wood or aluminum at WSU and were prepared for aerodynamic testing. Extensive experiments 
were conducted at the WSU 7- x 10-ft wind tunnel facility over a period of 10 weeks to generate 
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the required data for the clean and iced wing. Twenty ice shapes and eight roughness cases were 
tested along with the clean wing. The ice shapes included six IRT castings and seven smooth 
and seven rough LEW1CE shapes. The roughness cases included 120- and 150-grit sandpaper to 
simulate the effect of frost on wing aerodynamic performance. Lift, drag, pitching moment, 
hinge moment, and pressure distributions were obtained for all configurations tested. 

The experimental results obtained showed that the stall lift coefficients for the wing with the 
glaze ice shape castings were 1 1.5% to 93.6% less than the clean wing. For the 5-min rime ice 
shape, the stall lift coefficient was 3.4% higher than the clean wing. The IRT ice shape castings 
tested increased the minimum clean wing drag coefficient by 133% to 3533% and increased drag 
near stall by 17% to 104%. In general, the aileron remained effective in changing the lift of the 
clean and iced wings for all angles of attack and aileron deflections tested. Aileron hinge 
moments for the iced wing remained within the maximum and minimum limits defined by the 
clean wing hinge-moment data. In general, the trends in aerodynamic performance degradation 
of the wing with the simulated rough LEW1CE ice shapes were similar to those obtained with the 
IRT ice shape castings. However, in most cases, the ice shape castings resulted in greater 
aerodynamic performance losses than that obtained with the rough LEW1CE shapes. In most 
cases, the rough LEWICE ice shapes caused greater aerodynamic performance degradation than 
did their smooth counterparts. 
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1. INTRODUCTION . 


Aircraft certification requires the evaluation of the effects of ice accretions on aircraft 
aerodynamic components such as wings, tails, and control surfaces. Ice accumulation on 
aerodynamic surfaces can have a significant impact on aircraft performance, handling qualities, 
and thus, aircraft safety. A wide range of ice accretions is possible, depending on aircraft 
configuration, icing, and flow conditions. Potential ice accretions include glaze ice, rime ice, 
runback, and beak ice, as well as small ice shapes which can have considerable degradation in 
aircraft performance. In general, ice accumulation on aerodynamic surfaces can cause flow 
separation whose extent over the aerodynamic surface is a function of ice shape and wing 
geometry'. The term critical is often used to identify ice shapes responsible for large degradation 
in the aerodynamic performance of lifting surfaces. 

A number of experimental studies have been conducted over the years (refer to reference 1 for a 
comprehensive review) in an effort to assess the effect of various forms of ice accretions on 
aircraft aerodynamic performance and handling qualities. Most of these studies have been 
limited to two-dimensional (2D) airfoil sections due to tunnel and model cost constraints. Only a 
small number of investigations have addressed the impact of ice shapes on finite wuigs 

To address the lack of experimental data for three-dimensional (3D) iced wing configurations, 
scientists and engineers from the Federal Aviation Administration (FAA), the National 
Aeronautics and Space Administration (NASA), the general aviation aircraft industry, and 
Wichita Stale University (WSU) conducted a meeting to plan a tesearch program for a 
systematic evaluation of ice accretion effects on finite wings. The main objective of this 
research program was the development of a 3D experimental database of ice accretion effects on 
a swept finite wing. The database was needed for developing certification guidance material and 
for improving and validating simulation tools for aerodynamic analysis and design. A research 
grant was awarded to WSU to design and fabricate a wind tunnel model and to conduct icing and 
aerodynamic performance tests to develop the required database. The specific goals of the WSU 
research effort were to 

• investigate wing sensitivity to various forms of ice accretions. 

• compare the effects of actual ice shapes from icing tunnel tests with equivalent simulated 
ice shapes. 

• assess the effects of glaze ice shape features such as horn angle, horn height, and surface 
roughness on wing aerodynamic performance. 

• develop an experimental database of ice shape effects on aerodynamic performance of a 
sw ept finite wing with an aileron control surface. 

To accomplish the above objectives, wind tunnel tests were planned with a wing representative 
of modem business jet and regional jet aircraft wing planforms. The wing selection was based 
on input from the FAA, NASA, and the aircraft industry. An important consideration in the 
selection of the airfoil section for the finite wing was the availability of 2D iced airfoil 
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aerodynamic performance data for comparison with the 3D iced wing aerodynamic performance 
results. The airfoil selected was a GLC-305 section, which is representative of a modem 
business jet wing section. NASA has conducted extensive 2D wind tunnel tests with this airfoil 
in recent years using a range of ice accretions. 

The ice shapes selected for the experimental investigation included the following: 

• Po’ /urethane castings of ice accretions obtained from icing tests at the NASA Glenn 
Icing Research Tunnel (1RT). 

• Smooth LEW1CE ice shapes obtained for the same icing conditions as the ones used in 
the generation of the ice shape castings. 

• Rough LEW1CE ice shapes obtained by adoing grit roughness to the surface of the 
smooth LEW1CE ice shrpes. 

This report describes the icing tests performed in the 1RT and the aerodynamic tests 
conducted at WSU with the ice castings and the simulated LEWICE ice shapes. 
Experimental data presented includes lift, drag, pitching moment, aileron hinge moment, 
and surface pressure data. 
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2. EXPERIMENTAL FACILITIES AND PROCEDURES. 


The experimental i' 'esiigation consisted of ice accretion tests aad aerodynamic performance 
wind tunnel tests. The icing tests were conducted over a time period of 5 days at the IRT facility 
to obtain six ice shape castings for aerodynamic testing. The aerodynamic performance 
investigation took place at the 7- x 10-ft Low-Speed Wind Tunnel facility at WSU. The clean 
wing and 20 ice shape configurations were tested at the WSU wind tunnel over a time period of 
10 weeks. This section describes the test model and ice shapes tested, the experimental setup, 
test conditions, test measurements, and procedures used in the experimental investigations 
conducted. 


2.1 ICING TESTS VT THi- \ \SA G1.ENN IRT. 


The first part of the experiment was conducted at the NASA Glenn 6- x 9-ft Icing Research 
Tunnel to obtain castings of ice accretions formed on a swept wing model. Details of the icing 
tests are provided in the following sections. 

2.1.1 Test Facility. 

The IRT is a closed-loop refrigerated wind tunnel. Its test section is C ft (1.8 m) high, 9 ft 
(2.7 m) wide, and 20 ft (6.0 m) long. In the test section, the total air temperature can be varied 
between -20°F (-30°C) and +33°F (+1°C), within accuracy of ±1°F (±0.5°C), and a maximum 
velocity of 390 mph (160 m/s) can be attained. A spray system allows control of the liquid water 
content (LWC) between 022 to 3.0 g/m 3 and provides droplet median volumetric diameters 
(MVD) from 15 to 40 pm. Figure 2-1 shows the plan view of the NASA Glenn Icing Research 
Tunnel 



FIGURE 2-1. NASA GLENN ICING RESEARCH TUNNEL PLAN VIEW 
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2.1.2 Test Model. 


Aerodynamic consideration* and facility size limitations determined the overall size of the wing 
model. Details of the swept wing model are provided in figures 2-2 and 2-3. The model was a 
swept finite wing with a GLC-305 airfoil section aligned in the streamwise direction. The airfoil 
section had a maximum thickness to chord ratio of 0.087 at approximately 38% local chord. The 
airfoil section remained constant from the wing root to the wing tip. The wing had a 28° leading 
edge (LE) sweep, a 15.6° trailing-edge (TE) sweep, a 60-in. semispan, a 7.35 ft 2 area (half wing), 
an aspect ratio (AR) of 6.80 (left and right wing), a taper ratio of 0.4, and i. ' rometric twist of 0° 
at the root and -4° (washout) at the tip. The wing root and tip chords were 25.2 inches and 10.08 
inches respectively. In addition, the wing mean aerodynamic chord (MAC) was 18.72 inches 
and was located 25.74 inches from the wing root. The model was instrumented with 203 
pressure ports distributed chordwise at five spanwise locations corresponding to 15%, 30%, 
50%, 68%, and 85% semispan. 


AUlrOiL SECTION 

1. GLC-305 section in strtanwiM direction (l/c fnax =(X067 at 
jJcxQ .38). Airfoil section is constant from root to tip. 

2. Trailing edge thickness = 004 in (constant sps ra i sii) 





(X10 
(X06 
a oo 

-U06 

■ aw i ' i 1 i 1 i ' i 1 i 1 i 1 i 1 i 1 i 1 i 

0.00 0.10 0.20 0l30 0.40 CL5G 0.60 0,70 060 080 1430 

WING PLANFORM 

1. 5 It semspen (72% of WSU tunnel height) 

2. Roct chord = 25.2 in. Tip chord = 10.08 In, MAC = 16.72 in 

3. t/c vales linearly from root to tip 

4. PUnform Atm (half wing, S/2) « 7.35 II 2 ; wing arm m to tmml 
section area ratio * 7.357(7x10) = 0.105 

5w taper ratio = 0.4, trapezoidal planform, linear chord variation 
from root to tip 

6. AR « b z tS « 10* /(2*7.35) a 6.8027 

7. 25* quarter chord twaep, 28* leading edge sweep 

WING GEOMETRIC TWIST AND DIHEDRAL 

1. Geometric twist about 26% local chord (0* at root, -4* at tip) 

2. Dihedral « 0* AILERON CONTROL SURFACE 

1. 25% local chord (from 70% to 65% isniapa.i); 
Ctw* = 3.66 in, C«p * 2.712 in 

2. Elliptical leading edge (pivot at 60% local wing chord). 
Gap * 0.5% local chord 

3. Deflection range: -20* to +20* 

4. Aileron aree behind the hinge line, S. = 0 l 2647 ft 1 

6. AUaron average chord behind the tinge line, 
c. * [(2.1 63+2.82)72/1 23 - 021 16 ft 



FIGURE 2-2. WING PLANFORM FIGURE 2-3. AIRFOIL SECTION AND WING 

PARAMETERS 
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2.1.3 Test Model Hardware for the Icing Tests . 


The icing tests conducted for the purpose of generating ice castings for the swept finite wing 
required considerable pl annin g to maximize the number of test runs during each day of testing. 
In general, once an icing test was completed, the ice casting process required approximately 1 
day to provide a mold. Thus, the wing cannot be used until the mold was completed and 
removed from the wing. Anoincr limitation of the ice casting process was the spanwise length of 
the casting cannot exceed 25 inches. 

The goal of the icing tests conducted in the 1RT facility with the 60-in. semispan finite wing was 
to produce six ice castings approximately 68 inches long (the length of wing measured along the 
swept LE) in 5 days of testing. Due to the limitations of the ice casting process, only five 25-in. 
segments could be produced in 5 days, unless significant modifications were made to the wing 
model to allow multiple ice castings during each test day. This was accomplished by designing 
and constructing four wing leading edges. The first leading edge was for aerodynamic 
measurements only and extended the full span of the wing. Leading edges two through four, 
however, were divided into three spanwise segments each, as shown in figures 2-4 and 2-5. 
Thus, a total of nine leading-edge segments were produced in addition to the full-span leading 
edge. The removable leading edg^ were labeled top, middle, and bottom. The bortom 
removable leading edge (RLE) was 25 inches long and extended from the wing root to 36.79% 
semispan. The middle segment was also 25 inches long and extended from 36.79% semispan to 
73.78% semispan. The top RLE had a length of 17.95 inches and extended from 73.78% 
semispan to the wing tip. Note that the sum of the lengths of the three RLEs was 67.98 inches, 
which is the distance from the wing root to the wing tip measured along the swept leading edge. 
This distance is equal to the wing semispan divided by the cosine of the leading-edge sweep 

angle (i.e., — — — ). The advantage of having 1 8 leading-edge segments was that every time an 
cos 28° 

ice accretion test was completed, the leading-edge segment with the ice shape was removed for 
the ice moldin g process and a clean segment was placed on the wing. This permitted multiple 
icing runs to be conducted each test day. 

Additional model hardware had to be designed and fabricated to support the icing *ests at the 
NASA Glenn 1RT. The additional hardware included brackets and a 1-ft wing extension to 
permit the placement in the wing from 5.5 inches below the 1RT floor to 1 1.5 inches above the 
floor. This was required to place the wing inside the uniform LWC region of the icing cloud. 
By lowering the wing below the tunnel floor, LWC uniformity was maintained over the middle 
and top segments of the wing. With the wing in the high position, the middle and bottom 
segments of the wing leading edge were exposed to the uniform region of the icing cloud. For 
aerodynamic messure measurement, it was necessary to install the wing with its root section 
placed on the tunnel floor. This required additional hardware to be fabricated. The three wing 
placements with respect to the 1RT tunnel floor are shown in figures 2-6 to 2-8. 
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FIGURE 2-4, WING IN THE LOW FIGURE 2-5, SPECIFICATIONS OF THE 

POSITION SHOWING THE RLEs TOP, MIDDLE, AND BOTTOM RLEs 



a. Low wing position b. Normal wing position c. High wing position 

FIGURE 2-6. LOW, NORMAL, AND HIGH POSITIONS OF THE WING IN TP® IRT 
(Cuts wear© made to remove ice acaretion font wing at locations A. B, and C, / 
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WITH THE TOP LI 


Table 2-1 lists it * sk icing conditions selected to obtain the ice accretions for the aerodynamic 
investigation. These king conditions were selected to provide four glaze ice shapes with 
complete and incomplete scallop features, a 22,5-mln glaze ice accretion representative of an ice 
protection system (IPS) failure case, and a -jsin rime ice shape. Glaze ice accretions with 
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FOR TRACING PURPOSES 
(IRT-C522 at tip of bottom RL E) 




The process started when the LE was removed and taken to the cold room, The LE was placed 
inside a woodm box (figure 2-12) made specially for the melding process. The top of the b m 
mi me of its sides were removed to free up space foi handling the delicate ice accretion. Once 
the RLE was attached to the inside of die box with bolts, the removed side of the box was 
reinstalled. The box was left open at the top for pouring in the mold material. Figures 2-13 and 
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FIGURE 2-14. LEADING EDGE WITH ICE; 
ACCRETION PREPARED FOR MAKING 
THE MOLD— TOP VIEW 


; FIGURE 2** 13, MIDDLE RLE WITH IRT-i 
CS22 ICE ACCRETION ATTACHED TO 
THE INSIDE OF A WOODEN MOLDING; 
BOX— SIDE VIEW 







box coalmining the leading edge with the ice mcwtkm. The mold was left to cure overnight. In 
the mooring, the mold was separated from the lemiing edge. 

Once the molds were finished, they were mow to the model shop of the NASA Glenn Research 
Center, where the cartings were made. To make a easting, a tsold rv&§ placed k ‘U the same 
box. used to make it The wing RLE Hams was used to make the moid was placed in the same 
positww as when the mold was :tnade. This is illustrated in figure 2-15. A separation distance 
was left between the LE and the mold to create a thin wail for the casting, A polyu,rethane 
elastomer, CihaAMgy RP 6430, was poured in and left to cure. AH* curing, the LE was 
separated from the casting, and the mold material was removed from around the casting. 


FIGURE 2-15. THE MOLD AND LEADING EDGE WERE PLACED INSIDE THE 

BOX TO MAKE THE CASTING 





distance, The aileron surface behind the hinge line had a geometric mesa chord of 0,2 1 18 ft and 
« planfoim area of 0,264? ft 2 . 
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FIGURE 2-22. HOUSING FOR 
SERVOMOTOR, TORQUE CELL, AMD 
LVDT ^Photographed under the wing, near the 
trailing edge.) 


FIGURE 2-21. INSTALLATION OF 
CLEAN WING ON TUNNEL 
TURNTABLE— REAR VIEW 


Two Stage touminna (AS) brackets were designed to connect the wing to the tunnel external 
balance anti to raise the wing model 4,4 inches above the tunnel floor, Rasing the model above 
the tunnel floor was necessary to keep the wing root section above the floor boundary layer, 
which was approximately 1.5 inches high. Since the AJ brackets were exposed to the flow, the 
wing aerodynamic measurements would have been affected by the mmiymmc loads on the 
brackets. To eliminate this problem, a wooden fuselage-like body was used to shield the wing 
support brackets from the airflow. The streamlined body was 4.436 inches high and 55.6 inches 
long. Aerodynamic forces on the body did not affect the wing force and moment measurements 
i since it was not connected to the tunnel balance. Extensive computational flow dynamics (CFO); 
studies were performed prior to the experimental investigation to ensure that the streimlmetl 
body did net modify the wing flow field and that the flow over the body remained attached even 
ist high angles of attack (beyond wing stall). A small 0.25-tn. gap wan allowed between the wing 
root section a„d the stoeamlrted body to prevent contact between the wing and the body. The 
0,25-m, gap was sealed with a latex strip called dental dam, as shown in figures 2-23 and 2-24, to 
: prevent flow leakage through the gap flat could have modified the flow near the wing root. 
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Streamlined Body 




FIGURE 2-23. CLOSE-UP OF WING 
ROOT AND DENTAL 
DAM— PRESSURE SURFACE 


FIGURE 2-24. CLOSE-UP OF WING 
ROOT AND DENTAL 
DAM— SUCTION SURFACE 


2.23 Icing lese^ 


The ice shapes selected for the affodynwnic experiments included six polyurethane castings of 
actual ice accretions and seven simulated Ice shapes defined with the NASA Glenn LEWICE 2.0 
ice saretk® code. Icing conditions and ice shape notations are provided in tables 2-1 to 2-3. To 
assess foe effects of ice roughness on aerodynamic performance, the LEWICE ice shapes were 
tested with and without simulated roughness. A total of 20 ice shapes were investigated during 
foe exprirnentd study. 


2.2.3. 1 IRT Ice Shane Castings . 


The ice shapes toted were polyurethane castings of actual ice accretions obtained at the NASA 
Glenn IRT facility with foe GLC-305 wing model, as summarized in tables 2-1 and 2-2. For 
each icing condition, a set of three castinp were made from the top, middle, and bottom RLEs, 
and these castings were glued together using epoxy to make a foil span IRT ice shape casting. A 
total of six IRT ice castings were produced: IRT-CS10, IRT-IS10, IRT-SC5, IRT-CS2, IRT- 
CS22, and 1RT-IPSF22, corresponding to icing conditions I through 6, respectively. Figures 2- 
25 to 2-32 show foe installation of various IRT ice shape castings on foe wing model. Also, 
figure 2-33 presents close-up views of all six IRT ice shape castings. 




TABLE 2-3. ICING CONDITIONS FOR LEWICE ICE SHAPES 


Icing 

Cc®dMon 

Description 

AOA** 

(deg) 

V 

(mph) 

T»j*i 

(°F) 

LWC 

mvd: 

(pm) 

X 

(mm) 

1 

(Glaze)* 


4 

250.0 

; 25.0 

0.68 

20.0 

10.0 

HUH 

Incomplete Scallop Condition 
(ID: LS-IS10 and LR-ISI0) 

4 

150.0 

25.0 

0.65 

20.0 

10.0 

3 (Rime)* 

Scaled Condition from 2D tests 
(ID: LS-SC5 and LR-SCS) 

6 

201.3 

11.7 

0.51 

14.5 

5.0 

4 

(Glaze)* 


D 

250.0 

25.0 

0.68 

20.0 

2.0 

wlliiil 

mmmm 

n 

250.0 

25.0 ; 

0.68 j 

20.0 

22.5 

5 

(Glaze)* 

Complete Scallop Condition 
(ID: LS-CS22S and LR-CS22S) 

a 

250,0 



20.0 

22.5 

6 

(Glaze)* 

Failed Ice Protection 

(ID: LS-IPSF22 and LR-IPSF22) 

4 

150.0 

27.0 

0.46 

20.0 

22.5 


The methodology used to define the LEWICE ice shapes is discussed in section 2.23, 

Roughness was simulated with 36-stae loose gift (avenge roughness height is apfHroimately 0.021 J inch or 03336 
mm), 

’Computed using airfoil sections and t component of velocity normal to the wing tending edge (see section 223). 
^Ccmpitei using stieaxnwise velocity' and airfoil Factions (see mtim 223). 

**AQA = Angle of attack. 
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For ease of installation during testing, ice shape castings were made with cuffs tha* wrapped 
around the wing LE. The small ice shape castings were mounted to the wing using aluminum 
tape. The larger castings, IRT-CS10, IRT-CS22, and IRT-IPSF22, had aluminum brackets at 
selected spanwise locations that allowed these castings to be attached securely to the wing. This 
was necessary due to die considerable aerodynamic loads experienced by these large »ce shapes. 
The leading-edge cuffs and the brackets on the IRT castings provided the additional benefit of 
consistent ice shape attachment to the wing during repeated installations. 

2.2.3 .2 Smooth LEWiCE Ice Shapes . 

The simulated shapes were 3D ice shapes that were defined from a s .Ties of 2D ice sections 
obtained with the LEWICE ice accretion code [5]. Table 2-2 summarizes all icing conditions for 
LEWICE analyses. The procedure used to define the 3D LEWICE ice shapes is detailed in 
reference ( and is summarized below: 

1. Streamwise wing sections at 0% (wing root), 15%, 50%, 85%, and 100% (wing tip) 
semispan were selected for the development of the LEWICE ice shapes. 

2. Next, four additional wing sections were defined by taking the intersection of the wing 
with planes normal to the wing LE at 15%, 50%, 85%, and 100% semispan. The LE of 
the normal sections was at the same soanwisc location as the streamwise sections defined 
instep 1. 

3. 3D Navier-Stokes computations were performed at WSU with the clean wing. In the 
computations, the IRT walls were included to simulate the tunnel wall effects on the wing 
flow field. The geometric angles of attack (a) used in the analysis were 4° and 6° to 
match the angles of attack in the icing tests. Analysis pressures for streamwise sections 
at 15%, 50%, and 85% semispan were compared with experimental pressure distributions 
obtained in the IRT facility. Good correlation between experiment and analysis was 
demonstrated. From the computed flow fields, pressure distributions were obtained for 
the four wing sections normal to the wing LE defined in step 2 and for the streamwise 
section at the wing root. 

4. 2D ice accretion analyses were conducted with the LEWICE 2.0 computer code. The 
computations were performed using five sections of the GLC-305 swept finite wing. 
These included the streamwise section at the wing root and the four sections normal to 
the wing LE at 15%, 50%, 85%, and 100% semispan (see step 2 above). The icing 
conditions for the LEWICE analysis were identical to those used in the IRT icing tests 
(see tables 2-1 and 2-3). However, the LEWICE angle of attack for each wing section 
was adjusted to match the pressure distributions from the Navier-Stokes analysis 
described in step 3. The velocity for the ice accretion computations conducted with the 
four sections defined normal to the wing leading edge was set to equal the component of 
the ffee-stream normal to the wing LE. 

5. For each icing condition in table 2-3, five LEWICE ice shape sections (one streamwise 
and four normal to the wing leading edge) were obtained. The five ice shape sections 
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along with the wing geomehy were imported into a computer-aided design software 
package and were used to define 3D LEWJCE ice shapes. The 2D ice shape sections at 
the five spanwise locations were connecteo with spline surfaces and with plane surfaces. 
A comparison of the 3D ice shapes obtained using splice interpolation with 
corresponding ice shapes obtained with linear interpolation showed very small 
differences. Since the splme interpolation method provided a smooth transition in the 
spanwise direction between the five 2D sections, it was used to define the 3D LEWICE 
ice slopes, corresponding to icing conditions 1 to 4 (ice shapes LS-CS1U, LS-IS10, LS- 
SC5, and LS-CS2). Straight-line interpolation was used for the remaining two ice shapes 
(LS-CS22N and LS-1PSF22), which were the large 22.5-min glaze ice acc*etions. The 
reason for using straight-line interpolation for LS-CS22N nd LS-IPSF22 was to match 
the method used by some aircraft manufacturers in defining 3D LEWICE ice shapes from 
2D sections. 

The use of airfoil sections normal to the wing LE along with the normal component of the 
velocity vector for the LEWICE analysis was based on input from NASA, riowever, some 
airframe's prefer to use streamwise airfoil sections and the stream wise velocity in defining 3D 
ice shapes with the LEWICE computer code. Typically, streamwise ice accretion analyses 
produce ice shapes with larger horns (for glaze ice shapes) due to the higher water load resulting 
from the higher speed (streamwise velocity is greater than the velocity component normal to the 
LE). To compare the aerodynamic effects of LEWICE icc shapes based on normal and 
streamwise wing sections and flow velocities, one more (seventh) LEWICE ice shape was 
defined. This shape was obtained using the procedure discussed in steps 3 to o above. However, 
for this ice shape, streamwise wing sections at 0%, 15%, j0%, 85%, and 100% sero’spsn were 
used along with the streamwise free-stream speed to define five 2D LEWICE ice shares. The 
icing condition used for the seventh LEWICE ice shape (LS-CS22S) was icing condition ~ in 
table 2-1, which in the IRT icing tests produced a 22.5-min scalloped ice shape with Uuge horns. 

Close-up v’ews of all LEWICE ice shapes are shown in figure 2-34. The smooth LEWICE 
sections ant' . jrrcsponding IRT ice shape tracings are compared in figures 2-35 to 2-40. Section 
comparisons were made at spanwise locations labeled cut A, cut B, and cut C in figure 2-5. The 
sections at stations A and B wcie taken normal to the wing leading edge. The section at station 
C was in the streamwise direction. Section C was at the wing root, section B was 25 inches from 
the root (measured along the wing leading edge), and section A was at 50 inches from the root 
(measured along the wing leading edge). Note that the distance from the wing root to the wing 
tip measured along the wing leading edge was 67.95 inches. Figures 2-35 to 2-40 demonstrate 
that the icc accretions had a notable twist from the tip of the wing towards the root. This was 
due to the geometric twist of the wing. For all ice shapes obtained at a of 4°, the wing tip 
section was at abd of 0°. while the wing root was at +4°. Thus, the ice on the root accreted 
towards the lower surface of the wing; while at the wing tip, it was almost equally distributed on 
both surfaces. In comparing the LEWICF. sect’ons with the tracings of the IRT ice shapes, it is 
important to realize that the traces do not reflect the complex 3D features of the 1R1 icc shapes. 
A comparison between the 3D LEWICE and IRT ice casting shapes for icing condition 2 is 
provided in figure 2-41. Fig ures 2-42 and 2-43 compare the profiles of LEWICE icc shapes 
obtained for icing condition 5 using streamwise and normal sections to the wing LE as discussed 
above. 
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The effect of ice shape tm gkness was simalated fey adding 36-size (average roughness height m 
apprcaintafely 0,021 1 inches or 0,5356 mm) loose grit to the smooth LEWICE ice shapes. The 
rmighosed LEWICE (LR) ice shapes were named LR-CS10, LETS 10, LBL-5C5, LR-CS2, LR- 
CS22N, LR-C522S, and LR-BPSF22. Figures 2-44 to 2-51 show the installation of LEWICE ice 
shapes on the wing model 





FIGURE 2 - 47 . INSTALLATION OF LR- 
CS2 ROUGH LEWICE ICE SHAPE ON 
WING MODEL 


FIGURE 2-46, INSTALLATION OF LR- 
SC5 ROUGH LEWICE ICE SHAPE ON 
WING MODEL 
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A significant effort was directed at the start of the experimental irwestigation to ensure that all 
moasmmy hardware and software for the aenxlynamic experiments were in proper working order. 
Part of this effort Included a complete external balance calibration, testing of the pressure 
instiwnffltftkin, mi calibration of the aileron torque cell and LVDT. Next, a number of 
preliminary wind tunnel tests were conducted with the clean wing to evaluate test repeatability, 
assess tunnel flow angularity, and establish date-sampling rates and sampling period. Moreover, 
tests were jrerfbfmed to check force, moment, hinge moment, and pressure imtmm&mt&tim, and 
to evaluate the method for mounting the ice shapes onto the wing LE. During tire preliminary 
tests, the software for acquiring and reducing experimental data were tested for aecurecy using 
known static loads that were applied to the wing and the aileron. 

The third segment of the investigation was the production runs. The first set of production runs 
were conducted with the clean wing using fixed and free transition. Next, tests wore performed 
with all 20 ice shape- aMfipaaticms, which included six IRT ice shape castings, seven smooth 


i investigation included ADA sweeps for ten aileron deflections, m shown in table 24. The 
ileynolds number (Be) for all the tats was fixed at 1.8 million, and it was computed using tire 
wing MAC. For til runs, the complete set (6 ccmpcments) of forces and moments far tire whig: 
were obtained. Aileron hinge moments and wing surface pressures w«e also recorded. A 





limited Reynolds number study was also conducted with the clean and selected iced wing 
configurations for Reu\c of 0.5, 1.0, 1.5, and 1.8 million. A complete run log tor all the 
production runs is listed in appendix B. 


TABLE 2-4. TEST CONDITIONS FOR AERODYNAMIC PERFORMANCE DATA 


Reynolds number based on wing MAC (million) 

1.8 

Dynamic pressure, Q (psf) 

50 

Free-stream Mach number, M a 

0.185 

Free-stream velocity (V„) 

225.3 ft/s, 68.7 m/s, 153.7 mph, 133.5 kts 

Mach number normal to wing LE 

0.16 

Wing geometric a-range (deg.) 

-8 to +20 by 1° 

Aileron deflection, 5 a (deg.) 

-15,-10, -5, -2.5,0, 2.5,5, 10, 15,20 


Static tares were obtained for the clean wing and for all ice shapes investigated. There was no 
need to obtain dynamic tares for the wing model since the model mount was not exposed to the 
flow, and the streamline body fairing was not connected to the external balance. 

Flow angularity and wing downwash were measured for a portion of the runs. This was 
accomplished with two, seven-hole flow angularity probes, designed and calibrated by 
AEROPROBE. The seven-hole probes were installed on the clean wing at two stations 
corresponding to approximately 16% and 80% semispan. Seven-hole probes were selected 
because they can measure the three components of velocity, the total pressure, and the static 
pressure at a point in the flow. The seven-hole flow probes provided results with high accuracy 
for flow angles as high as 75° [7], The data from the seven-hole probes were used to estimate 
local a near the inboard and outboard sections of the wing. The tests with the flow probes 
included a-sweeps for three aileron deflections (0°, -15°, and +20°). The Reynolds number in all 
cases was 1.8 million, based on wing MAC. Installation of the seven-hole flow probes is 
depicted in figures 2-52 to 2-55. 

At the end of production runs, flow visualization tests were performed for selected 
configurations, which included all 1RT castings, two smooth LEW1CE ice shapes (LS-CS22S 
and LS-IPSF22), and all rough LEWICE ice shares. Visualization of the clean and iced wing 
flew fields was accomplished using white yam tufts attached to both surfaces of the wing model 
and the wing fairing. Each flow visualization run was performed at a Reynolds number of 
1 .5 million based on wing MAC and consisted of an a-sweep (-8° to +16°, increment of 1 °) with 
the aileron in the neutral position. Figures 2-56 and 2-57 show example installations of yam 
tufts on the clean and iced wing. Three video cameras were used to monitor the flow pattern on 
both surfaces of the model. One of the cameras, which was positioned outside the south wall of 
the tunnel, was used to capture the flow pattern on .he suction (upper) surface of the wing. The 
other two cameras were located above the window at the tunnel ceiling and outside the tunnel’s 
north wall. These cameras recorded the flow pattern on the pressure side of the wing. Separated 
flow patterns near the wing tip of an iced wing configuration are illustrated in figures 2-58 and 
2-59. 
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FIGURE 2-52., INSTALLATION OF 
SEVEN-HOLE FLOW PROBES ON 
WING MODEL-— FRONT VIEW 


FIGURE 2-53. INSTALLATION OF 
SEVEN-HOLE FLOW PROBES ON 
WING MODEL— SIDE VIEW 




FIGURE 2-56. TUFTS INSTALLATION FIGURE 2-57. INSTALLATION OF 
ON CLEAN WING AND STREAMLINED YARN TUFTS ON ICED WING 

BODY-PRESSURE SURFACE AND STREAMLINED 

BODY— SUCTION SURFACE 



FIGURE 2 58. TUFTS SHOW THAT 

FLOW DOWNSTREAM Of ICE SHAPE FIGURE 2-59. TUFTS SHOW COMPLETE 
WAS SEPARATED NEAR WING TIP FLOW SEPARATION AT WING TIP 
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2.2,5 Test Measmianents , 


Aerodynamic forces and moments were measured with a six-component truncated pyramid-type 
exterml balance. Aileron hinge moments were measured with a lO-k-lb temperature 
compmsated torque sensor. The resolution of the hinge-moment torque cel! was 0.005 in-lb. 
Force balance measurements were resolved as follows: ±0.75 lb for lift, ±03 lb for drag, and 
±0.2 ft-Ib for pitching moment. The resolution of the aerodynamic coefficients obtained at a 
tunnel dynamic pressure of 50 pf was 0,002 for lift coefficient (Cl), 0.0008 for drag coefficient 
(Co), 0,0003 for pitching-moment coefficient (Cm), and 0.0001 for hinge-moment coefficient 
(Cm). 

Pressure measurements were obtained wild a Pressure Systems Inc. 8400 Industrial System 
processor. Pressures were resolved to ±0,0025 psi corresponding to » pressure inefficient (Cp) 
of 0.0072 for a test dynamic pressure of 50 psf. 

Aerodynamic coefficients for each a were computed by taking the average of two sets of 
measwemeots. Each set consisted of 1024 data points (per channel) obtained over a time 
interval of 4 seconds. Thus, the value of each of the force and moment coefficients was 
computed by taking the average of 2048 data points per a. 
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Conditions for the IRT icing teste are summarized to table 2-1 and include five glaze conditions 
with ice accretion times (x) of 2, 10* and 22,5 minutes and a 5-min rime ice case. Ice accretion 
characteristics for each icing condition are discussed below. 



3,1,1 Icing Condition 1— Complete Scallop Condition (7= 10 min) . 


At this glaze icing condition, the ice accretion on the top RLE segment (figure 3-1) shows well- 
defined scallop tips beginning at the attachment lines. The space between the scallop tips shows 
irame accuBiiilMiott of ice. The scallop tips are covered with roughness elements. The side view 
of the ice accretion on the pressure side of the wing (lower surface), shown in figure 3-2, 
i demonstrates the feather formation of the scallop tips. Figure 3-3 provides an o verall view of the 
ice accretion obtained at the top RLE. The casting photograph in figure 3-4 shows how the ice 
i accretion increases in size away from the tip of die wing. The ice accretions on the middle and 
; bottom RLEs exhibited the same characteristics m the ones observed with the ice on the top 
RLE. 




FIGURE 3-2. ICING CONDITION 1, SID1 
VIEW OF ICE ACCRETION ON 
PRESSURE SIDE OF THE TOP RLE 


FIGURE 3-1. ICING CONDITION 1, 
FRONT VIEW OF ICE ACCRETION ON 
THE TOP ft | p 
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This condition was obtained by scaling a reference icing condition [9] tested witl» a 2D 36-in. 
chord GLC-305 wing model. The icing condition for the 2D wing was a = 6®, V* * 201 .3 mph, 
T««i:= H.7°F, LWC - 0.4O g/m J , MVD = 20 pm, and t = 16.7 minute. Aerodynamic k M with 
the 2D airfoil and the ice shape obtained with the reference icing condition were conducted in the 
Low-Turbulence Pressure Tunnel (LTPT) at NASA Langley Research Center [10). Comparison 
of the aerodynamic performance data from the LTPT tests with c»ri»p>nding results front the 
WSU finite wing tests using the scaled ice shape will help evaluate differences between 2D and 
3D iced wing mmdymmim. 
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FIGURE 3-8, ICING CONDITION 2* 
CASTING FROM THE ICE ACCRETION 
ON THE BOTTOM RLE 


FIGURE 3-7. ICING CONDITION 2, 
FRONT VIEW OF THE COMPLETE ICE; 
ACCRETION ON THE BOTTOM RLE 





FIGURE 3*41, ICING CONDITION 3. 
FRONT VIEW OF THE COMPLETE ICE 
ACCRETION ON THE BOTTOM RLE 


FIGURE 3-12. ICING CONDITION 3, 
CASTING FROM THE ICE ACCRETION 
ON THE BOTTOM RLE 





FIGURE 3-15. ICING CONDITION 4, FIGURE 3-16. ICING CONDITION 4, 
FRONT VIEW OF THE COMPLETE ICE CASTING FROM THE ICE ACCRETION 

ACCRETION ON THE MIDDI E RLE ON THE MIDDLE RLE 
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FIGURE 348. ICING CONDITION 5, SIDE 
VIEW OF ICE ACCRETION ON SUCTION 


FIGURE 347. ICING CONDITION 5* 
i FRONT VIEW OF ICE ACCRETION ON 


FIGURE 3-20. ICING CONDITION Si 
CASTING FROM THE ICE ACC1F 
ON THE MIDDLE RT £ 


FIGURE 349. ICING CONDITION 5* 

front view of the complete ice 

ACCRETION ON THE MIDDLE RLE 






When ttti icing condition was run with the wing model in the high position, ft piece of the ice 
acacias was asitfiiteallf shed off the whig. For tMi reason, the vdtodiy for icing awltttai 5 
was reduced to 225 mph with the whr* in the high position. To miintiiiii the mam wrier mass 
flax as in the case at 250 tnpb, the LWv .«§ ina»a»d to 0.89 gfm*. Photopiphle data from the 
ice accf edon at 250 tnpb and the one ©turned at 225 mpt* tiiicsted that the two foe acaet»r« 
were very close is overall ciaursctaristics. Tie ice wrings made from the foe amttfon cAtaiaed 
»i 225 tnph and 0.89 ffrn 1 were the middle and totem RUEs. The foe casting for the top ELI 
segment was obtained flan the foe atxreifoct fats conducted at a speed of 250 mpfa and LWC of 
0.M gfrn s . 












where Hui is the hinge-moment and 5, and c a are the aileron planform area and the aileron mean 
chord behind the hinge line respectively. (The values for S a and c a are provided in figure 2-3.) 
The sign convention for C H is positive for hinge moments causing the aileron trailing edge to 
move down. Aileron deflections corresponding to trailing edge down are also considered 
positive. Figure 3-25 depicts the sign convention for Ch and 5 a. 

In the following discussion, the hinge-moment curve is divided into three regions (A, B, and C), 
corresponding to the linear, near-stall, and poststall ranges of the lift curve, as shown in 
figure 3-26. Typically, in region A, the decrease of Ch was linear, small and gradual, while in 
region B, the growth was much greater due to increasing flow separation over the control 
surface. Note the anomaly or break from region A to region B with increase Cl and the 
differences between the clean and iced configurations. In region C, the hinge moment remained 
nearly constant, since beyond stall the pressure distribution o.er the aileron upper surface did not 
vary significantly for positive a. 



I 


| 



(a) Clean configuration. Cl > 0 



t 


(b) Clean configuration. Cl < 0 



(c) Iced configuration. Cl > 0 


(d) Iced configuration. Cl < 0 


FIGURE 3-26. REGIONS A, B, AND C FOR HINGE-MOMENT CURVE OF CLEAN AND 

ICED CONFIGURATIONS 
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All pressure data are presented in the form of Cp versus chordwise distance (x/c). In all cases, 
negative Cp coefficients corresponding to the suction side (upper surface) of the wing are plotted 
upward. 

3.2.1 Data Quality and Repeatability . 

The WSU wind tunnel facility has been used for commercial and research testing over the years. 
Facility personnel have considerable expertise in wind tunnel testing, instrumentation, data 
acquisition, and data processing. To obtain quality' data, careful planning and well thought out 
test were implt-m^nted during the aerodynamic investigation as discussed below. 

• An analysis was performed to determine instrumentation sensitivity needed to provide the 
required resolution of force, moment, and pressure coefficients for the primary test 
condition, which was a-sweeps at Reynolds number based on MAC of 1.8 million. 

• Prior to WSU wind tunnel entry, a complete external balance calibration was performed. 
Instrumentation such as pressure transducers, torque tubes, load cells, flow probes, and 
model hardware were checked to ensure proper operation and accuracy. 

• The repeatability of the experimental data is a function of instruments don, daia-sampling 
rates, tlow quality', flow unsteadiness, and model setup procedures. To verify the 
repeatability of the experimental setup, tests with the baseline configuration were 
repeated, and the results were compared. Repeatability of lift, drag, pitching-moment, 
and hinge-moment coefficients are provided in figure 3-27. Also, figures 3-28 to 3-30 
show pressure distributions of 15%, 50%, and 85% semispan locations at a of -4°, 0°, 4°, 
8°, 12°, and 16°. The maximum average percentage difference of aerodynamic 
coefficients from individual test runs from the average of all test repeats was 1% for CL, 
2% for Cd, 1% for Ch, and 3% for Cp. Most of the variation occurring near wing stall 

• At regular intervals during the wind tunnel tests, the balance, the hinge-moment system, 
and the pressure transducers were tested with known inputs to verify that they were 
working properly. For the balance tests, known forces were applied to the model, and 
force and moment data were obtained for a complete a -sweep at zero airspeed. 
Figures 3-31 and 3-32 show the setup for balance and hinge-moment system checking. 

• To ensure repeatability in the installation of the ice shape tested, all ice shapes were 
attached to cuffs designed to fit the wing LE. For the large ice shapes, brackets were also 
installed at selected span wise locations to provide additional support and minimize 
deflections due to the aerodynamic loads. 

• Preliminary tests were conducted with all models prior to the start of the production runs 
to verify tunnel and model instrumentation and data acquisition hardware and software. 
Static tares were obtained for all ice :ihapes to account for the weight effects of each ice 
shape on the aerodynamic measurements. 
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• A limited number of tests were conducted prior to the production runs to investigate data- 
sam p lin g periods, as shown in figure 3-33. All the force and moment coefficients 
obtained from wind tunnel tests are average values based on multiple measurements over 
a period of time. Sampling rate and sampling time depend on tunnel facility, test model, 
and airspeed, and are usually established experimentally for each tunnel facility. For the 
WSU tests, 1024 readings were taken per a for each force and moment coefficient. The 
sampling time was 4 seconds. For a dynamic pressure of 50 psf, the airspeed in the WSU 
wind tunnel was approximately 222 ft/s. Given that the MAC of the wing is 1.56 ft, 1 
second of data acquisition was equivalent to averaging the flow over the model 142 
times, in 1 second, the air had traversed a distance equal to 142 chord lengths. 


• All tests were conducted at constant Reynolds number, i.e., the speed of the tunnel was 
adjusted to maintain constant Re^AC- 




(a) Lift coefficient 



FIGURE 3-27. REPEATABILITY OF C L , C* C M , 



-10 •* 5 10 15 20 25 

Aagto of attack (dag) 


(d) Hinge-moment coefficient 
AND C H ; CLEAN CONFIGURATION; 


Re= 1.8xl0 6 ; 5 A = 0° 


3-11 






Ptmsufi coefficient pressure coefficient Pressure eosfflder* 
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(a) From view of balance check setup 


(b) Rear view of balance check setup 


(b) Side view of hinge- moment load cell check setup 




(a) Lift coefficient 


(b) Orag coefficient 



(c) Pitching moment coefficient 


(d) Hinge-moment coefficient 


FIGURE 3-33. DATA CAMPLING STUDIES OF C L , Co, C m , AND C H ; CLEAN 
CONFIGURATION; Re= 1.8xl0 6 ; 6 A = 0° 


In most 'and tunnel studies, aerodynamic performance measurements rely on standard wind 
tunnel instrumentation, which typically include external or internal balances and, in some cases, 
pressure instrumentation. Near stall, extensive flow separation and vortex shedding increase 
model vibration and cause large wake unsteadiness which can affect the accuracy of the 
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measured aerodynamic coefficients. Considering the lift coefficient, the objective in most wind 
tunnel tests is to determine lift curve slope, near-stall and poststall behavior, and C for 
^nginpffrin^ purposes. In determining the degree of measurement accuracy needed, an analysis 
should be performed based on the acceptable level of accuracy with respect to the aircraft 
performance. For example, assuming an error in Cl of 0.0 i for the wing used in this study, the 
corresponding error in lift at a speed of 131 kts (151 mph, Q = 50 psf) was approximately 3.7 lb. 
Given a Cua»x of 0.87, the total lift generated by the wing at the same flow condition was 
3 k' 7 lb. For engineering purposes, error of 3.7 lb in lift out of a total lift of 319.7 lb should 
be considered to determine if better accuracy is needed in measuring Cumi. 

Precise measurement of aerodynamic properties near stall is a very difficult task, which is 
beyond the scope of the study described in this report. The experimental difficulties stem from a 
number of factors such as the ones listed below: 

• The exact dynamic responses of the balance with the model installed (e.g., response 
versus frequency) should be known and appropriate daia-sampling rates and sampling 
periods should be established from experiments. The sampling periods should be long 
enough to account for both model and tunnel flow unsteadiness, 

• Model deflections and variations in a due to model vibration and flow pulsing are 
difficult to accurately measure or monitor in real time. 

• Tunnel blockage effects under stall conditions can be very significant, depending on 
model and tunc - 1 .ize, and their impact on tunnel dynamic and static pressures are 
difficult to determine and correct 

Generally, experimental data for unsteady conditions require special tunnel facilities and are 
usually obtained with simple geometries such as airfoil sections. In most cases, the best way to 
obtain such data is through the use of extensive surface pressure time histories, which can then 
be integrated to provide the required coefficients. 

In summary, the experimental data provided in this report have been obtained under carefully 
controlled conditions. Data near stall should be used with the understanding that it may be 
subject to some uncertainty. 

3.2.2 Compressibility Effects . 

For the clean GLC-305 airfoil section, the critical free-stream Mach number for 2D flow was 
determined to be 0.22. This value was obtained from "He intersection of the two curves defined 
by equations 3-2 and 3-3 [ to - • 1 14], 


cp* = 


r-M; 


r + 1 




7-1 


-1 


(3-2) 
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Cp = 


(3-3) 
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CPo 
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Cp, 


In the equations above, Cpa- is the critical pressure coefficient for which the local Mach number 
is 1, Cpa is the maximum suction press ore coefficient for incompressible flow, and Cp is the 
compressible pressure coefficient obtained using Cpo and Lai tone’s compressibility correction 
given by equation 3-3. The value of Cpa for the GLC-305 section was -10.1 and corresponded to 
a of 13 5°. This value was obtained from a 20 incompressible viscous flow analysis using the 
XFOIL [15] computer code. 


Note that the critical Mach number of 0.22 obtained from the 2D analysis and equations 3-2 and 
3-3 is a conservative estimate. Typically, in 3D flow, the value of the free-stream Mach number 


(Af*) required to achieve sonic flow over the wing for fixed a will be higher due to wing sweep 
and 3D flow relief effects. For the swept wing tested, the critical Mach number lor 3D flow was 

f t\ n A 


0249 


0.22 
cos 28° 


Since all tests at the WSU wind tunnel facility were conducted at of 0.185, 


rtf rrtmprrcsihility rtn arTrvtynumir raftffir.ipnts. uim nnl ^ignifirjinr 


3.2.3 Clean and Iced Wing Performance . 

Aerodynamic performance for the clean and iced wing is presented in figures 3-34 to 3-45 and in 
tables 3-1 to 3-9 for 5 a of 0°. The data presented demonstrates the effects of all the ice shapes 
(six 1RT ice shapes, seven smooth LEW1CE ice shapes, and seven rough LEWICE ; ^e shapes) 
investigated on lift, drag, pitching moments, hinge moments, and on the surface pressing 
distributions. For the purpose of discussion, the percentages in tables 3-1, 3-2, 3-4, 3-5, 3-7, and 
3-8 have been calculated using the following formula: 


Change in property X (%) = AX = Ussi — xlOQ (3-4) 

X cfeaa 

and, increase in X means that a negative X becomes more negative and a positive X becomes 
more positive. 
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(c) Drag coefficient (d) ACo vs a 



(e) Pitching -moment coefficient 


(f) Hinge-moment coefficient 


FIGURE 3-34. EFFECT OF 1RT ICE SHAPES ON C L , C D , C M , AND C H ; IRT-CS10* 
IRT-IS10, IRT-SC5, IRT-CS2, IRT-OS22, AND IRT-IPSF22 CONFIGURATIONS; 

Re= 1.8x10 6 ;6 a = 0° 
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Drag coafflrtarft 















(e) Cp vs x/c (a • 12°) 


(f) Cp vsx/c (a - 16°) 


FIGURE 3-39. EFFECT OF SMOOTH LEWICE ICE SHAPES ON PRESSURE 
DISTRIBUTIONS AT 15% SEMISPAN; LS-CS10, LS-1S10, LS-SC5, LS-CS2, LS-CS22N, 
LS-CS22S, AND LS-IPSF22 CONFIGURATIONS; RE= 1.8X10 6 ; S A = 0° 
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FIGURE 3-42. EFFECT OF ROUGH LEWICE ICE SHAPES ON C L , C D , C M , AND C u ; 
LR-CS10, LR-1S10, LR-SC5, 1.R-CS2, LR-CS22N, LR-CS22S, AND LR-IPSF22 
CONFIGURATIONS; Re= 1.8xI0 5 ; 6 A = 0° 
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TABLE 3-1. EFFECT OF 1RT ICE SHAPES ON C L AND Cd, Retuc= 1.8xl0 6 ; S A = 0° 


a (deg) 

Configuration 

C{_ 

AC L 

Cd 

AC d 


CLEAN 


0.0% 

0.008 

0.0% 


1RT-CS10 

0.19 

-13.6% 

0.074 

825.0% 


IRT-ISlo 

0.20 

-9.1% 

0.050 

525.0% 

3 

IRT-SC5 

0.22 

0.0% 

0.016 

100.0% 


irt-cs: 

0.21 

-4.5% 

0.021 

162.5% 


IRT-CS22 

0.05 

-77.3% 

0.224 

27u0.0% 


LRT-IPSF22 

0.18 

-18.7% 

0.081 

- 

912.5% 


CLEAN 

0.38 

0.0% 

0.01' 

0.0% 


IRT-CS10 

0.35 

-7.9% 

0.089 

535.7% 


IRT-ISIO 

0.37 

-2.6% 

0.064 

357.1% 

5 

1RT-SC5 

0.38 

0.0% 

0.023 

643% 


IRT-CS2 

0.37 

-2.6% 

0.029 

107.1% 



0.05 

-86.8% 

022° 

1535.7% 

4 1RT-IPSF22 

0.34 

-10.5% 

0.097 

592.9% 

"LEAN 

0.76 

0.0% 

0.053 



IRT-CS’.O 

0.54 

-28.9% 

0.165 

211.3% 


IRT-ISlU 

> . \.r"* 

i -15.8% 

0.154 

190 5% 

10 

IRT-SC5 

0.74 

-2.6% 

0.084 

583% 


IRT-CS2 

0.72 ! 

-5.3% 

0.100 

88.7% 


IRT-CS22 

0.21 

-72.4% 

0.261 

392 5% 


IRT-IPSF22 

0.5! 

-30.3% 


2243% 


CLEAN 

0.86 

0.0% 

0.198 

0.0% 


IRT-CSIO 

0.6G 

-30.2% 

02138 

2021% 


IRT-ISIO 

0.61 ! 

-29.1% 

02124 

13.1% 

15 

IRT-SC5 

0% 

4.7% 

C^21 

1 1.6/o 


IRT-CS2 

0.76 1 

-11.6% 

02109 

5.6% 



0.41 

i ••2.3% 

0.323 

63.1% 


IRT-IPSF22 

i 0.56 

-34.9% 1 

02134 

18.3% 


TABLE 3-2. EFFECT OF 1RT ICE SHAPES ON Q^u, <W. AND 
Cd^,; Rbkiac- 1.8xl0 6 ; & = 0° 


Config'^ralion 

Cl^uJ j 



Q. at 
a = 13. R 3 

ACLat 

a = 13.S’ 


ACdjtuh 

Q) at 

a = 13.S 3 

ACd { 

a =13.8° 1 

Clean 

0.87 j C.0% 

13.8° 



0.0% 

0.006 

0.0% 

0.147 

0.0% 

LkT-CSIo 



103° 

-3.3° 

0.56 

-35 6% 

0.072 

1100.0% 

0.209 

42.2% 

un-isio 

O.o4 

-26.4% 

10 .6° 

o 

r O 

i 

0.59 

-32.2% 

C.047 

683.3% | 0.198 

34.' 7 % 

IRT-SC5 

0.90 

3.4% 

15.8° 

2 ~ 

v,,86 

-1.1% 

0.0 \* 


17.0% 

IR C CS2 r 11 

-113% 

. - J 



0.76 

-12.6% 

0.018 

200.0% 

M74 

18.4% 

\ 1RT-CS22 

0.06 

-93.6% 

A0° 

S3 

0.36 

-58.6% 

02>18 

3533.3% 

0.300 

104.1% 

1R7-IPSF22 ! 0 5> -j9.1 10.5° 

-3 3’ 

033 

-?9.1% 

0.0*8 

1200.0% 

ii i 

0C10 

42.9% 
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TABLE 3-3. EFFECT OF IRT ICE SHAPES ON HINGE-MOMENT COEFFICIENT SLOPE; 

R(!uac= 1.8x10 6 ;8a = 0° 


Ice shape 

Region A 

Region B 

Region C 

dCH/da 
(per deg) 

a-f8nge 

(deg) 

dCH/da 
(per deg) 

a -range 
(deg) 

dCH/da 
(per deg) 

a-range 

(deg) 

Clean 

-0.0026 

OtolO 

-0.0155 

10 to 14 

-0.0001 

14 to 20 

IRT-CS10 

-0.0038 

0 to 6 

-0.0091 

6 to 10 

-0.0019 

10 to 20 



0to6 

-0.0106 

6to 10 

-0.0016 

10 to 20 

IRT-SC5 

-0.0022 

0 to9 

-0.0143 

9 to 15 

0.0004 

15 to 20 

IRT-CS2 

-0.0023 

0 to9 

-0.0137 

9 to 13 

-0 0016 

13 to 20 

IRT-CS22 

-0.0057 j 

0 to6 

-0.0065 

6to 11 

-0.0021 

11 to 20 

IRT-IPSF22 1 

-0.0041 

0 to6 

-0.0100 

6 to 10 

-0.0017 

10 to 20 


TABLE 3-4. EFFECT OF SMOOTH LEWICE ICE SHAPES ON C L AND Co, 

Reu 4c=1.8x10 6 ;6 a = 0° 


a (deg) 

Configuration 

C L 

AC l 

C D 

ACd 

3 

CLEAN 

0.22 

0.0% 

0.008 

0.0% 

LS-CS10 

0.18 

-18.2% 

0.079 

887.5% 

LS-IS10 


0.0% 

0.040 

400.0% 

LS-SC5 

H5EE955 

0.0% 

0.012 

50.0% 

LS-CS2 

0.22 

0.0% 

0.020 

150.0% 

LS-CS22N 

0.16 

-27.3% 

0.135 

1587.5% 

LS-CS22S 

0.12 

-45.5% 

0.155 

1837.5% 

LS-IPSF22 

0.22 

0.0% 

0.071 

787.5% 

5 

CLEAN 

0.38 

0.0% 

0.014 

0.0% 

LS-CS10 

0.34 

-10.5% 

0.093 

564.3% 

LS-IS10 

0.39 

2.6% 

0.058 

314.3% 

LS-SC5 

0.39 

2.6% 

0.017 

21.4% 

LS-CS2 

0.39 

2.6% 

0.030 

114.3% 

LS-CS22N 

0.25 


0.156 

1014.3% 

LS-CS22S 

0.21 

-44.7% 

0.168 

1100.0% 

LS-IPSF22 

0.36 

-5.3% 

0.088 

528.5% 

10 

CLEAN 

0.76 

0.0% 

0.053 

0.0% 

LS-CS10 

0.61 

-19.7% 

0.171 

222.6% 

LS-IS10 

0.63 

-17.1% 


* 79.2% 

L8-SC5 

0.79 

3.9% 

535 


LS-CS2 

0.6y 

-9.2% 

0.120 

126.4% 

LS-CS22N 


-43.4% 

0.233 

339.6% 

LS-CS22S 

0.43 

-43.4% 

02.38 

349.1% 

■ " ll "1 1 1 — 

0.61 1 

-19.7% 

0.169 

218.9% 

15 

CLEAN 

0.86 

0.0% 

0.198 

0.u% 

LS-CS.10 

07] 1 

-17.4% 

0.264 

33.3% 

LS-IS 0 

0.&7 

-22.1% 

0.231 

16.7% 

LS-SC* 

0.93 

0.1% 


-4.5% 

LS-CS2 

0.74 

-14.0% 

0.217 1 

9.6% 

LS-CS22N 

0.50 

-41.9% 

0.313 

58.1% 

i 

i LS-CS.2S 

0.55 

-36.0% 

0.319 

61.1% 

LS-IPSF22 

0.68 

-20.9% 

0.1*56 

29.3% 














































































































































































































TABLE 3-5. EFFECT OF SMOOTH LEW1CE ICE SHAPES ON C . cw, AND Ci w 

ReuAC^ 1.8x10 {, ;5 a = 0° 


Configuration 


ACu*: 

CWl 

Aa,ai: 

Cl at 
a=13.8° 

AC L at 
0=13.8° 


AC D m ir. 

Coal 

a-13.8° 

6Cq ax 
a- 13.8° 

Clean 

0.87 

0 0% 

13.8° 

0.0% 

0.87 

0.0% 

0.006 

0.0% 

0.147 

0.0% 

LS-CS10 

0.71 

-18.4% 

14.6° 

5.8% 

0.70 

-19.5% 

0.078 

1200.0% 

0.229 

55.8% 

I -S-IS J P 

0.6 

-21.8% 

13.6° 

-1.4% 


-21.8% 

0.036 

500.0% 

0.202 

37.4%, 

LS-SC5 


14.8° 

12 % 

0.93 

6.9% 

0.010 

66.7%) 

0.124 

-15.6% 

LS-CS2 

0.75 

-13.8% 

14.7° 

6.5% 

0.74 

-14.9% 

0.017 

1833% 

0.182 

23.8% 

LS-CS22N 

0.49 

■ 43 . 1 % 

13.5° 

-2.2% 

0.49 

-43.7% 

0.128 

20333% 

0389 

%.6% 

LS-CS22S 

0.54 

-37.9% 

13.4° 

-2.9% 

0.54 

-37.9%, 

0.148 


0.294 

100.0% 


0.69 

-20.7% 

14.6° 

5.8% 

0.68 

-21.8%o 

0.066 

1000.0% 

0.225 



TABLE 3-6. EFFECT OF SMOOTH LEWICE ICE SHAPES ON HINGE-MCM. \ . 
COEFFICIENT SLOPE; Rbuac^ L8x10 6 ; 6 a « 0° 


Ice Shape 

Region A 

Region B 

Region C | 

dC»da 
(per deg) 

a-range 

(deg) 

dQyda 
(per deg) 

a-range 

(deg) 

dCu/da 
(per deg) 

a-range 

(leg) 

Clean 

-0.0026 

0 to 10 

-0.0155 

10 to 14 

-0.0001 

14 to 20 

LS-CS10 

-0.0134 

0 to 7 

-0.0014 

7 to 15 

-0.0024 

14 to 20 

LS-IS10 

-0.0067 

0 to 7 

-0.0033 

7 to 13 

-0 0023 

13 to 20 

LS-SC5 

-0.0027 

Ota 11 

-0.0175 

11 to 14 1 

.'.•5 

14 to 20 

LS-CS2 

-0.0031 

0to8 

-0.0068 

8 to 14 

-S.0OI7 

14 to 20 

LS-CS22N 

-0.0089 

0 to 6 

-0.0039 

6 to 14 

-0.0060 

14 to 20 

LS-CS22S 



-0.0045 

7 to 13 

-0.0016 

13 to 20 

LS-IPSF22 

-0.0109 

0 to 7 

-0.0014 j 

7 to 14 

-0.0021 

14 to 20 
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TABLE 3-7. EFFECT OF ROUGH LEW1CE ICE SHAPES ON C L AND Coi 

Remc =1.8x10''; 8 A = C° 


a (deg) 


ConliguraUoo 


Cl 


AC l 


Cd 


AC d 


CLEAN 


LR-CS10 


LR-IS10 


LR-SC5 


LR-CS2 


LR-CS22N 


LR-CS22S 


LR-IPSF22 


0.22 


0 . 0 % 


0.18 


0.22 


0.21 


0.21 


0.17 


0.15 


0.21 


-18.2% 


0 . 0 % 


-4.5% 


-4.5% 


-22.7% 


-31.8% 


-4.5% 


0.008 


0.077 


0.043 


0.0:4 


0.023 


0.189 


0.214 


0.070 


0 . 0 % 

*2.5%' 


437.5% 


75.0“ 


187.5% 


2262.5% 


775.0* 


CLEAN 


LR-CS10 


LR-IS10 


LR-SC5 


LR-CS2 


LR-CS22N 


LR-CS22S 


LR-IPSF22 


0.38 


0 . 0 % 


0.33 


-132% 


0.38 


0 . 0 % 


0.38 


0 . 0 % 


0.38 


0 . 0 % 


0.24 


-36.8% 


0.21 


-44.7% 


0.36 


-5.3% 


0.014 


0.091 


0.060 


0.019 


0.032 


0.202 


0.225 


0.089 


0 . 0 % 


550.0% 


328.6% 


35.7% 


128.6% 


1342.9% 


1507.1% 


535.7% 


10 


CLEAN 


LR-CS10 


LR-IS10 


LR-SC5 


LR 2S2 


LR-CS22N 


LR-CS22S 


LR-IPSF22 


0.76 


0 . 0 % 


0.55 


-27.6% 


0.61 


-19.7% 


0.76 


0 . 0 % 


0.69 


-9.2% 


0.40 


-47.4% 


0.34 


-55.3% 


0.57 


-25.0% 


0.053 


0.168 


0.147 


0.058 


0.121 


0254 


0266 


0.167 


0 . 0 % 


217.0% 


177.4% 


9.4% 


128.3% 


3792% 


401.9% 


215.1° 


15 


CLEAN 


LR-CS10 


LR-IS10 


LR-SC5 


LR-CS2 


LR-CS22N 


LR-CS22S 


LR-IPSF22 


0.86 


0 . 0 % 


0.61 


-29.1% 


0.63 


-26.7% 


0.98 


14.0% 


0.72 


-16.3% 


0.41 


-52.3% 


0.35 


-59.3% 


0.61 


-29.1% 


0.198 


0249 


0226 


0212 

0218 “ 


0.307 


0.311 


0246 


0 . 0 % 


25.8% 


1«.H 


7.1% 


10 . 1 % 


55.1% 


242% 


TABLE 3-8. ETFECT OF ROUGH LEWICE ICE SHAPES ON C ■»», 
(Xcuii» AND Q^; Re^AAc = LSxlO 6 , 6* = 0° 


ConiiRuralion 

I9RR1 



IMS® 

1 C L at 

.. - P H° 

ACl at 

a = 13 8° 1 


m 

Cd at 
a - 1?C° 

AC D at 

a - 13.3° 

Clean 

0.87 

0.0% 

13.8° 

O 

o 

o 

0.87 ! 

0.0% 

0.006 1 

0.0% 

0.147 

0.0% 

LR-CS10 

0.61 

-29.9% 

13.6° 

i 

o 

to 

o 



0.076 1 

1 1 66.7% 

0.221 

50 3% 

LR-IS10 

0.65 

-25.3% 

13.6° 

- 0 . 2 ° 

| 0,64 

-26.4% 

0.C39 

550.0% 

0200 

36.1% 

LR-SC5 

0.98 

12.6% 

15. 8** 


MEEH 

5.7% 

0.013 


0 146 

-0.7% 

LR-CS2 

0.72 

BUI 

13.7° 

-0.1° 

022 

-172% 

0,020 

mmm 


252% 

LR-CS22N 

0.41 

-52.9% 

11.4° 

-2.4° 

C.40 

-54.0% 

c.’si ! 

— - - -- J 

2933.3% 

0.291 

98.0% 

LR-CS22S 

0.34 

-60.9% 

io.:° 

-3.5' 

0.33 

-62.1% 

r 0.205 

3316.7% 

0.295 

100.7% 

LR-1PSF22 

C.ol | 

-29.9% 

n.e 

L^_ 

0.61 

-29.9% 1 

0.065 1 

983.3% | 

| 0.2 IS 

48.3% 
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TABLE 3-9. EFFECT OF ROUGH LEWICE ICE SHAPES ON HINGE-MOMENT 
COEFFICIENT SLOPE; ReuAC= 1.8xl0 6 , 6 A = 0° 



Region A 

Region B 

Region C | 







Clean 

-0.0026 

0to 10 

-0.0155 

10 to 14 

-0.0001 

14 to 20 

LR-CS10 

-0.0164 

0 to 6 

-0.0017 

6to 13 

-0.0023 

13 xo 20 

LR-IS10 

-0.0062 

0 to 6 

-0.0040 

6 to 13 

-0.0021 

13 to 20 

LR-SC5 

-0.0026 

Gto 11 

-0.0157 

11 to 15 

0.0005 

15 to 20 

LR-CS2 

-0.0032 

0 to 8 

-0.0101 

8to 12 

-0.0016 

12 to 20 

LR-CS22N 

0.0007 

0 to 7 

-0.0036 

7 to 1 1 

-0,0100 

11 to 20 

LR-CS22S 

0.0007 

0 to 5 

0.0009 

5 to 10 

-0.0120 

10 to 20 

LR-IPSF22 

-0.0160 

0 to 6 

-0.0013 

6to 12 

-0.0025 

12 to 20 


3.2,3. 1 Clean Wing . 

Maximum lift coefficient for the clean wing was 0.87 and occurred at an a of 13.8°, as shown in 
figure 3-34(a) and in table 3-2. The lift dropped gradually after stall and wa» reduced to about 
0.81 at a of 20°. The slope of the linear portion of the lift curve was 4.41 per radian or 0.077 per 
degree- This slope correlates well with the slope of 0.064 obtained from equation 3-5 [16]. 



■ AR 



(— 1 -(Afi'Wj 2 


(3-5) 


where ao is the 2D lift slope of the GLC-305 airfoil and is equal to 0.084 per degree based on the 
data of reference 17, K\a is the wing sweep angle at 50% chord and is equal to 22°, AR is the 
wing aspect ratio (6.8) given in figure 2-3 and AL is the ffee-stream Mach number, which fb; 
the WSU wind tunnd tests was 0.185, as shown in table 2-4. 

From examination rf the pressure distributions and the flow visualization data obtained, flow 
separation was initiated at the wing LE and was combined with a leading-edge vortex. Flow 
separation was first observed for a of 9° near the 30% semispan station. The region of flow 
separation was very small in both the span wise and chord wise directions (2% chord and about 
3% semispan). As a was increased to 10°, a large trapezoidal region of rough and separated 
flow was observed between the 30% and 85% cemispan stations and between 12% and 60% of 
wing chord. At a of 1 1 °, considerable TE separation occurred near the wing trailing edge over a 
chord length of about 20%. Trailing edge separation extended from about 55% semispan to the 
wing tip. At a of 12°, flow separation was observed over the wing tip (outboard 15% of 
semispan) and near the TE of the inboard (5% to 30% semispan) portion of the wing. Finally, at 
a between 13° and 14°, complete flow separation occurred over most of the wing upper surface 
which led to wing stall. 

Drag performance for the clean wing is presented in figure 3-34(c) and in tables 3-1 and 3-2. 
The minimum drag coefficient was 0.006 at a of 1°. The drag increased to a maximum value at 
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about 0.28 at a of 20°. The maximum lift to drag ratio for the swept wing was 27.5 at a of 3.2°. 
Pitching moment about the 25% MAC location was nearly flat, as shown in figure 3-34(e), until 
about stall where a considerable increase in negative pitching moment (leading edge down) was 
observed as the load center moved downstream of the quarter-chord point due to flow separation. 
Pitching-moment coefficient ranged from -0.187 to 0.15. The pitching moment is slightly stable 
and there is a tendency for pitch up at a of 10° to 12°. The pitching moment then breaks with 
positive stability, airplane nose down, following stall and ct in the range 15° to 17°. 

Hinge-moment coefficients are presented in figure 3-34(f) and in table 3-3. The maximum slope 
of the hinge-moment coefficient occurred in region B and was -0.0155 per degree. This 
corresponds to a change in hinge moment of 0.52 in-lb per degree for the wing model tested. For 
a control deflection of 0°, the hinge-moment coefficient was negative (i.e., leading edge down) 
for positive C as expected. The graph of hinge moment coefficient versus a is linear in the 
region corresponding to the lin ar part of the lift curv e and breaks as maximum lift is approached 
and the hinge moment slightly reverses at a larger than 15°. 

3.2.3.2 Icinc Research Tunnel Ice Shane Castings. 

3 .2 .3 .2.1 Lift Coefficient . 

Ice shapes caused significant changes in the wing lift characteristics, as shown in figures 3-34(a) 
and 3-34{b). With the exception of the IRT-SC5 ice shape, which increased G h.i, and clh.ii with 
respect to the clean wing, all ice shapes reduced lift throughout the a -range. In addition, the iced 
wing ctaaii and the linear lift s'ope were reduced with respect to the clean configuration. 
Table 3-1 provides lift performance for the clean and iced wing cases for angles of attack 
corresponding to the linear and nonlinear portions of the lift curve. Tabic 3-2 compares 
and clh.ii for all six ice shapes tested The Ci^oi reduction was in the range of 1 1.5% to 93.6% 
and the corresponding reduction in a*ua ranged from 8% to 56%. The largest degradation in lift 
performance was obtained with the 1RT-CS22 ice shape. For this ice shape, the large upper and 
lower horns near the wing leading edge caused extensive flow separation, even at low a. The 
improved lift performance observed with the IRT-SC5 ice shape was mainly due to the small 
leading-edge droop (LE flap effect) caused by the ice shape. The pressure distributions at the 
50% S'-mispan station presented in figure 3-36 demonstrate that for a of 8 , the suction peak for 
the IRT-SC5 case was lower than the clean wing, indicating a lower effective a. In addition, 
figure 3-36(f) shows that at a of 16°, the baseline experienced extensive flow separation, while 
flow separation for the wing with the IRT-SC5 ice shape was considerably reduced. 

The aerodynamic performance data for the IRT-CS2, IRT-CS10, and IRT-CS22 complete 
scallop glaze ice shapes indicate that the ice s&apec with the larger horn size resulted in greater 
penalties in lift characteristics, as shown n Agues 3-34(a) and 3-34(b). The iced wing 
performance in terms of ou.n and Cu.uii was progressively reduced in the following ice shape 
sequence: IRT-CS2, IRT-CS1Q, and IRT-CS2? 
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3.2.3.2.2 Drag Coefficient . 


The increase in the drag coefficient due to the ice shapes tested is demonstrated in 
figures 3-34(c) and 3 -34(d) and in tables 3-1 and 3-2. In general, the increase in minimum drag 
due to the ice accretions ranged from 13 to 35 times that of the clean wing. The smallest drag 
rise was caused by the IRT-SC5 ice shape, while the largest drag increment was due to the 
IRT-CS22 wing. Near the clh.m of the clean wing, the drag due to the ice shapes was 17% to 
104% greater than the clean wing. Furthermore, as the ice shape horn size increased (from 
1RT-CS2 to 1RT-CS22 ice shapes), the drag increment of the iced wing increased, dire to 
increased flow separation downstream of the ice shape. 

3.2.3.23 Pitching-Moment Coefficient . 

All pitching moment data presented in figure 3-34(e) are about the 25% MAC point The 
addition of the ice shapes caused considerable changes in Chmaca due to the shift in the load 
distribution caused by the separated flow downstream of the ice shapes. In general, for positive 
a, the ice shapes caused more positive pitching moment prior to stall compared to the clean wing 
case. Thus, die lift vector for the iced wing was upstream of the 25% MAC location. The ice 
shap .auses unstable behavior, with a break at a of 6°. Stability increases for a greater than 6°. 
Note that a 0.01 change in the value of the pitching-moment coefficient corresponds to a change 
of 5.7 ft-lb in pitching moment about the MAC/4 point 

3.2 1 .2.4 Hinge-Moment Coefficient . 

The differences observed between the clean and iced wing cases in figure 3-34{f) were mainly 
due to the increased separation over the aileron upper surface caused by the ice shapes. In 
general, the ice shapes moved the start of region B to the left (lower a) and increased the hinge- 
moment coefficients over region B (Ch became more nega,.vs). In all cases, the maximum hinge 
moment for the iced wing in region C was bounded by the maximum Ch of the clean wing, as 
shown in figure l*-34(f). Note a shift in the entire Ch versus a plot A control force reversal 
(i.e., change in Ch from positive to negative) was observed between a of 0° and 3° for the 10- 
and 22.5-min ice shapes. This was caused by increased flow separation over the lower surface of 
the wing that resulted in greater suction over the control lower surface. Thus, the aileron had the 
tendency to move trailing edge down for a between 0° and 3°. 

The effect of hom height on aileron hinge moments can be assessed by reviewing the results 
obtained with the 2 -min (IRT-CS2), 10-min (IRT-CS10), and 22.5-min (IRT-CS22) complete 
scallop glaze ice shapes. As demonstrated in table 3-3, the slope of the Ch curve in region B 
decreased as the ice hom height was decreased. However, in regions A and C, the Ch slope 
increased as the ice shape hom height was increased. 

? y jJU Pressure Distributions . 

Wf *► -ibutions corresponding to wing sections at the 15%, 50%, and 85% semispan 
'v .■ presented in figures 3-35, 3-36, and 3-37, respectively. The pressures arc for a of 
.«'■ '2°. and 16°, which cover the linear and nonlinear lift range. The results indicate 
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that, in most cases, the addition of ice shapes resulted in a dramatic change in the clean wing 
pressure distribution. Near the wing LE, surface pressures for the iced wing featured a region of 
separated flow, which was characterized by a flat pressure distribution followed by pressure 
recovery. The region under the flat curve in the pressure data indicates the presence of a 
separated flow bubble. The extent of bubble was a function of ice shape and a. In many cases, 
particularly near the 50% and 85% span wise stations, massive flow separation was observed as a 
was increased. Massive flow' separation was associated with flat pressure distributions that 
extended to the wing TE. At high Positive a, separated flow was observed over the wing upper 
surface. For ice shape with the large horns, extensive flow separatum occurred over both wing 
surfaces at low a. Flow separation caused considerable changes in the load distribution over the 
wing and the control surface, which was the main reason few the observed changes in 
aerodynamic performance and aileron hinge moment with respect to the clean wing. Pressure 
trends for the iced configurations and for each a presented in figures 3-35, 3-36, and 3-37 are 
summarized below: 

• a = -4°: In general, suction was observed over the lower surface of the wing and was 
associated with a leading-edge bubble. Few all ico shapes, the bubble extent increased 
towards the outboard sections of the wing. In most cases, the highest suction occurred 
near the inboard sections of the wing (i.e., 0% to 30% of semispan). The IRT-CS22 ice 
shape caused massive flow separation at practically all spanwise stations. 

• a = 0°: Suction remained higher over the lower surface for practically all ice shapes. 
Once again, LE bubbles were evident at all three spanwise locations for most of the ice 
shapes tested. The bubble extent was a function of hom size. Large horns resulted in 
longer bubbles. For the IRT-CS10, IRT-IS10, 1RT-CS22, and IRT-IPSF22 ice shapes, 
the maximum suction occurred at spanwise stations loetted between 15% and 50% 
semispan. For the small ice accretions, namely 1RT-SC5 and IRT-CS2, the maximum 
suction took place over the outboard part of the wing. Note that at a of 0°, the wing tip 
was at a geometric a of -4° due to the wing twist 

• a = 4°: For this angle of attack, suction for most of the ice shape cases was increased 
over the wing upper surface. The Cp on the lower surface varied from negative near the 
LE, indicating die presence of separation bubbles, to low negative and, in some cases, 
positive near the TE. The extent of the upper surface bubbles varied with spanwise 
distance from root to t'~ H e to the negative wing twist which resulted in a lower auxat 
over the outboard sect ns. For the 10- and 22.5-min ice shapes, the extent of the 
separation bubbles was, in most cases, greater over the middle and outboard sections of 
the wing. For these ice shapes, maximum suction did not vary significantly with 
spanwise location. 

• a = 8°: Upper surface suction and bubble extent was increased at this a at all spanwise 
stations for ice shapes 1RT-CS10, IRT-IS10, IRT-CS2, IRT-CS22, and IRT-IPSF22. For 
the 10- and 22.5-min ice shapes, extensive flow separation was observed at the 85% 
semispan station. The 5-min rime ice shape, IRT-SC5, did not exhibit bubble formation 
over the two inboard stations. However, at the outboard station (85% semisp-n), a small 
bubble extending to about 20% chord was observed. Maximum upper surface suction 
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took place over the wing near 15% semispan. The only exception was the 1RT-CS22 
case tor which the maximum suction occurred at the 85% semispan station. With the 
exception of the 1RT-CS22 ice shape, all other ice shapes had positive or low negative 
pressure coefficients over the lower surface of the wing. 

• a = 12°: Large bubbles were observed at the 15% semispan station and complete flow 
separation was evident over the wing upper surface at the 50% and 85% semispan 
locations with the 10- and 22.5-min ice shapes tested. The flow over the 1RT-SC5 and 
IRT-CS2 ice shapes exhibited smaller bubbles compared to the larger ice shapes. The 
bubbles for these two ice accretions extended from about 20% to about 75% chord 
length, depending on ice shape and span wise station. For all ice shapes, maximum 
suction was observed at the inboard station located at 15% semispan. 

• a = 16°: At this angle of attack, massive flow separation was evident over the wing 
upper surface at all span wise stations for all 10- and 22.5-min ice shapes. In all cases, 
maximum suction took place at the 15% semispan station. The 2- and 5 -min ice 
accretions exhibited long bubbles over the inboard stations. However, at the outboard 
stations, the flow for these two ice shapes was completely separated over the wing upper 
surface. 

The impact of the horn size on aerodynamic performance can be explained by reviewing the 
pressure distributions. Specifically for ice shapes with large horn heights, LE separation bubbles 
were observed at low a. These bubbles became progressively longer in the chordwise direction 
ana eventually burst to form a region of massive flow separation. Ice shapes of smaller horn size 
formed bubbles with smaller chordwise extent. These bubbles did not appear until a was higher 
with respect to where the bubbles were observed with the larger horn ice shapes. With the IRT- 
CS22 configuration, extensive flow separation occurred over the wing suction surface even at 
low a. Although the flow over the upper surface of wing with the IRT-CS22 ice shape remained 
separated throughout the positive a -range, the lift increased as a was increased beyond about 6°. 
This was due to the increase in the pressure on the lower surface of the wing where the flow 
remained mostly attached. Pressure on the upper surface did not vary significantly with a once 
the wing had stalled. 

3.2.3.3 Smooth LEWICE Ice Shanes . 

The effects of the smooth LEWICE ice shapes on wing performance are demonstrated in 
figures 3-38 to 3-41 and in tables 3-4 to 3-6. 

3.2.3.3.1 Lift Coefficient . 

With the exception of the LS-SC5 ice shape, which increased C L ,*»u and ewi, all ice shapes 
reduced lift performance with respect to the clean wing, as demonstrated in figures 3-38(a) and 
3-38(b) and in tables 3-4 and 3-5. In assessing these reductions note that they could be larger if 
clean wing data were available for higher Reynolds numbers. The iced wing lift cu’ ves exhibited 
reduced lift coefficient and lower lift slope with respect to that of the clean wing. The reduction 
of C i rf .n ranged from 18.4% to 43.7%. The largest loss in lift was obtained with the LS-CS22N 
ice shape. As shown in table 3-5, smooth LEWICE ice shapes LS-CS10, LS-SC5, LS-CS2, and 
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LS-IPSF22 increased by 5.8% to 7.2%, while LS-IS10, LS-CS22N, and LS-CS22S ice 
shapes reduced 0^.11 by 1 .4% to 2.9%. The 5 -min rime ice shape LS-SC5 increased Cl^uu by 7% 
and increased a«uii by 1°. The observed performance increment in the iced wing lift was due to 
the leading-edge droop (i.e., leading-edge flap effect) caused by the 5-min rime ice shape. 

The effect of horn height on the lift performance can be observed by comparing the iced wing lift 
performance with the LS-CS2, LS-CS10, LS-CS22N, and LS-CS22S ice shapes. The Cl^uu 
reduction with respect to the clean wing was -13.8%, -18.4%, -43.7%, and -38% for the LS-CS2, 
LS-CS10, LS-CS22N, and LS-CS22S ice shapes respectively. The corresponding change in ou.n 
was 6.5%, 5.8%, - 22 %, and -2.9% for the same sequence of ice shapes. 

3-2.3.3.2 Drag Coefficient . 

In general, all smooth LEWICE ice shapes increased wing drag considerably, as demonstrated in 
figures 3-38(c) and 3 -3 8(d) and in tables 3-4 and 3-5. The results presented show that the 
increase in C n m . due to the ice shapes ranged from 66.7% for ice shape LS-SC5 to 2366.7% for 
ice shape LS-CS22S with respect to the clean wing. For a near the clean wing a ^.n. drag rise for 
the iced wing was in the range of -15.6% for the LS-SC5 to 100% for LS-CS22S with respect to 
the clean case. The effects of the LEWICE ice shape horn height on Co was similar to that 
obtained with the IRT ice shapes, i.e., ice shapes with larger horn heignts produced greater drag 
increments. Moreover, the larger horns of the LS-CS22S ice shape produced a larger drag 
increment with respect to the clean wing at lower a in comparison to the LS-CS22N ice shape. 
As a was increased, the difference in drag performance between these two 22.5-min ice shapes 
diminished. 


3.2.3.3.3 Pitching-Moment Coefficient . 

Fran figure 3-38(e), it is observed that the pitching moments about the 25% MAC point of the 
iced wing were considerably more positive or more negative than the clean wing, depending cm 
a. This was the result of a shift in the center of the aerodynamic load due to flow separation. In 
general, flow separation downstream of the iced shape resulted in a front-loaded pressure 
distribution caused by the leading-edge separation bubbles. Another contributing factor to the 
observed increase in iced wing pitching moment was the aerodynamic load on the ice shape, 
which increased as a was increased. The largest increment in Cm was obtained with the LS- 
CS22N ice shape. 

Pitching moment for the LS-SC5 configuration was comparable to the clean wing, particularly 
for a in the range of -4° to approximately 12°. As a was increased above 12°, however, the LS- 
SC5 ice shape produced more positive or more negative Cm than the clean case. The similarities 
and differences in the Cm curves of the clean and LS-SC5 cases can be explained by reviewing 
the Cp distributions shown in figures 3-39 to 3-4 1 . For a of 4° and 8°, LS-SC5 produced a lower 
suction peak than the baseline at all three spanwise stations. However, pressure recovery 
downstream of the suction peak and the Cp distributions over the lower surface were similar to 
the clean case. As a result, iced wing Cm at 4° and 8° a was comparable to the chan wing. On 
the other hand, as a was increased to 12°. the iced configuration (LS-SC5) ban a greate. suction 
peak than the clean case *4 50% semispan and over the outboard wing stations. As a was further 
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increased to 16°, while both configurations (clean and LS-SC5) experienced extensive flow 
separation near the wing tip, the LS-SC5 ice shape resulted in a separation bubble on the wing 
suction surface and produced higher leading-edge suction than the clean case, as shown in 
figures 3-39(f) and 3 -40(f). Consequently, the LS-SC5 case resulted in higher Cm for a greater 
than 12°. 

3.2.3.3.4 Hinge-Moment Coefficient . 

Hinge-moment coefficients for all smooth LEW1CE ice shapes tested are presented in 
figure 3-38(0- Hinge-moment coefficient slope (Chq) for the clean and iced wing are compared 
in table 3-6. In general, hinge moments for the iced wing were more positive or more negative 
with respect to the clean case, particularly in regions A and B of the hinge-moment curve. The 
observed increase or decrease in hinge moment was the result of flow separation over the aileron 
due to the ice shapes. In most cases, the start of region B in the iced wing C H curve was at a 
lower a, and the slope of Ch was less than the clean wing. Note that the maximum and 
minimum hinge moments for all ice shape cases were bounded by the maximum and minimum 
Ch of the clean wing. 

3.23.3.5 Pressure Distributions. 

Cp distributions are presented in figures 3-39 to 3-41 for the 15%, 50%, and 85% semispan 
stations. Near the wing tip, extensive flow separation was observed ~ver the wing suction 
surface for most iced configurations, as shown in figure 3-41. Considering the wing pressure 
distributions with the glaze ice shapes LS-CS2, LS-CS10, LS-CS22N, and LS-CS22S, ice shapes 
with larger horns produced longer bubbles, as was the case with thei r IRT ice shape counterparts. 
In general, the LS-CS22N and LS-CS22S cases had comparable Cp distributions at all stations. 
However, for 12° and 16° a, the LS-CS22S had higher pressure over the wing lower surface and 
hence generated higher lift 

3.23.4 Rough LEWICE Ice Shanes . 

The effects of the rough LEWICE ice shapes on wing performance are demonstrated in 
figures 3-42 to 3-45 and in tables 3-7 to 3-9. The results showed that in most cases, die rough 
LEWICE ice shapes resulted in greater reductions m Cmuu and a— » and larger drag increments 
compared to the smooth LEWICE ice shapes. Roughness has a significant inpact on the 
location of flow separation, particularly for large glaze ice chapes where the flow at high 
Reynolds numbers usually separates near the horn tip. The size of the separated region 
downstream of the ice shape is a function of the locution of the separation point. Rough and 
smooth ice shapes can result in considerably different aerodynamic effects, particularly when the 
radius of the hom tips is large. For smali ice shapes that do not exhibit horn features such as 
rime ice shapes, roughness may actually increase aerodynamic performance, depending on the 
Reynolds number. In this case, the location of flow separation is affected by the state of the 
boundary layer. At low Reynolds numbers, roughness can increase the energy levels in 
boundary layer and delay flow separation, resulting in improved aerodynamic performance such 
s higher Cl^uu and cu.ii. 
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3.2.3.4.1 Lift Coefficient . 

Note that effects of higher Reynolds numbers representative of flight on the dean wing are not 
known. In most cases, the rough LEWICE ice shapes decreased Cl %• *ughr-ut the a-range, as 
demonstrated in figures 3-42(a) and 3 -42(b). One exception to this observation was the 5-min 
rime ice shape LR-SC5 that resulted in improved lift performance with respect to the clean wing. 
Tables 3-7 and 3-8 provide lift coefficients for selected angles of attack for the clean and iced 
wing cases. The reduction in Cl^uu was 30%, 25%, -12.6% (increase), 17%, 53%, 61%, and 
30% for the LR-CS10, LR-IS10, LR-SC5, LR-CS2, LR-CS22N, LR-CS22S, and LR-IPSF22, 
respectively. The reduction in a«aii for the glaze-iced configurations ranged from 0.1° to 3.5°, 
depending mi ice shape. The 5-min rime ice shape increased the clean wing Oaui by 2°. With 
the exception of the 5-min rime and 2-min glaze ice shapes, all other rough LEWICE ice shapes 
reduced the linear lift slope of the clean wing. 

Observe that ice shapes with the larger horn size resulted in greater penalties m lift performance, 
Le., LR-CS2 ice shape had higher Ostui and Cl^uii than LR-CS10, LR-CS22N, and LR-CS22S 
cases. Effects of horn height on Cl were similar to the trend of IRT ice shapes. 

3.2.3.4.2 Drag Coefficient . 

Table 3-7 lists Cd of clean and iced configurations at angles of attack of 3°, 5°, 10°, and 15°, 
while table 3-8 shows the effect of rough LEWICE ice shapes on Cn m ,n. Graphically, the Cd and 
ACd curves are illustrated in figures 3 -42(c) and 3 -42(d), respectively. All roughened LEWICE 
ice shapes increased Cd throughout the a-range. The increase in C n mm with respect to the clean 
case ranged from 1 16.7% for LR-SC5 to 3316.7% for LR-CS22S. In terms of the horn height 
effects cm Cd, the trends were similar to that obtained with the IRT ice shapes, i.e., ice shapes 
with larger hom heights produced greater drag increments. Moreover, the larger horns of LS- 
CS22S ice shape produced greater drag at lower a in comparison to the LS-CS22N ice shape. 
But, as a increased, the difference in drag performance diminished. 

3.2.3.4.3 Ph ^ mg-Moment Coefficient . 

The *rends 1 m for the rough LEWICE ice shapes were similar to that obtained with the smooth 
LEWICE iu. hapes, as demonstrated in figure 3 -42(e). For positive a less tl an about 15°, flow 
separation downstream of the ice shapes caused the lift center to move upstream c r *he 25% 
MAC location, which led to more positive (leading edge up) pitching moment than ne clean 
wing. However, for greater a, most of the iced wing cases experienced a large change, from 
positive to negative, in pitching moment. This was attributed to the flat pressure distribution 
associated with massive flow separation over the suction surface ot the wing. For the glaze ice 
shapes LR-CS22N and LR-CS22S, the pitching moment remained positive throughout the 
positive a-range. For these two large ice shapes, extensi\e flow separation was present over 
both upper and lower wing surfaces for practically all positive a and the lift center did not move 
upstream of the 25% MAC location. 
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3.2.3.4.4 Hinge-Moment Coefficient . 


AiLron hinge moments for the iced wing configurations are depicted in figure 3 -42(f). Table 3-9 
provides the slopes and a-range for regions A, B, and C of the Ch curves. The results show that 
the LR-SC5 ice shape had a Ch curve similar to the clean case for practically the complete range 
of angles of attack. The ice shapes LR-CS22N and LR-CS22S resulted in large changes from the 
clean wing aileron hin g e moments. In feet, for a in the rang^ of 6° to 8° for the LR-CS22N case 
and 6° to 1 1° for the LR-CS22S ice shape, the aileron hinge moment was nearly zero, indicating 
zero stick force and possibly aileron float. Force reversal was the result of nearly equal pressure 
distributions over the suction and pressure sides of the aileron surfaces. The other four rough 
LEWICE ice shapes (LR-CS10, LR-IS10, LR-CS2, and LR-IPSF22) cause mc.^ positive or 
more negative hinge moments with respect to the clean wing over regions A and B of the Ch 
curve. Once again, the iced wing hinge moments were bounded by the clean wing Eh limits. 
One exception was the 5 -min rime ice shape that resulted in slightly higher hinge moments with 
respect to the clean wing in region C of the Ch curve for positive a. 

3.2.3A5 Pressure Distributions. 

Effects of rough LEWICE ice shapes on pressure distributions are presented in figures 3-43 to 
3-45 for the 15': \ 50%, and 85% semispan stations. All glaze ice shapes produced separation 
bubbles downstream of the ice horns, which reduced the suction peak near the leading edge of 
the clean wing. Among all glaze ice cases, LR-CS2 caused the least penalties on aerodynamic 
performance due to the feet that its ice horns were smaller. Consequently, for positive a, the 
separation bubble on the suction wing surface had shorter extern. Pressure distributions aft of the 
bubble and on the pressure surface were comparable to the e’eae wing. As the ice shape horn 
height increased (such as LR-CS10, LR-CS22N, and LR-Ch?2S configurations), the size of the 
separation bubble was increased. 

3.2/ Compan/on of IRT and LEWICE Ice Shapes 

Figures 2-33 to 2-38 demonstrate that the differences in the geometric features of the LEWICE 
and IRT ice shape sections were considerable for the CS10, IS 10, CS22, and IPSF22 cases. 
These differences inchided overall shape of ice section, height of the upper and lower ice horns, 
hom angle with respect to the ho-izontrt and location of hom tip, and root with respect to the 
leading edge of *he wing. For example, at stations A and B, the upper and lower hom tips of 
IRT-CS10 were turther upstream of the LEWICE ice saape (LS-CS10) hom tips, as 
demonstrated in figure 2-33. In addition, it these two stations, the hom angles with icspect to 
the horizontal were, in most cases, greater for the LEWICE ice shape. Note, however, that the 
upper and lower hom heights of the LEWICE and the IRT ir* shape sections were similar for 
both span wise stations. At station C, the LEWICE ice shape upper and lower horns extended 
further upstream and downstream with respect to the horns of the IRT casting. In addition, the 
height (distance of hom tip from airfoil LE) of the LEWICE horns was greater. Other significant 
differences between the simulated and the IRT ice shapes included the intricate roughness and 
feather-like features seen on the surface of the IRT ice shape castings present? ! in figure 2-33. 
To assess the effect of these geometric differences in the actual and simulated ice shapes on 
aerodynamic performance, experimental data obtained with the castings and the cor expending 
smooth and rough LEWICE ice shapes are compared in figures 3-46 to 3-69. 
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FIGURE 3-58. EFFECT OF IRT-SC5, LS-SC5, AND LR-SC5 ICE SHAPES (ICING 
CONDITION 3) ON PRESSURE DISTRIBUTIONS AT 15% SEMISPAN; 

Re= 1.8x10 s ; 6 A = 0° 














FIGURE 3-59. EFFECT OF IRT-SC5, LS-SC5, AND LR-SC5 ICE SHAPES (ICING 
CONDITION 3) ON PRESSURE DISTRIBUTES AT 50% SEMISPAN; 

Re= 1.8x10 6 ;5 a = 0° 
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FIGURE 3-60. EFFECT OF IRT-SC5, LS-SC5, AND LR-SC5 ICE SFIAPES (ICING 
CONDITION 3) ON PRESSURE DIST 7 WTIONS AT 85% SEMISPAN; 

Re= 1.8xlC - 0° 
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FIGURE 3-63. EFFECT OF IRT-CS2, LS-CS2, AND LR-CS2 ICE SHAPES (ICING 
CONDITION 4) ON PRESSURE DISTRIBUTIONS AT 85% SEMiSPAN; 

Re= 1.8x10 6 ;5 a = 0° 
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FIGURE 3-66. EFFECT OF IRT-CS22, LS-CS22N, LS-CS22S, LR-CS22N, AND LR-CS22S 
ICE SHAPES (ICING CONDITION 5) ON PRESSURE DISTRIBUTIONS AT 85% 

SEMISPAN: Re= 1.8x10°; S A = 0° 






































3.2.4.1 Icing Condition 1 . 


Lift, drag, and hinge-moment < ~ efficients for this case are presented in figure 3-46. The results 

indicate the following: 

• Lift coefficients for the 1RT, smooth, and rough LEWICE ice shapes were in good overall 
agreement for a in the range of -7° to 7°, as shown in figure 3-46(a). All ice shapes 
resulted in a lower lift slope with respect to the clean wing. Over the stall region, the lift 
curves for the rough LEWJCE ice shape (LR-CS10) and that of the IRT casting (IRT- 
CS10) were in good correlation. The smooth LEWICE ice shape (LS-CS10), however, 
resulted in higher C t.^ i and cwi. The reduction in G^uii with respect to the clean wing 
for the IRT-CS10, LS-CS10, and LR-CS10 ice shapes was 38%, 18%, and 30% 
respectively. The change in a*uu with respect to the clean wing case was -24%, 5.8%, 
and -1.4% for ice shapes IRT-CS10, LS-CS10, and LR-CS10 respectively. 

• All ice shapes caused a significant increase in Cd (e.g., 1 100% to 1200% in Cr> m in) with 
respect to the clean wing, as demonstrated in tables 3-1, 3-2, 3-4, 3-5, 3-7, and 3-8. Drag 
coefficients obtained with all three ice shapes were in good overall agreement. However, 
the drag due to LEWICE ice shapes was, in general, higher than the IRT ice shape 
casting, as shown in figure 3-46(c). For a greater than 10°, the smooth LEWICE ice 
shape produced higher drag compared to the rough LEWICE shape. This was the result 
of more extensive flow separation downstream of the smooth LEWICE ice horns. In 
general, the size and extent of flow separation downstream of large glaze ice shapes are 
affected by the roughness of the hom tip. 

• Pitching-moment characteristics of clean and iced configurations are presented in 
figure 3 -46(e). The results show that all three ice shapes produced pitching moments that 
were considerably higher or lower than the clean wing, depending on a. The trends in 
pitching moment with a for the ice shapes IRT-CS10, LS-CS10, and LR-CS10 were 
similar. As a was increased from -8° to a*uu, pitching moment was increased from 
negative (nose down) to positive (nose up) for all iced wing cases. This was mainly due 
to the leading-edge bubbles that resulted in a forward shift (ahead of the 25% MAC 
point) of the lift vector. When the bubble was on the lower surface (negative a), lift was 
negative, and the lift vector was ahead of the 25% MAC, resulting in negative Cm. 
Positive a had the opposite effect on lift and Cm- Past stall, Cm exhibited a downward 
trend as the load center moved aft of the 25% MAC point. This was due to massive flow 
separation over the wing, which resulted in a flat pressure distribution over the upper 
surface of the wing. Notable differences in the magnitude of the pitching moment of the 
IRT-CS10 ice shape compared to the LEWICE ice shapes were observed for a in the 
range of 6 to about 15°. These differences were, again, a function of the extent of the 
flow separation over each ice shape. 

• Aileron hinge-moment coefficients for the IRT and LEWICE ice shapes are presented in 
figure 3-46(f). In all cases, the hinge moments due to the ice shapes were within the 
range defined by the minimum and maximum hinge moments of the clean wing. In 
general, the increase in hinge moments cursed by the LEWICE ice shapes was greater 
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than the 1RT ice shape for a in the range of -3° to 20°. Note that the Cm* in region A of 
the LEW1CE ice shape curves was considerably higher than the 1RT ice shape. 

3. 2.4.2 Icing Condition 2 . 

Lift, drag, and hinge-moment coefficients for xhis ice shape configuration are presented in 

figure 3-47. The results indicate the following: 

• Lift coefficients for the 1RT, smooth, and rough LEWICE ice shapes were in good overall 
agreement for a in the range of -7° to 12°, as shown in figure 3-47(a). The lift slope for 
all three ice shapes was the same as the clean wing. The reduction in Cl.suii with respect 
to the clean wing for the IRT-IS10, LS-IS10, and LR-IS10 cases was 26.4%, 21.8% and 
25.3% respectively. The reduction in a^nii with respect to the clean wing was 23%, 1.4%, 
and 1.4% for ice shapes IRT-IS10, LS-IS10, and LR-IS10 respectively. 

• Drag coefficients for all three ice shapes were in good overall agreement for a in the 
range of 6° to 20°, as shown in figure 3-47(c). For a less than 6°, the IRT-IS10 ice shape 
resulted in higher drag loan that obtained with the smooth and rough LEWICE ice 
shapes. The increase in Cn,mm for the iced wing with respect to the clean wing case was 
683%, 500%, and 550% for the IRT-IS10, LS-IS10, and LR-1S10 ice shapes respectively. 

• Pitching-moment coefficients for the IRT-IS10, LS-IS10, and LR-1S10 ice shapes are 
shown in figure 3-47(e). For a in the linear lift range, the results show that all ice shapes 
exhibited similar trends in Cm behavior. Specifically, a maximum Cm, of approximately 
0.03 was attained at a of about 7°. However, beyond 7° AOA, the wing with the IRT- 
IS10 ice shape exhibited a drop followed by a small increase in Cm, indicating that the lift 
vector was moving back and forth. The LEWICE ice shapes resulted in a flat pitching- 
moment curve with positive Cm values for a in the range 7° to about 14°. After wing 
stall, all ice shapes manifested decreasing Cm (more negative) with a. 

• Aileron hinge-moment coefficients for the 1RT and LEWICE ice shapes are presented in 
figure 3-47(f). The ice shapes increased aileron hinge moments for a in the range of -6° 
to 13°. In all cases, however, the hinge moments due to the ice shapes were within the 
range defined by the minimum and maximum hinge moments of the clean wing. The 
increase in the iced wing hinge moment with respect to the clean wing for the LEWICE 
cases was greater than that obtained with the IRT ice shape for a in the range of 5° to 
10 °. 

3.2.4.3 Icinu Condition 3 . 

The IRT and LEWICE ice shapes were very similar in shape and size, as demonstrated in 

figure 2-35. Aerodynamic performance of the GLC-305 wing with the IRT and LEWICE rime 

ice shapes is presented in figure 3-48. The data presented indicate the following: 

• Note that effects of higher Reynolds numbers representative of flight on the clean wing 
are not known. The lift coefficients for the clean and all iced wing cases were in good 
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overall agreement for a in the range of -8° to 12°, as shown in figure 3-48(a). The ice 
shapes tested resulted in 3.4% to 12.6% higher Cl,^ and 7.2% to 14.5% higher a«uu 
with respect to the clean wing. The main reason for the observed increase in Cl,««u and 
cl*. ii for the iced wing cases was the wing leading edge droop (see figure 2-35) feature of 
these three ice shapes. Note that the rough LEWICE shape resulted in the largest lift 
performance increment. This was probably due to delayed boundary layer separation 
caused by the grit roughness. 

• The 1RT and rough LEWICE ice shapes resulted in higher drag with respect to the clean 
wing for most of the a-range tested. The smootn LEWICE ice shape (LS-SC5) increased 
drag for a up to 10°, as shown in figure 3-48(c). At higher angles of attack, however, the 
drag of LS-SC5 was in most cases lower than that of the clean wing. Th i increase in the 
iced wing Cry™- with respect to the clean wing case was 133%, 67%, and 117% for the 
IRT-SC5, LS-SC5, and LR-SC5 ice shapes respectively. 

• From figure 3-48(e), it is observed that Cm was gradually increased with AOA prior to 
stall. But, as the IRT-SC5 and LR-SC5 configurations stalled. Cm decreased 
considerably and eventually became considerably more negative compared to the clean 
wing. As for LS-SC5 ice shape, its Cm distribution maintained more positive Cm than 
that of the clean case for angles of attack greater than 7°. 

• Aileron hinge-moment coefficients for the 1RT and LEWICE ice shapes were very 
simil ar to that of the clean wing for angles of attack greater than -4°, as shown in 
figure 3 -48(f). Between a of -8° and -4°, the ice shapes resulted in lower aileron hinge- 
moment coefficient. 

3.2.4.4 Icing Condition 4 . 

The IRT and LEWICE ice shape sections presented in figure 2-36 were in good overall 
agreement. Note that the cross section area of LEWICE ice shape at station A was larger than 
that of the IRT ice shape casting. Aerodynamic performance degradation due to the IRT and the 
LEWICE ice shapes is presented in figure 3-49. The experimental lift, drag and, hinge-momem 
coefficients presented indicate the following: 

• The lift coefficients of the clean and iced wing cases were in good agreement over the 
linear lift range corresponding to angles of attack between -6° and 9°, as shown in 
figure 3-49(a). The reduction in the Cl,suk with respect to the clean wing caused by the 
IRT, smooth LEWTCE, and rough LEWICE ice shapes were 1 1.5%, 13.8%, and 17.2% 
respectively. The corresponding change in a,t«u with respect to the clean wing was -8%, 
6.5%, and -1%. Note that the rough LEWICE shape resulted in the largest lift reduction 
near stall, while the IRT casting caused the largest reduction in ou.n. 

• The IRT and rough LEWICE ice shapes resulted in higher drag with respect to the clean 
wing throughout the a-rangc, as shown in figure 3-49(c). The increase in the iced wing 
C >,min with respect to that of the clean wing was 200%, 183%, and 233% for the IRT- 
CS2, LS-CS2, and LR-CS2 ice shapes respectively. 
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• Pitching moment trends for all three ice shapes were very similar. In general, the iced 
wing Cm was more negative or more positive than that of the clean wing throughout the 
a -range. The difference between he clean and the iced wing Cm ranged from -0.04 to 
0 . 02 . 

• The aileron hinge moment behavior for the clean and IRT-CS2 cases were very similar 
for angles of attack in the range of -7° to 14°. The LEW1CE ice shapes increased aileron 
hinge moments with respect to the clean case over region B of the Ch curve. For all ice 
shapes, the hinge-moment coefficients were within the range defined by the hinge 
moments of the clean wing. 

3.2.4.S Icing Condition 5 . 

The IRT and LEWICE ice shape sections presented in figure 2-37 exhibited large differences in 
bom shape, size, and overall section area. In addition, the IRT-CS22 ice shape had complex 
roughness and scallop features not present in the LEWICE ice shapes. Two sets of smooth and 
rough LEWICE ice shapes were tested for this icing condition as discussed in section 2.2.3.2. 
Cne was based on computations using sections normal (N) to the wing leading edge while the 
other was computed using stream wise (S) sections. The wing aerodynamic performance data for 
the IRT, the smooth, and rough LEWICE ice shapes presented in figure 3-50 indicate the 
following: 

• The iced wing performance exhibited large reductions in lift, and in lift slope with 
respect to that of the clean wing. The largest degradation in lift was caused by the IRT- 
CS22 ice shape. The rough LEWICE ice shapes caused larger reductions in lift than their 
smooth counterparts. The rough LEWICE ice shape LR-CS22S obtained using 
streamwise wing sections and free-stream velocity resulted in larger lift degradation than 
that obtained with the rough LEWICE 'ce shape LR-CS22N, which was computed with 
the normal component of the free-stream velocity and airfoil sections normal to the wing 
leading edge. The difference in lift reduction due to the IRT-CS22 ice shape and that 
caused by the rough LEWICE ice shapes was considerable, particularly for angles of 
attack in the range of 0° to approximately 13°. The reduction in the Cl,suu with respect to 
the clean wing caused by the IRT-CS22, LS-CS22N, LS-CS22S, LR-CS22N, and LR- 
CS22S ice shapes was 93.6%, 44%, 38%, 53%, and 61% respectively. The 
corresponding reduction in a...n with respect to the clean wing was 56%, 2%, 3%, 17%, 
and 25%. 


• All ice shapes tested caused very large drag increments, as shown in figure 3-50(c). The 
increase in the CD,mm for the iced wing was 3533%, 2033%, 2367%, 2933%, and 3317% 
for the ice shapes 1RT-CS22, LS-CS22N, LS-CS22S, LR-CS22N, and LR-CS22S 
respectively. The LEWICE ice shape that provided the best drag correlation with the 
1RT-CS22 ice shape was the LR-CS22S. 

• Figure 3-50(e) shows pitching moment characteristics of the wing w ith all five CS22 ice 
shapes. Separated flow downstream of both upper and lower ice horns on the 1RT-CS22 
ice shape led to small changes in Cm compared to the clean case, despite the large 
differences in the wing pressure distributions obtained with the 1RT-CS22 ice shape and 
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the clean wing. The four LEWICE ice shapes produced significantly more negative or 
more positive Cm than that of the clean wing. Also, the difference in CM.ueu.-C’ xwict 
ranged from 0.07 at negative AOA to -0.085 at positive AOA. With increasing a jsitive 
a, the separation bubbles behind the glaze ice horns shifted the load center upstream of 
the quarter-chord location and hence mcreased Cm- Conversely, beyond cuau, separated 
flow over the upper surface and attached flow over the lower surface moved the center of 
pressure downstream and led to decreased Cm- It is worth noting that both LS-CS22N 
and LR-CS22S ice shapes had higher maximum Cm than their CS22S counterparts. Also, 
the break from positive to negative values in the Cm curve occurred at a higher a for the 
CS22N ice shapes. 

• The uleron hin g e moment behavior for the clean and iced wing cases exhibited large 
differences, aj shown in figure 3-50(f). However, for all ice shapes, the hinge-mcmem 
coefficients ’Acre practically within tl . range defined by the minimum and maximum 
hinge moments of the clean v ‘^g. The LEWICE ice shape that provided the best Ch 
correlation with the 1RT casting was the LS-CS22N. 

3.2.4.6 Icing Condition 6 . 

The IRT and LEWICE ice shape sections presented in figure 2-38 demonstrated notable 
differences in size and shape. However, the hom angles and hom heights of the IRT ice shape 
castings and the LEWICE ice shapes were similar. Aerodynamic performance degradation due 
to the IRT and the LEWICE ice shapes is presented in figure 3-51. The experimental lift, drag 
and hinge-moment coefficients presented in these figures indicate the following: 

• The lift coefficient for the iced wing exhibited lower lift and linear lift slope than that of 
the clean wing, as evident in figure i -5 1(a). The reduction in the Ci^tui of the iced wing 
with respect to the clean wing was 39%, 21% and 30% for the IRT, smooth LEWICE, 
t~i rough LEWICE ice shapes respectively. The corresponding change in asuii with 
respect to the clean wing was -24%, 6%, and -9%. Note that the IRT-1PSF22 ice shape 
resulted in the largest reduction of Cl^uii and a*uu with respect to the clean wing. The lift 
curves for the IRT and rough LEWICE ice shapes were in good agreement except near 
stall where the IRT shape caused approximately 1% to 8% more reduction in lift with 
respect to the clean wing. 

• All ice shapes resulted ui significant drag increments with respect to the clean wing drag, 
as shown in figure 3-5 1(c). The increase in Cu.mm of the iced wing with respect to the 
clean configuration was 1200%, 1000%, and 983% for the 1RT-1PSF22, LS-IPSF22 and 
LR-1PSF22 ice shapes respectively. The trends in the IRT and LEWICE ice shapes drag 
curves were similar. However, the IRT ice shape caused higher drag at angles of attack 
in the range of -8° to 1 1°. For a greater than 1 1°, the LEWICE ice shapes resulted in 
higher drag. 

• The trends in Cm for the 1PSF22 ice shapes were quite similar to that obtained with the 
CS10 ice shapes. In general, the LS-IPSF22 and LR.-1PSF22 cases produced similar 
pitching moment characteristics, as demonstrated in figure 3-5 1(c). The observed 
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differences in Cm between 1RT and LEWICE ice shapes i^iuded more positive Cm & 
the LEWICE ice shapes for angles of attack between 6° ana 18° and higher dCM/da afte 
a of approximately 12°. The iced wing cases resulted in considerably more positive Ci 
than the clean v/ing for angles of attack in the range of 1 ° to 11° for the IRT-IPSF22 an 
1° to 15° for the LEWICE ice shapes. 

• The iced wing aiieron hinge moment was within the range defined by the minimum an 
maximum hinge moments of the clean wing for all ice shape cases. For most angles c 
attack, the LEWICE ice shapes caused a larger increase (more negative or more positive 
in hinge moments than that caused by the IRT ice shape. 

2.2. 4.7 Summary . 

With the exception of icing condition 5, the rough LEWICE ice shapes caused similar loss in lit 
compared with the IRT ice shape castings. The difference in iced wing C i «t.n (IRT Cl^uii 
rough LEWTCE C i t.n 'i obtained with the rough LEWICE and IRT ice shapes was -0.07, -0.0) 
-0.08, 0.05, -0.35, -0.28, and -0.08 for CS10, IS10, SC5, CS2, CS22N, CS22S, and 1PSF2 
respectively. In all but one case, the IRT ice shape castings resulted in higher reduction in Cl^h 
compared to that obtained with the rough LEWTCE ice shapes. Again, with the exception c 
icing condition 5, the observed differences in Cl,suii reduction between the IRT and LEWTCE ic 
shapes were in the range of 1% to 9% with respect to the clean wing Q which was equal t 
0.87. In addition, the behavior of the lift curves obtained witn the IRT and LEWICE ice shape 
was very similar. 

In general, the incremeL ' in wing drag caused by the LEWTCE and IRT ice shapes was in goo 
ccrr^lalior. The difference in iced wing Cd,i»j (rough LEWTCE Cn m.„ - IRT for th 

rough LEWTCE and IRT cases w r as 0.004, -0.008, -0.001, 0.002, -0.C36, -0.013, and -0.013 fc 
CS10, IS10, SC5, CS2, CS22N, CS22S, and 1PSF22 respectively. 

With the exception of icing condition 5, the trends in aileron hinge moments obtained with th 
rough LEWTCE and IRT ice shapes were in good agreement. For icing conditions 2 and 3, th 
IRT and LEWTCE Ch magnitudes were in good correlation throughout the a-range. For icin 
conditions 1, 4, and 6, the LEWICE ice shapes resulted in larger hinge moments over regions j 
and B of the Ch curve compared with the IRT ice shapes. For icing condition 5, the behavior c 
Ch for the rough LEWTCE and IRT ice shapes was considerably different. This was mainly du 
to higher suction over the wing upper surface caused by the rough LEWICE ice shapes. 

In view of the differences in the geometric features of the large glaze IRT and LEWTCE ic 
shapes tested, the LEWICE ice shapes provided a good engineering approximation to th 
aerodynamic effects of the IRT ice shapes for five out of the six cases tested. 

5.2.4.8 Pressure Distrib utions— Clean and All Iced Configurations. 

Pressure distributions for the clean wing - .u all 20 ice shape configurations tested are presentc 
in figures 3-52 to 3-69 for a in the range of -4° to 1*° j»- tpai. ise locations correspondin 
to 15% semispan (near wing ir.ot), 50% emis; _„mispan (near wing tip). Th 



pressure data provided insight into the change in wing load distribution with a and the stalling 
behavior of the clean and iced wing. In the discussion below, reference is made to three angles 
of attack, namely the wing geometric angle of attack, the local geometric angle of attack, and the 
effective angle of attack These angles are defined as follows: 

• Geometric angle of attack (a): this is the w ing angle of attack as set on the wind tunne l 
turntable and is the a used in all force and moment plots. 

• Local geometric angle of attack (Otocu): this is the geometric angle of attack at a given 
span wise station, which can be determined from a and the local geometric twist 

• Effective angle of attack (a«): this angle of attack is the sum of 0Uoc*i and the local 
up wash or downwash angle. Upwash or downwash is caused by the wing tip vortex and 
depends on wing lift 

For negative a, the wing upper surface w as at a higher pressure than the lower surface. Due to 
wing twist, which was 0° at the root and -4° at the tip, the tHod varied from root to tip. Thus, at 
negative a, the wing root was at a lower auxai than the wing tip. In other words, the atod 
increased from root to tip. At positive a, the wing tip was at a lower ajod than the wing root. 
However, for positive a, the wing experienced upwash, which increased as the a was increased. 
Thus, for positive lift, the wing upwash was higher near the wing tip than near the wing root 

The term suction used in the discussion of the clean and iced wing pressures refers to a region of 
flow where the static pressure is lower than in the free steam. Suction is associated with 
negative pressure coefficients and is usually observed in regions of the wing where the 
combination of wing curvature and a accelerates the flow to speeds higher than the free stream 

In most cases, the pressure data presented in figures 3-52 to 3-69 indicated that the addition of 
icc shapes resulted in a dramatic change in the pressure distribution mainly due to a combination 
of leading-edge separation bubbles and extensive flow separation over the wing. Flow 
separation is a function of ice shape and a. The length of the separated region increased as a, or 
the height of the ice shape horn, were increased. The large changes in surface pressures 
downstream of the ice shapes were responsible for the observed behavior of the iced wing force 
and moment coefficients. 

Separated flows contain vortex structures, bubbles, and free shear layers that have a significant 
impact on the flow field. Vortex structures may be fixed in location as in the case of cavity 
flows and large stationary bubbles, or they can convect with the flow and interact with other 
vortices, free shear layers, or wall-bounded flows. Vortices can also stretch, depending on 
external forces (normal and shear forces) imposed on the vortex. The shape and strength of 
vortices in close prex unity to the airfoil surface can produce considerable changes in surface 
pressures and, therefore, in aerodynamic loads. 

The flow over iced wings, particularly for cases involving glaze ice shapes with large horns, 
features large regions of flow separation downstream of the icc shapes, which are typically 
referred to as leading-edge bubbles or long bubbles. These bubbles are regions of viscous flow 
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bound by the inviscid flow streamlines. The streamlines detach from the wing surface upstream 
of the bubble and reattach downstream of the bubble. If reattachment takes place over the wing 
surface, then the bubble is closed. The location of the reattachment point is a function of a and 
horn size. As a is increased, the bubble becomes progressively longer and eventually massive 
flow separation takes place once the bubble reaches the wing TE. In some cases, a long bubble 
can combine with separated flow initiated at the TE, and massive flow separation can occur 
before the bubble reaches the wing TE. 

For swept wings, LE separation bubbles may include a combination of LE flow separation and 
LE vortex flow. When sweep is incorporated in a wing, a conical vortex lying on the wing 
surface can occur, as discussed in reference 18. This LE vortex results from both the LE 
separation bubble and the spanwise pressure gradient due to sweep. The vortex cross section 
grows nearly perpendicular to the wing LE (conical vortex) in the spanwise direction. The 
growth of vortex diameter as the wing tip is approached is the result of a combination of vortex 
flow with the thick boundary layer transported by the spanwise flow from the inboard sections. 
With increased a, the vortex becomes stronger near the inboard sections but diffuses near the 
outboard sections. "Near the wing tip, the diffused conical LE vortex can combine with the wing 
tip vortex, and complex flow separation patterns are often t ^served. LE vortices were ohserved 
with both the clean and iced wing configurations tested as a was increased. 

In general, pressure distributions for wing flows with LE bubbles exhibit a considerable 
reduction in LE suction and a rounding of the LE pressure distribution. As the a is increased, 
the bubble becomes longer in the streamwise direction and the suction over the bubble is 
typically reduced. This led to a flat pressure distribution over the wing surface. In general, LE 
bubbles shift the load towards the LE. However, as the bubble becomes longer and massive flow 
separation takes place, the wing load shifts in the downstream direction. Massive flow 
separation is an unsteady phenomenon that involves convection of vortices, which cause 
significant changes to the load distribution over the wing arid, thus, to wing pitching moments. 
Flow s, 7 oration near the outboard portion of the iced wings could also cause large changes in the 
load over aileron controls, and in some cases, it could reduce or even eliminate the control 
effectiveness. Flow separation over the aileron is also responsible for control force reversal and 
potential large changes in control stick forces. Thus, bubble formation, growth, and bursting can 
cause large changes in the wing pitching moments and in aileron hinge moments. 

Another control issue regarding rolling moments is that ice accretions on the left and right wings 
arc not typically symmetric. As a result of the ice asymmetry, flow separation in the spanwise 
direction could be different between the two sides of the wing, thus leading to considerable 
rolling mnments- 

In reviewing the pressures distributions of the clean and iced wing configurations presented in 
figures 3-52 to 3-69, the following trends and features in the pressure data should be examined. 

• Location of the area centroid of the chordwise pressure distribution with respect to the 
25% MAC point of the wing. This determines pitching-moment behavior. 
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• Change in pressure distribution in the spanwise direction. This affects pan wise load 
distribution and provides insight into spanwise flow separation. The soanwise load 
distribution also affects wing rolling moment. 

• The growth of separation bubbles with a. Sudden changes in bubble features are 
undesirable because they cause sudden changes in wing loads. 


• The shape of the pressure distribution over the upper and lower control surface. This is 
responsible for the observed behavior of control Hin g e moments. 

3 .2.4.9 Clean Wins ( Figures 3-52 to 3-541 . 

• a = -4°: High suction was observed over the lower surface of the wing. Suction 
increased with spanwise distance and reached a maximum suction Cp of -3.53 at the 85% 
semispan station. The pressure data indicated that the flow was attached. 

• a = 0°: The flow remained attached. Due to the wing geometric twist, the average a was 
negative, resulting in negative lift. In all cases, the lower surface of the wing experienced 
suction that was lower near the root and higher near the tip. 

• a = 4°: For this a, aioc«i was 4° at the root and 0° at the tip. The pressure over the wing 
upper surface was lower than over the lower surface, resulting in positive lift. Upper 
surface suction was hi gher near the wing root than near the wing tip. The flow remained 
atta c h ed 

• a = 8°: Hign suction was observed over the wing upper surface at all spanwise stations. 
The suction Cp increased from -3.3 at 15% semispan to -3.5 at 50% semispan, and then it 
decre&wd t? -2.7 at 85% semispan. The pressure data indicated that flow separation was 
not present over the wing for this a. 

• a = 12°: LE separation bubbles are evident for this near stall a. The extent of the 
separation bubbles increased from root to tip, while the suction over the bubble 
decreases. 

• a = 16°: A large LE separation bubble was observed near the wing root. The flat Cp 
curves corresponding to the wing upper surface at the 50% and 85% sen span stations 
are indicative of massive flow separation. The pressure data presented for a of 12° and 
16° show that flow separation (wing stall) was initiated ut the tip of the wing and 
progressed toward the root. 

3.2.4.10 Iced Wing . 

Pressure distributions tor the iced wing with the IKT castings and the simulated smooth and 

rough LEWICE ice shapes are discussed below. The experimental pressure distributions for the 

iced configurations tested exhibited leading-edge bubbles over the upper, lower, and in many 

cases, over both surfaces of the wing. As the a was increased, the bukble(s) grew Longer in the 
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downstream direction and eventually burst into massive flow separation as the streamlines failed 

to reattach to the wing surface. 

3.2.4.10.1 Icing Condition 1— Ice shapes: IRT-CS10. LS-CS 10. LR-CS 10 (Figures 3-52 to 3-541 . 

• a = -4°: At 15% semispan, LE bubbles were evident over the wing lower surface for all 
ice shape cases. The suction over the bubble region was greater for the LEW1CE ice 
shapes compared to the IRT ice shape casting. At the mid-semispan location, extensive 
flow separation was observed over the lower wing surface for the t-ses presented. The 
flow remained separated over the lower surface near the tip station, and the suction was 
red' ’ced. .An LE bubble was present over the lower surface for the LEW1CE ice shapes. 

• a = 0°: Flow separation in the form of LE bubbles w as present over both wing surfaces 
for all ice shapes tested. Suction over the lower surface was higher tnan over the upper 
surface. In most cases, the LEWICE ice shapes resulted in higher suction over the lower 
surface. The bubble extent was less near the wing root than at the outboard stations. 
Near the wing tip, flow separation extended to the wing TE. 

• a = 4°: For this angle of attack, suction was increased over the wing upper surface. The 
Cp on the lower surface varied from high negative near the LE, indicating the presence of 
separation bubbles to low negative and, in some cases, slightly positive near the TE. The 
extent of the upper surface bubbles decreased (separation was reduced) with spanwise 
distance from root to tip due to the negative wing twist, which lowered the aiocai over the 
outbc rd sections. 

• a = 8°: Long bubbles were observed over the inboard wing upper surface for all the ice 
shapes. The bubbles grew longer (chotdwise) over the outboard sections of the wing and 
eventually, complete flow separation was observed near the wing tip. Suction over the 
bubbles for the LEWICE ice shapes was higher than obtained with the IRT ice shape 
casting. In most cases, the pressure distribution over the wing lower surface obtained 
with the LEWICE shapes exhibited bubble formation between 0% and approximately 
50% chord. 

• a = 12°: Large bubbles were observed at the 15% station, and complete flow separation 
was evident over the wing upper surface at the 50% and 85% semispan locations with all 
ice shapes tested. The IRT ice shape had lower suction over the two inboard stations than 
the LEWICE ice shapes. Bubbles were also present over the lower surface of the wing 
tor the two LEWICE shapes. 

• a = 16°: At this angle of attack, extensive flow separation was observed over the upper 
surface of the wing at all spanwise stations. The only exception was the smooth 
LEWICE ice shapes for which the flow over the inboard (15% semispan) station was 
partially attached. 
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3.2.4.10.2 Icing Condition 2— Ice Shapes: 1RT-IS10. LS-1S10. LR-IS10 (Figures 3-55 to 3-571 . 


• a = -4°: At 15% semispan, LE bubbles were evident over the wing lower surface for all 
ice shape cases. In addition, at the mid-semispan location, the bubbles grew longer. 
Furthermore, reduced suction and extensive flow separation were observed over the wing 
lower surface near the tip station. Note that the LE bubbles of the IRT-IS10 and LR-1S10 
cases were comparable at the 50% and 85% semispan stations, while the bubble of the 
LS-IS10 configuration produced higher suction. 

• a = 0°: Flow separation in the form of LE bubbles was present over both wing surfaces 
for all ice shapes tested. Suction over the lower surface was higher than over the upper 
surface. The IRT and LEWICE ice shapes resulted in comparable suction over the lower 
surface. The bubble extent was less near the wing root than at the outboard stations. 
Near the wing tip, flow separation extended to the wing TE. 

• a = 4°: For this angle of attack, suction was increased over the wing upper surface. The 
Cp on the lower surface varied from high negative near the LE, indicating the presence of 
separation bubbles, to low negative and, in some cases, slightly positive near the TE. 
Note that at the 15% semispan station, the suction on the wing upper surface, due to the 
IRT-1S10 ice shape, was lower than that produced by the LEWICE ice shapes. However, 
the difference in suction between the IRT and LEWICE ice shapes diminished with 
span wise distance from root to tip. 

• a = 8°: Long bubbles were observed over the inboard wing upper surface for all ice 
shapes. The bubbles grew longer (chordwise) over the outboard sections of the wing and 
eventually complete flow separation was observed near the wing tip. Suction over the 
bubbles for the LEWICE ice shapes was higher than that obtained with the IRT ice shape 
casting near the wing root. In most cases, the pressure distribution over the wing lower 
surface obtained with the LEWICE ice shapes exhibited bubble formation between 0% 
and approximately 50% chord. 

• At the higher a of 12° and 16°, extensive flow separation was observed over the upper 
surface of the wing at all spanwise stations. The only exception was the flow at the 15% 
semispan station, which was partially attached for all ice shape configurations. 

3.2.4,10.3 Icing Condition 3— Ice Shapes: IRT-SC5. LS-SC5. LR-SC5 (Figures 3-58 to 3-601 . 

• a = -4°: At 15% semispan, small LE bubbles were observed over the wing lower surface 
for all ice shape cases. The suction over the bubble region was comparable between the 
LEWICE and IRT ice shapes casting. The bubbles grew longer (chordwise) at the mid- 
semispan station. At the 85% semispan location, extensive flow separation was observed 
over the wing lower surface for all cases presented. An LE bubble was present over the 
lower surface for all ice shapes. 

• a = 0°: Flow separation in the form of LE bubbles was present over both wing surfaces 
for all ice shapes tested. Suction over the lower surface was higher than over the upper 
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surfc'c. The bubble extent was less near the wing root than at the outboard stations. 
Near the wing tip, flow separation extended to the wing trailing edge. 

• For angles of attack of 4° and 8°, Cp distributions of all ice shape configurations were in 
good agreement with the clean case. At the two inboard locations, the minimum pressure 
near the LE for the iced wing was lower than the clean wing, indicating lower a «. The 
reduction in a« was due to the LE droop caused by the rime ice accretion. 

• a = 12°: At the 15% semispan station, only the IRT-SC5 ice shape produced a LE 
bubble on the wing upper surface. At the same time, observe that suction was increased 
for both LEWICE ice shapes. As the flow progressed to the 50% semispan location, the 
Cp distributions exhibited a bubble near the LE for all three iced wing configurations, 
where the suction over the bubble for the LS-SC5 case was higher than the other two SC5 
ice shapes. The LE bubbles grew in chordwise length from the 50% semispan station to 
the wing tip. 


• a = 16°: At this angle of attack, all ice shapes maintained partial attached flow over the 
wing suction surface at the two inboard stations. However, extensive flow separation 
was observed over the upper surface of the wing at the 85% location. 

3.2.4. 10.4 Icing Condition 4 — Ice Shapes: 1RT-CS2. LS-CS2. LR-CS2 (Figures 3-61 to 3-631 . 

• a = -4°: At 15% semispan, LE bubbles were evident over the wing lower surface for all 
ice shape cases. From root to mid-semispan stations, the bubble region extent increased. 
Extensive flow separation was observed over the wing lower surface near the tip station 
for all cases presented. 

e a = 0°: Flow separation in the form of LE bubbles was present over both wing surfaces 
for all ice shapes tested. Suction over the lower surface was higher than over the upper 
surface. Note that the bubble extent was lesr near the wing root than at the outboard 
stations. However, the LEWICE ice shapes produced separation bubbles of greater 
extent than the IRT ice shape. 

• a = 4°: For this angle of attack, suction was increased over the wing upper surface. 
Separation bubbles were observed over the lower surface between 0% and approximately 
20% chord. The extent of the upper surface bubbles decreased (separation was reduced) 
with span wise distance from root to tip due to the negative wing twist, which lowered the 
(Xiocai over the outboard sections. Moreover, suction over the LE bubbles of the LEWICE 
ice shapes was higher than the 1RT-CS2 case at all spanwise stations. 

• a = 8°: Once again, bubbles were observed over the inboard wing upper surface for all 
the ice shapes. The bubbles grew longer (chordwise) over the outboard sections of the 
wing. For the LEWICE ice shapes, suction over the bubbles was higher than obtained 
with the IRT ice shape casting at the 15% semispan station. However, at the mid- 
semispan and near tip stations, the smooth and rough LEWICE ice shapes had lower 
suction over the separation bubbles than the IRT-CS2 ice shape. 
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• a = 12°: Large LE bubbles were observed at the 15% semispan station, and complete 
flow separation was evident over the wing upper surface at the two outboard locations 
with all the ice shapes, except the 1RT-CS2 ice shape for which the flow was partially 
attached near the wing tip. 

• a = 16°: At this angle of attack, extensive flow separation was observed over the upper 
surface of the wing at all span wise stations. 

3.2.4. 10,5 Icing Condition 5— Ice Shanes: 1RT-CS22. LS-CS22N. LS-CS22S. LR-CS22N. and 

LR-CS22S (Figures 3-64 to 3-661 . 

• a - -4°: At 15% semispan, LE bubbles were evident over the wing lower surface for all 
ice shapes. Suction over the bubble region was greater for the LEW1CE ice shapes 
compared to the 1RT-CS22 case. Comparing the smooth and rough LEWICE ice shapes, 
it was observed that the addition of roughness to the ice horns produced bubbles with 
lower suction but of greater chordwise extent than obtained with the smooth LEWICE ice 
shape. At the mid-semispan location, extensive flow separation was observed over the 
wing lower surface for all cases, except the LS-CS22S and LR-CS22S ice shapes. The 
flow remained separated over the lower surface near the tip station. 

• a = 0°: Flow separation in the form oi LE bubbles was present over both wing surfaces 
for all ice shapes tested. Suction over the lower surfac vas higher than over the upper 
surface. In most cases, the LEWICE ice shapes resulted in hi ghe r suction over the lower 
surface. Near the wing tip, flow was completely separated. 

• a = 4°: For this angle of attack, LE bubbles were present over both surfaces near the 
wing root. For the rough LEWICE ice shapes, suction over the bubbles was higher than 
their smooth counterparts. At the mid-semispan location, extensive flow separation was 
observed for the IRT and rough LEWICE ice shapes, while for the LS-CS22N and LS- 
CS22S cases, the flow was partially attached. An LE bubble was evident for the IRT- 
CS22 ice shape at the 85% semispan station. Flow remained separated near the wing tip 
for all LEWICE cases 

• a = 8°: Long bubbles were observed over the inboard wing upper surface for all 

LEWICE ice shapes, whereas flow downstream of the IRT-CS22 horns was completely 
separated over the wing upper surface. For the LEWICE cases, the LE separation 
bubbles failed to reattach over the wing sections at the 50% and 85% semispan stations, 
and complete flow separation was observed. For the wing with LEWICE ice shapes, the 
pressure distribution over the lower surface exhibited bubble formation between 0% and 
approximately 50% chord. 

• a = 12°: Aft of 1RT-CS22 upper ice horn, complete flow separation occurred. For all 
LEWICE ice shapes, large bubbles were observed at the 15% semispan station, and 
massive flow separation was evident over the wing upper surface at the 50% and 85% 
semispan locations. Note that flu. IRT ice shape had lower suction than the LEWICE icc 
shapes. Bubbles were also present over the lower surface of the wing for all ice shapes 
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from root to tip. Yet, only the LS-CS22S ice shape maintained a lower surface bubble 
near the wing tip. 

• a = 16°: At this angle of attack, extensive flow separation was observed over the upper 
surface of the wing at all span wise stations. Once again, all ice shapes produced an LE 
bubble over the wing pressure surface at the three span wise stations, as shown in the 
figures. 

3.2.4.10,6 Icing Condition 6— Ice Shanes: IRT-IPSF22. LS-1PSF22. LR-IPSF22 t Figures 3-67 

to 3-691 . 

• a = -4°: At 15% semispan, LE bubbles were evident over the wing lower surface for all 
ice shape cases. The suction over the bubble region was greater for the LEWICE ice 
shapes compared to the IRT ice shape casting. At the mid-semispan location, bubbles of 
greater extent were observed over the wing lower surface for all cases presented. Over 
the lower surface, the flow was separated near the tip station for the IRT-1PSF22 
configuration, whereas an LE bubble was present for the LEWICE ice shapes. 

• a = 0°: Flow separation in the form of LE bubbles was present over both wing surfaces 
for all ice shapes tested. Suction over the lower surface was higher than that over the 
upper surface. In most cases, the LEWICE ice shapes result id in higher suction over the 
lower surface. The bubble extent was less near the wing root than at the outboard 
stations. Near the wing tip, flow separation extended to the wing TE. 

• a = 4°: For this angle of attack, suction was increased over the wing upper surface. The 
Cp on the lower surface varied from high negative near the LE, indicating the presence of 
separation bubbles, to low negative and, in some cases, slightly positive near the TE. 

• a = 8°: Long bubbles were observed over the inboard wing upper surface for all ice 
shapes. The bubbles grew longer (chordwise) over the outboard sections of the wing, and 
eventually, complete flow' separation was observed near the wing tip. Suction over the 
bubbles for the LEWICE ice shapes was higher than obtained with the IRT ice shape 
casting. Separation bubbles were also observed on the lower wing surface for all ice 
shapes. The lower surface bubble increased its chordwise length from root to tip. 

• a = 12°: Large bubbles were observed at the 15% station, and complete flow separation 
was evident over the wing upper surface at the 50% and 85% semispan locations for all 
ice shapes tested. The IRT ice shape had lower suction over the two inboard stations than 
the LEWICE ice shapes. Bubbles were also present over the lower surface of the wing 
for all ice shapes. 

• a = 16°: At this angle of attack, extensive flow separation was observed over the upper 
surface of the wing at all spanwisc stations. The only exception was the smooth 
LEWICE ice shapes for which the flow over the inboard (15% semispan) station was 
partially attached. Once again, LE bubbles were observed for all ice shapes over the 
wing lower surface. 
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3.2.5 Effect of Aileron Deflection . 

The effect of aileron ceflection on the aerodynamic performance of the clean and iced wing 
is presented in figures 3-70 to 3-73 and in table 3-10. The main performance parameters 
presented include: 

• Cl versus a for all aileron deflections (5a) 

• Ch versus a for all 5a 

• Cl versus 5a for selected a 

• Ch versus 5 a for selected a 

Another important parameter for assessing aileron performance is the change in rolling moment 
with Sa and a. Rolling moment data were obtained during the experimental investigation at 
WSU. 

The experimental lift and hinge-moment data presented in figures 3-70 to 3-73 were obtained 
with a reflection plane wing model having a single aileron. Aircraft wings, however, have left 
and right ailerons that are not independent but are connected via cables to the control wheel. 
Thus, the control forces experienced by the pilot are due to the difference in the hinge moments 
generated by the left and right ailerons. To estimate the net aileron control force from the 
experimental data presented, the hinge moment for positive and negative aileron deflections must 
be combined. For example, using the results presented in figure 3-70(f) for the clean wing, a of 
0° and 6 a of 5°, the difference between for aileron deflection of -5° and Ch for aileron 
deflection of 5° provides the net aileron control forc e. This procedure assumes that the up and 
down aileron deflection angles are the same. However, in some aircraft, left and right aileron 
deflection angles are not the same, and in such a case, the actual deflections for each aileron have 
to be used to compute the control force. Another assumption made regarding the 
superimposition of the experimental results is that ice accretions for the left and right wings are 
symmetric. In general, the left and right wing ice shapes will not be exactly the same, and this 
could cause different separation patterns over die left and right sides of the wing and the ailerons. 
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FIGURE 3-70. EFFECT OF AILERON DEFLECTiON ON LIFT AND AILERON HINGE 
MOMENT COEFFICIENTS; CLEAN CONFIGURATION; Re= 1.8xl0 6 


3-83 












Lift coefficient Hng* moment coefficient 





(c) I-., ige-moment coefficient (negative 5 a) 



006 




-020 | t I I 1 p I M p I I 1 p t I T p 1 1 I p t I 1 p M I ) 

-10 -6 0 5 10 15 20 25 

Angle of attack (deg) 


(d) Hinge-morr.ent coefficient (positive 6 a) 



-16 -10 -6 0 5 10 16 20 “ t 6 -10 -5 0 5 10 16 20 

Altar on Deflection (dag; AU«ron Diction (deg) 

(e) Lift coefficient vs 6a (f) Hinge-moment coefficient vs 6a 

FIGURE 3-71, EFFECT OF AILERON DEFLECTION ON LIFT AND AILERON HINGE- 
MOMENT COEFFICIENTS, 1RT-CS22 ICE SHAPE; Re= 1.8x10° 
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(a) Lift coefficient (negative 6a) 
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(e) Lift coefficient vs 6a 
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FIGURE 3-72. EFFECT OF AILERON DEFLECTION ON LIFT AND AILERON H1NGJ 
MOMENT COEFFICIENTS; 1RT-IPSF22 ICE SHAPE; Re= 1.8xl0 6 
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(c) Hinge-moment c ^efficient (negative 5a) (d) Hinge-moment coefficient (positive 5a) 




(e) Lift coefficient vs 5 a (?) Hinge-moment coefficient vs 6 a 

FIGURE 3-73. EFFECT OF AILERON DEFLECTION ON LIFT AND AILERON HINGE- 
MOMENT COEFFICIENTS; LR-1PSF22 ICE SHAPE; Re= 1.8xl0 6 
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TABLE 3-10. C L ,suii OF ALL CONFIGURATIONS WITH 
AILERON DEFLECTION; R6uac= 1-8x10° 


^Conll^ 

-15 e 

-10° 

-5° 

E9 


2.5° 

5° 

10° 

15° 

fO 

o 

u 

Clean 

0.829 

0.845 

0.858 

0.868 

0.874 

0.883 

0.887 

0.904 

0.913 

0.929 

LRT-CS10 

0.491 

0.506 

0.522 


0.536 

0.544 

0.533 

0.573 

0.593 

0.613 

IRT-IS10 

0.591 

0.603 

0.618 

0.628 

0.646 

0.648 

0.662 

0.t76 

0.688 

0.716 

IRT-SC5 

0.862 

0.874 

0.888 

0.893 

0.896 

0.908 

0.918 

0.930 

0.948 

0.961 

IRT-CS2 

0.737 

0.752 

0.770 


0.775 


0.788 

0.807 


0 83 4 

IRT-CS22 

0.005 

0.016 

0.025 


0.033 

0.036 

0.037 

0.039 

0.039 

o.o: 

1RT-IPSF22 

0.506 

0.515 

0.527 


0.528 

0.527 

0.531 

0.546 

0.563 

0.579 

LS-CS10 

0.669 

0.679 

0.692 

0.705 

0.710 

0.719 

0.723 

0.744 

0.760 

0.772 

LS-IS10 

0.640 



0.672 

0.679 

0.687 

0.697 

0.722 

0.735 

0.747 

LS-SC5 

0.900 

0.915 

0.927 

0.933 

0.939 

K&m 

0.965 

0.980 

0.984 

1.000 

LS-CS2 

0.714 

0.722 

0.737 

0.741 

0.748 

0.759 

0.769 

0.780 

0.800 

0.812 

LS-CS22N 

0.450 1 


0.476 

0.476 

0.484 

0.490 

0.496 

0.505 


0.520 

LS-CS22S 

0.495 ! 

0.507 

0.516 

0.517 


0.542 

0.548 


0.575 

0.601 

LS-LPSF22 

0.633 

0.648 

0.660 

0.667 

0.676 

0.683 

0.697 

0.713 

0.737 

0.746 

■Q2SD3H 

0.564 

0.584 

0.597 

0.607 

0.613 

O.o2 4 

0.635 

0.647 

0.669 

0.683 

LR-IS10 

0.589 

0.605 



0.644 

0.651 

0.661 

0.678 

0.694 

0.711 

LR-SC5 

0.921 

0.939 

0.955 

0.962 

0.980 

0.982 

0.994 

1.012 

1.027 

1.044 

LR-CS2 

0.677 

0.691 

0.703 

0.709^ 

0.717 

0.730 

0.736 

0.754 

0.770 

0.783 

LR-CS22N 

0.406 

0.405 

0.409 

0.406 

0.407 

i»glU 

0.413 

0.416 

0.424 1 

0.427 

LR-CS22S 


0.351 

0.349 

0.346 

0.342 " 

0.344 

0.343 

0.352 

0.351 

0.356 


0.566 

0^77 

0.590 

0.602 

0.612 

0.617 

0.626 

0.646 

0.659 

0.680 


3.2.5. 1 General Comments Regarding Aileron Controls . 

Ailerons are designed to meet aircraft roll performance criteria in terms of roll rate (typically 1 
to 15 degrees per second) and attitude. In addition, the ailerons should be able to provit 
sufficient roll control to hand!' special aircraft situations such as an engine-out case or steal 
heading sideslip. Aileron forces should remain within the capabilities of the pilot as defined 1 
the certification requirements [11], which state that the maximum permissible aileron contr 
force is 50 lb with two hands on the control wheel and 25 lb for a single hand on the contr 
wheel (autopilot settings vary with design but often a 20 lb limit is used). 

The change in aileron hinge moments with 5 a and a has a direct impact on aileron control force 
The slope of Ch with 8a, Ch.sa, also known as the control heaviness parameter, is of prima 
importance in evaluating control surface behavior and control forces. In general, aileron desig 
have negative Ch,sa so that the aileron has the tendency to return to its undeflected position. T1 
hinge-moment slope, Ch,<i, also known as the control floating parameter, affects the change in tl 
control force resulting from the response of the aircraft to the control movement. W hen tl 
ailerons are deflected and wing roll is initiated, the local a of the upward and downward movii 
wings is changed (typically by a couple of degrees for small deflection), and the control for 
required to maintain a steady maneuver is cither greater or less than the control force required 
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iniiiatp the maneuver, depending on the sign and magnitude of Cra- Control surface designers 
us uall y aim at small positive or negative values for Ch, q to avoid large differences in control 
forces during maneuver [18 and 19]. The effect of Crq on control forces is of secondary 
importance, particularly for small aileron deflection and roll rates that cause small changes in a. 
For large deflections, however, which can produce large roll rales and therefore a large change in 
a, C h „ could have a notable effect on control forces. 

Once the clean wing aileron design requirements have been satisfied, the effect of ice shapes on 
aileron performance should be assessed for unwanted changes in control behavior. In general, 
the presence of ice on the wing leading edge can result in early flow separation over the wing 
and control surfaces, and the iced wing lift with 5 a could change with respect to the clean wing. 
It is important that the lift of the iced wing exhibits the expected behavior with positive and 
negative aileron deflections, although some loss in aileron effectiveness is expected. 

The effect of ice shapes on aileron hinge moments (i.e., control forces) and on the overall 
behavior of the hinge-moment curve should be evaluated. With iced wings, the growth in hinge 
moment with a (region B as defined in figures 3 -26(a) and 3-26(b)) takes place at a lower a with 
respect to the clean wing, as shown in figures 3-26(c) and 3-26(d). For such a case, a large 
amplitude roll maneuver could change the a sufficiently to place the aileron on the down moving 
side of the wing in the steep region of the hinge-moment curve. This could cause a considerable 
increase in control force. If, in addition to the increase in control force, the aileron effectiveness 
is reduced because of flow separation over the ailerons caused by the ice shape, then a situation 
could develop where recovery may be difficult. 

in uc«essing aileron performance for iced wings, the magnitude of Ch and the slope of Ch with 
respect to a and 5 a should be considered. Ideally, the magnitude of Ch should be balanced; 
.small enough to avoid excessively hi gh forces, yei large enough to avoid very small forces. The 
actual permissible values of Ch will be a function of the aerodynamic and mechanical design of 
the control, as well as certification requirements. Any change in the slope of Ch with a or 5 a 
should al-sn be small and gradual. 

3.2. 5,2 Clean Wing . 

For a reflection place model, positive aileron deflections increase wing lift and reduce the a.«..n. 
while negative deflections have the opposite effect. Note that large positive aileron deflections 
at high angles of attack should be evaluated carefully since they increase the an gl e of attack over 
the outboard portion of the wing and can result in early tip stall. 

For the clean wing tested, the aileron remained effective for all a and aileron deflections, i.e., 
positive aileron deflection (TE down) increased wing lift, while negative deflection resulted in 
lift reduction, as expected. In addition, as the aileron deflection was increased, the lift increment 
(for aileron TE down) or decrement (for aileron TE up) increased, as shown in figures 3-70(a) 
and 3-70(b). 

The aileron hinge moments presented in figures 3-70(c) and 3-70(d) exhibited gradual change 
with a for all aileron deflections tested. The slope, Ch^i, for all three regions A, B, and C (as 
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defined in figure 3-26) of the Ch curve were negative, indicating that the hinge moment became 
less positive or more negative as a was increased from -8° to 20°. On a number of cases, a 
change in the sign of Ch was observ ed in the hinge-moment data. Consider, for example, the 
case in figure 3-70(c), corresponding to 8 a of -15° (TE up). For this case, a change in a from 
10° to 16° resulted in a change in aileron hinge moments from 0.07 to approximately -0.03. The 
change in sign of Ch occurred near wing stall, corresponding to a of approximately 13.5°. For a 
between 10° and 13.5° and for 8 a of -15°, the hinge moment was positive, indicating that the 
aileron trailing edge had the tendency to return to its undeflected position. However, at a greater 
than 13.5°, the aileron had the tendency to maintain its deflected position due to extensive flow 
separation over the upper surface of the aileron. Flow conditions which force a control surface 
to maintain its deflection status result in control force reversal. However, with aileron controls, 
such a condition over one of the ailerons may not be a problem if the other aileron (the one 
deflected TE down in this case) remains effective. Considering the 15° aileron deflection case in 
figure 3-70(d), it is observed that for a of 13.5°, the hinge moment for the TE down aileron was 
approximately -0.18, i.e., the TE down aileron had a strong tendency to return to its neutral 
position. Thus, the combined hin g e moment and, therefore, the net control force due to both 
ailerons maintained the correct behavior (i.e., the controls had the tendency' to return to the 
neutral position). 

The change in Cl with 8 a depicted in figure 3-70(c) was practically linear for a in the range of - 
8° to 16°. The change in lift as S A was increased from -15° to 20°, while maintaining a constant 
a. was in the range of 0. 1 to 0.2, depending on a. The higher lift increments were obtained for a 
in the linear lift range where the flow over the control surface was mostly attached. 

The change in control hinge moment with 8a is demonstrated in figure 3-70(0- The aileron Ch 
was a function of the prevailing flow conditions over the aileron. Consider the line in 
figure 3-70(f) corresponding to a of 0°. At this a, the flow was attached on both surfaces of the 
wing for all aileron deflections. The results show that for negative 8a (aiieron TE up), Ch was 
positive, that is the aileron had the tendency to return to its undeflected position. For positive 
aileron deflections (TE down), Ch is negative, which once again indicates that the aileron had the 
tendency to return to the neutral position. The linear behavior of the curve shows that as 8 a was 
increased, the moment and, therefore, the control force were increased in a linear fashion. Note 
that for a full-span wing, the net aileron control force can be obtained from the results presented 
in figure 3-70(0 by computing the change in Ch, corresponding to the same positive and negative 
aileron angular deflections. For example, at a of 0°, the change in the magnitude of Ch 
corresponding to aileron deflections of ±5°, was approximately 0.075; Ch for TE up aileron was 
0.032, while Ch for TE down aileron was -0.043. 

Now consider the Ch versus 8 a curve corresponding to a of -8°. At this a, the upper surface of 
the wing was the pressure side, while the lower surface was the suction side. For negative 
aileron deflections, the hinge moment was positive as the aileron upper surface experienced 
higher pressure than the lower surface. For small positive aileron deflections, the pressure over 
the aileron upper surface remained higher than over lower surface, which now had lower suction, 
and the hinge moment remained positive. As the aileron approached maximum positive 
deflection, the pressure over the upper surface did not change significantly but the suction over 
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the lower surface w as considerably reduced. Thus, the magnitude of the positive hinge moment 
was reduced. The reason for the nonlinearity in Cu versus 5 a curve was mild flow separation 
over the aileron lower surface. 

The results presented in figure 3-70(f) demonstrate approximately linear change in Ch with 5a 
for a corresponding to the linear lift range. For these a, the change in aileron hinge moment was 
approximately 0.21 as 8 A was increased from -15° to 15°. At high negative or positive a, the 
change in hinge moment with aileron deflecti.cn was not linear. This is demonstrated in figure 
3.70(f) for a of 16°, 12°, and -8°. 

3.2.5.3 1RT-CS22 Ice Shape . 

The effects of the 22.5-min glaze ice accretion with the large scallop features on the aerodynamic 
performance of the swept wing for all aileron deflections tested is presented in figure 3-71. The 
experimental results obtained indicate the following: 

• Lift increase or decrease with 5 a was as expected for most a as demonstrated in 
figures 3-71 (a) and 3-7 1(b). For a-range of 1° to 6°, however, negative aileron 
deflections increased lift while positive deflections had practically do effect on lift. This 
was mainly due to massive flow separation over the wing and control surfaces caus'd by 
the large size (2- to 3-in.) ice horns. 

• The maximum positive Ch was 0.1 and occurred for 5a of -15° and a of -8° as shown ir 
figure 3-71(c). The maximum negative hinge-moment coefficient was -0.19 at a of 20° 
and 5 a of 20° as evident from figure 3-7 1(d). Corresponding maximum and min imum 
values for the clean wing were 0. 1 65 and -0.22. Thus, the iced wing hinge moments were 
within the maximum and minimum limits of the clean wing hinge moments. It is 
important, however, to point out that for some aircraft, the maximum hinge moment for 
the clean wing are typically not reached either due to the aircraft natural angle of attack 
limits or due to an artificial stall barrier. Id such a case, the iced wing hinge moments 
should be compared to the actual clean wing hinge moment limits. 

• A change in the sign of Ch was observed for all negative aileron deflections and for 
a-range of 2° to 5°, as shown in figure 3-71(c). For positive aileron deflections in the 
range of 0° to 5°, the change in the sign of Ch occurred between a of 1° and 3°. as 
demonstrated in figure 3-71(d). For larger positive aileron deflections, the sign of Ch 
was reversed twice near a of 0°. As discussed previously, for the clean wing case, a sign 
change in Ch indicates that the aileron has the tendency to maintain its deflected position. 
To determine the net aileron control hinge moment and its direction for the combined left 
and right aileron system, the Ch for both positive and negative aileron deflections should 
be examined. 

• As demonstrated in figure 3-7 1(e), the change in wing lift with 5 a for constant a was 
practically linear for all a presented. The lines had a small positive slope for most a. 
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However, for a-range of 0° to 4°, the slope was nearly zero and, in some cases, slightly 
negative, indicating that the aileron was not effective in changing lift 

• The behavior of Ch with 8 a was notably nonlinear for ail a presented in figure 3-71(f). 
Given an a and ass uming attached flow, Ch will be positive for negative aileron 
deflections (TE up) since the aileron will tend to move to its undeflected position. For 
positive aileron deflections, the sign of Qi will eventually become negative. For fixed 
positive aileron deflection, the pressure over the aileron lower surface and the suction 
over the upper surface should increase as the a is increased from -8° to 0° (note that such 
a change in a with fixed 8 a is not practical under flight conditions). The Ch should 
become progressively more negative and the hinge-moment curves will move further 
apart. This was not the case with the 1RT-CS22 ice shape, as shown in figure 3-7 1(f) for 
a-range of -8° to 0°. The intersection in the hinge-moment curves corresponding to these 
a near 8a of 5° is indicative of extensive flow separation over the control surface. 

3 .2.5 .4 1RT-IPSF22 and LR-IPSF22 Ice Shapes . 

Aerodynamic performance with aileron deflection for the 22.5-min glaze ice shape casting IRT- 
IPSF22 and the corresponding rough LEWICE ice shape LR-IPSF22 are presented in 
figures 3-72 and 3-73 respectively . These ice accretions were representative of ice protection 
system failure cases for the wing model tested. The results presented in figures 3-72 and 3-73 
indicate the following: 

• Increase or decrease in lift with Sa was as expected for all aileron deflections and a, as 
demonstrated in figures 3-72(a), 3-72(b), 3-73(a), and 3-73(b). 

• The maximum positive Ch for both ice shapes was approximately 0.14 and occurred for 
8a of -15° and a of -8°, while the maximum negative Ch was -0.2 at a of 20° and 8a of 
20°, as shown in figures 3-72 (c), 3-72(d), 3-73(c), and 3-73(d). Once again, the iced 
wing hinge moments were within *he maximum and minimum hinge-moment bounds of 
the clean wing. 

• The change in lift with 8 a for the IRT-IPSF22 and LR-IPSF22 ice shapes was practically 
linear for all a as shown in figures 3-72(e) and 3-73(e). For all cases presented in these 
figures, the Cl versus 8a curves had a small positive slope indicating a monotonic 
increase in lift with aileron deflection. 

• For most of the cases presented in figures 3-72(0 and 3-73(0, the change in the hinge- 
moment coefficient with aileron deflection exhibited nonlinear behavior mainly due to 
flow separation over the aileron. 

3.2.6 Effect of Large Scallop Features . 

In general, ice accretions on swept wings can be classified as complete scallops, incomplete 
scallops, and no-scallop, as discussed in references 2, 4, 20, and 21. Schematics of these three 
types of ice accretions are provided in figure 3-7< . 
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Complete scallops are glaze ice accretions that appear only on swept wings and are characterized 
by the presence of scallop tips extending from the attachment line. The scallop tips have a 
particular shape, height, and spacing. Incomplete scallops are glaze ice shapes that also appear 
only on swept wings and where scallops tips form beginning at some distance from the 
attachment line. No-scallops are ice accretions where scallop tips are not present. 

Scallop tips are made of glaze ice feathers, which develop from roughness elements that form 
dtiring the ice accretion process as shown in figure 3-75(a). The feathers have a preferred 
direction of growth that is perpendicular to the external streamlines, as shown in figures 3-75(b) 
and 3-75(c). Ice feathers join along the preferred direction to form ridges, witn incipient scallop 
tips formed by the feathers at the end of each ridge. As the ridges grow, they form scallop tips. 
As the ice accretion grows, the scallop tips mergw with adjacent scallop tips by joining at the top. 
This increases both the size of the scallop tips and the gap betw een them. Detailed descriptions 
and experimental data for scallop ice accretions can be found in references 2, 4, 20, and 21. 

The type of ice accretion that will form on a swept wing depends on the airfoil geometry, sweep 
angle, flow, and tunnel conditions. Icing conditions that lead to scalloped ice features can 
produce ice shapes with complex surface characteristics such as the IRT-CS10 and IRT-CS22 
c*ses shown in figure 2-33. 

Ice shapes with large scallop formations contain gaps between scallops, as demonstrated in 
figure 2-33(b) for the IRT-CS22 case. Aerodynamicists have debated over the years the effects 
of scallop features on aerodynamic performance. This is of practical interest because state-of- 
the-art ice accretion codes are not able to produce scalloped ice shapes. Simulated ice shapes for 
swept wings developed with the use of ice accretion codes or other empirical means have solid 
horns, as shown in figure 2-34. If the scallop-gap features in the ice horns have a significant 
impact on aerodynamic performance, then methods should be explored for incorporating their 
effects in simulated ice shapes. 

To investigate the effect of scalloped ice shapes on aerodynamic performance, a limited study 
was conducted with the 1RT-CS22 ice shape. In this study, the gaps between scallops were 
progressively filled with a modeling compound to produce an ice shape with solid horns, as 
shown in figure 3-76. To reduce the aerodynamic loads on the modeling compound, the 
experimental study with the baseline and modified 1RT-CS22 ice shapes was conducted at a 
lower airspeed for which the Reynolds number based on MAC was 1 .0 million. Experimental 
results from these tests are presented in figure 3-77. The first set of tests was conducted with the 
middle 25-in segment of the 1RT-CS22 ice shape filled with the modeling compound (case 1, 
R516 in figure 3-77). Next, the gaps in the 18-in tip segment were filled (case 2, R517 in 
figure 3-77), and finally, the gaps in all three ice segments were filled (case 3, R518 in figure 3- 
77). Each time a new set of force, moment, and pressure data was obtained. 
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(a) Complete scallop 


(b) Incomplete scallop 
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(a) Attachment line zone, glaze ice feathers zone 
and critical distance (view from direction 
normal to leading edge) 
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(b) Attachment line zone, glaze ice feathers zone 
and critical distance (overall view) 





(c) No scallops 

FIGURE 3-74. ICE ACCRETION ON 
A SWEPT WING AT GLAZE ICE 
CONDITIONS [21] 


(c) Feather detail 

FIGURE 3-75. DEVELOPMENT OF 
ROUGHNESS ELEMENTS AND ICE 
FEATHERS DURING ICE ACCRETION 
ON A SWEPT WING [20 AND 21] 
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Lift, drag, hinge-moment, and pressure coefficients for the baseline 1RT-CS22 ice shape (R1 13 
in figure 3-77) and for the three cases where the gaps between scallops were partially or 
completely filled with modeling compound are compared in figure 3-77. The results presented 
in this figure show a progressive increase in Cl for a-range of -2° to 9° as more gaps in the ice 
horns were filled with the modeling compound. The maximum gain in Cl was obtained for 
case 3 and ranged from 0.022 to 0.073 with respect to that of the baseline 1RT-CS22 case. The 
change in Cd with respect to the baseline ice shape was small but measurable for all three cases 
tested. For case 1, the change in the Ch with respect to the baseline IRT-CS22 ice shape was 
small. However, for cases 2 and 3, the hinge moments were less negative (i.e., the aileron 
trailing edge had less of a tendency to move up). 

The exact mechanism for the observed aerodynamic effects caused by the gaps between the ice 
scallops is not known. However, the pressure data at 15%, 50%, and 85% semispan stations for 
case 3 (R518) presented in figure 3-77 for a of 5° offer some clues. The pressure distributions 
for the case where all the gaps between scallops were filled with the modeling compound 
exhibited increased suction on both the upper and lower wing surfaces with respect to the 
baseline 1RT-CS22 ice shape. 

For large glaze ice accretions, the region between the horns is typically a high-pressure region 
because the air flow slows down within the cavity formed by the horns. The regions downstream 
of the upper and lower boms are in general low-pressure regions with small or large separation 
bubbles. With scalloped ice shapes, the gaps between the scallops allow the high pressure 
between the horns to leak to the low-pressure region downstream of the horns, thus reducing the 
suction near the wing LE. This is clearly evident in the 15% and 50% semispan pressure 
distributions presented in figures 3-77(d) and 3-77(e). At the 3C% semispan station, the IRT- 
CS22 ice shape with modeling compound produced lower suction peaks near the wing LE than 
the baseline ice shape. This was mainly due to flow separation that was more extensive for the 
case of the modified ice shape. Since Cl depends on the pressure difference between the upper 
and lower wing surfaces, as this difference was reduced due to flow leakage in the case of the 
baseline 1RT-CS22 ice shape, the net lift was also decreased. 

The observed lift increase with the modified 1RT-CS22 ice shape (case 3, R518) in figure 3-77(a) 
can be explained by direct examination of the pressure distributions presented in figures 3-77(d) 
to 3-77(f). The pressure distributions show that, for the 15% semispan station, case 3 produced 
72% more positive lift than the baseline case. At the 50% semispan station, case 3 resulted in 
18% less negative lift than the baseline ice shape. At the 85% semispan location, both 
configurations resulted in about the seme amount of negative lift. Thus, for case 3 where all 
gaps between the scallops were closed, the increase in suction over the wing resulted in a net 
increase in Cl with respect to the baseline ice shape. 

In summary', the limited study conducted with the glaze 1RT-CS22 scalloped ice shape showed 
that the scallop features and in particular the gaps between scallops can result in greater loss of 
lift than ice shapes with solid horns. It must be stressed, however, that these results should not 
be generalized until further studies are conducted with more scalloped type ice shapes. 
Furthermore, the results presented here may only be applicable to ice shapes with large scalloped 
features and for ice shapes with small scallops, the eff ects may not be as significant. 
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3.2.7 Reynolds Number Effects . 


Reynolds number effects on the clean and selected iced wing configurations is presented in 
figure 3-78. Reynolds numbers included 0.5, 1.0, 1.5, and 1.8 million were computed based on 
the wing mean aerodynamic chord of 1.56 ft. Corresponding tunnel airspeeds were 42, 84, 128, 
and 152 mph. 

For the clean wing, experimental results for Reynolds numbers in the range of 0.5 to 1.8 million 
are presented in figure 3-78. For the low Reynolds number of 0.5 million, the Cl,suii was 
approximately 0.87 and a*uii was approximately 16°. For higher Reynolds numbers of 1.0 and 
1.5 million, both Cl^uu and dsuii decreased and were approximately 0.84 and 15° respectively. 
However, as the Reynolds number was increased to 1.8 million, C i rt .n increased to 0.87 while 
cu.n decreased to 13.8°. 

Reynolds number effects on C L of the 1RT-1PSF22, LS-IPSF22, LR-1PSF22, IRT-SC5, and LR- 
SC5 configurations art presented in figures 3-78b to 3-78(0. h* general, the Reynolds number 
effects on Cl of the IRT-IPSF22 ice shape w as small, as demonstrated in figure 3-78(b). The 
effects of Reynolds number on the nea. stall lift performance of the LS-IPSF22 and LR-IPSF22 
ice shapes were more pronounced than for the IRT-IPSF22 case, as demonstrated in figures 3- 
78(c) and 3-78(d). For the LR-IPSF22 configuration, as the Reynolds number was increased 
from 0.5 to 1 .8 million, C r ...n decreased while c u.n remained approximately the same. 

For the 5-min rime 1RT-SC5 and LR-SC5 ice shapes, the effects of Reynolds number on Cl are 
demonstrated in figures 3-78(e) and 3-78(f) respectively. Most of the change in lift with 
Reynolds numoor was observed in the near-stall and poststall flow regimes. In general, Ci^un 
decreased as the Reynolds number was increased from 1.0 to 1.8 million. The behavior of CL 
for the low Reynolds number case of 0.5 million was similar to that for the 1 .5 million Reynolds 
number. The Reynolds number effects observed with the 5-min rime ice cases were mainly due 
to the change in the location of flow separation over the rough ice shapes. At low Reynolds 
numbers, roughness can enhance or reduce performance, depending on the value of Reynolds 
number, by changing the state of the viscous boundary layer. At flight Reynolds number, 
however, even small levels of roughness can result in considerable performance losses, as was 
demonstrated in reference 18. The Reynolds number range explored in this study was small and 
considerably lower than that experienced by full-scale wings at flight condition?, it is not known 
how much different flight test results would be. 
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4. CONCLUSIONS . 

Wind tunn el tests were conducted to investigate the effect of 20 ice shapes on the aerodynamic 
performance of a swept finite wing model. The wing consisted of an 8.7% thick airfoil section, 
which remained constant from root to tip. The wing was tapered and had a leading-edge (LE) 
sweep of 28°, a trailing-edge sweep of 15.6°, aspect ratio of 6.8, and -4° geometric twist at the 
tip (washout). Tests were performed with the clean wing, six ice shapes castings obtained from 
ice accretion experiments at the NASA Glenn Icing Research Tunnel (IRT), and seven smooth 
and seven rough LEWICE ice shapes. The LEWICE ice shapes were obtained for the same icing 
conditions as those used in the IRT ice accretion tests. One LEWICE ice shape was defined for 
each icing condition using airfoils sections . crmal to the wing LE at five spanwise stations. For 
icing condition 5, an additional LEWICE ice shape was defined using stream wise airfoil 
sections, as discussed in section 2.2.3 of this report. Roughness effects for the LEWICE ice 
shapes were simulated with 36-size grit. Tests were performed with the clean and iced wing in 
the low-speed 7- x 10-ft wind tunnel at Wichita State University. A complete set of force and 
moment coefficients were obtained along with aileron hinge moments and pressure distributions 
for a range of test conditions. Test conditions included Reynolds number of 1 .8 million based on 
the wing mean aerodynamic chord yielding test Reynolds numbers well below typical flight 
Reynolds numbers, angles of attack in the range of -8° to 20°, and aileron deflections of -15°, 
-10 3 , -5°, -2.5°, 0°, 2.5°, 5°, 10°, 15°, and 20°. From the results presented, the following 
conclusions are drawn. 

4, ! PERFORMANCE OF CLEAN WING . 

The clean wing C r ...h was 0.87 and occurred at a of 13.8°. The linear lift slope was 0.077 per 
degree or 4.41 per radian. Drag coefficient (Cd) varied from a minimum value of 0.0063 at a of 
I°to0.15ataofl3.8°. The maximum lift to drag ratio for the swept wing was 27.5 at a of 3 .2°. 
The hinge-moment coefficient (Ch) with the aileron in the neutral position varied from 0.079 at 
a of -8° to -0.02 at a of 0°, to -0.098 at a of 15.8°. TK Ch increased gradually through the 
nonlinear portion of the lift curve ud attained a maximum slope of -0.0195 per degree. The 
clean wing experienced a leading edge long bubble/vortcx stall, which progressed from the 
center of the wing to the tip and then to the root as a was increased. 

4 2 EFFECTS OF IRT ICF SHAPE C iSTINGS. 


The six IRT ice shape castings included 2-, 10- and 22.5-min glaze ice accretions with 
incomplete and complete scallop features and a 5 -min rime ice shape. The glaze ice snapes 
resulted in 11% to 93.6% reduction in G „.n ana 8% to nearly 56.5% reduction in a.uu with 
respect to the clean wing. These ice shapes increased ♦he clean wing CD,mm by 200% to 3533%. 
In many cases, the behavior of the aileron Ch for the clean wing was considerably altered by the 
glaze ice accretions. However, in all cases, Ch of the iced wing remained within the maximum 
and minimum limits defined by the clean wing aileron Ch. The 5-min rime ice shape improved 
the wing G ...i: by approximately 3% and the cuuu by 14.5%. The iced wing Cn,m.n was 133% 
greater than that of the clean wing. Aileron Ch for the 5-min rime case were similar to that of the 
clean wing. 
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For the IRT-CS22 ice sbtpe with the large scallops, the results presented showed that the gaps 
between the scallops increased ’ <1 ^gradation. This was due to flow leakage from the near 

stagnation region between the ice shape horns to the low-pressure region dc wnstream of the ice 
horns. This is a preliminary finding and requires further investigation before more general 
conclusions can be drawn regarding the effects of gaps in scalloped ice shapes. 

4.3 DEFINITION OF LEW ICE ICE SHAPES . 

A methodobgy for defining three-dimensional simulated ice shapes for a swept finite wing using 
the two-dimensional LE\ ICE ice accretion code was presented. 

4.4 EFFECTS OP SMOOTH i TWICE ICE SHAPES . 

The glaze ice shapes resulted in 14% to 44% reduction in The a«uii for the iced wing 

ranged from 13.4° to 14.8°. The iced wing Cu.mm was 183% to 2367% greater than that of the 
clean wing. In most cases, the aileron Ch for the wing with the glaze LEWICE ice shapes were 
more positive or more negative than that of the clean wing for a-range of -7° to +13°. However, 
in all cases, the Ch of the iced wing remained w ithin the maximum and minimum limits defined 
by the clean wing aileron Ch. The 5 -min rime LEWICE ice shape improved C r f 1l ll by 
approximately 8% and the clean wing out.n was increased by 1°. The Cd.^u for this ice shape 
was 67% greater than that of the clean wing. Aileron Ch for the o-min rime LEWICE ice shape 
were similar to that of the clean wing. 

4.5 EFFECTS OF ROUGH LEWICE ICF SHAPES . 

The rough glaze ice shapes resulted in 17% to 61% reduction in C i ...n . The a*uu for the iced 
wi”.g ranged from 12.6° to 13.7°. The iced wing Cu min was 233% to 3317% greater than that of 
the clean wing. In most cases, the aileron Ch for the wing with the rough glaze LEWICE ice 
shapes were more positive or more negative than those of clean wing for a-range of -7° to +13°. 
However, in all cases, Ch of the iced wing remained within the maximum and minimum limits 
defined by the clean wing aileron Ch- The 5-min rime LEWICE ice shape improved the clean 
wing Cl^uu by approximately 13% and increase by 2°. The Cu^un for this ice shape was 
1 17% greater than that of the clean wing. Aileron Ch for the 5-min rime LEWICE ice shape 
were in general similar to that of the clean wing. 

4.6 ROUGH VERSUS SMOOTH 1..FWIC1: ICE SH \FES. 


In gcr the smooth and rough LEWICE ice shapes produced similar lift curves for all but two 
of the cn ice shapes tested. The two cases where considerable differences in lift t.uavior 
were observed were the CS22N and CS22S cases. With the exception of the 5-min rime ice 
shape, the addition of roughness resulted in lower lift coefficients at stall. The change in Cl (i.e., 
rough LEWICE Ci^uii - smooth LEWICE was as follows: -0.1, -0.03, +0.D4, -0.03, -0.08, 
-0.2, and -0.08 for the CS10, lb 10, SC5, CS2, CS22N, CS22S, and 1PSF22 cases. With the 
exception of the CS22N and CS22S cases, the rough and smooth LEWICE ice shapes resulted in 
similar Cd and Ch. 
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The trends in aerodynamic performance losses observed with the smooth and rough LEW1CE i© 
shapes tested were consistent with results from other experimental studies involving LEWICi 
ice shapes. 

4, ICING RESEARCH VERSUS ROUGH L EV» H'lL ICE SH vPl.S . 

With the exception of icing condition 5 (22.5-min complete scallop case), the rough LEWICE i© 
shapes caused similar loss in lift compared with the 1RT ice J>ape castings. The dfference u 
iced wing (’• defined by (AC t = 1RT casting CLsau - rough LEWICE CL^u) was -0.07 
-C.01, -0.08, h- 6.05, -0.35, -0.28, and -0.08 for CS10, IS 10, SC5, CS2, CS22N, CS22S, am 
IPSF22 respectively. 

In most cases, the wing drag increase caused by the LEWICE and IRT ice shapes was in goo< 
correlation. The difference in iced wing Cn^. „ (i.e.. Rough LEWICE Cr» — - IRT ice shapi 
Cd^h) was 0.004, -0.008, -0.001, -K3.0C2, -0.036, -0.013, and -0.013 for CS10, IS10, SC5, CS2 
CS22N, CS22S, and 1PSF22 respectively. 

With the exception of icing condition 5, the trends in aileron Ch obtained with the rougl 
LEWICE and IRT ice shapes were in good agreement. Fcr icing conditions 2 and 3, the IRT am 
LEWICE Ch magnitudes were in good correlation throughout the a -range. For icing condition 
1, 4, and 6, the LEWICE ice shapes resulted in larger Ch over regions A and B of the Ch curv< 
compared with the IRT ice shapes. 

In nummary, the IRT and the rough LEWICE ice shapes produced similar aerodynamic effect 
for five out of the six icing conditions tested. This is an important finding considering th 
observed differences between the actual and the simulated ice shapes. It must be stressed 
however, that considerable more experimental work, is needed with a range of ice accretions ti 
determine the generality of this finding. 

*4 £ AILERON PERFORMANCE . 

For the clean win° and, for practically, all iced wing rases presented, the aileron remain© 
effective in increasing and decreasing lift with aileron deflection. The ice shapes cans© 
considerable changes to the Ch of the clean wing. However, in all cases, the iced wing C| 
remained within the maximum and minimum limits of the clean wing Ch. 

4/: RFYNQl PS NUMBER EFFECTS . 

In general, for the low Reynolds number (Re) range of 0.5 to 1.8 million used in this study, th 
effects of Re on iced wing lift performance were small. The only notable Re effect was in th 
behavior of Q ^ for the 5 -min rune ice shapes. For these ice shapes, lift near stall decreased a 
Re was increased from 1.0 to 1.8 million. 
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APPENDIX A— RUN LOG FOR ICING TEST AT NASA GLENN ICING 

RESEARCH TUNNEL 
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APPENDIX B— RUN LOG FOR AERODYNAMIC INVESTIGATION AT 
WICHITA STATE UNIVERSITY 


6a (deg) 


0 



Confinuration ,.e 


Clean 


Clean 


Clean 


Clean 


Clean 


Clean 


Clean 


Clean 


Clean 


Clean 


Clean 


IRT-CS10 (Irt-icel) 


IRT-CS10 (irt-tcel) 


IRT-CS10 (irt-icel) 


IRT-CS10 (irt-icel) 


IRT-CS10 (irt-icel) 


IRT-CS10 (Irt-icel) 


IRT-CS10 (Irt-icel) 


IRT-CS10 (irt-icel) 


IRT-CS10 (irt-icel) 


IRT-CS10 (irt-icel) 


IRT-CS10 (Irt-icel) 


Clean 


Clean 


Clean 


IRT-IS1U (irt-ice2) 


IRT-IS10 (irt-ice2) 


IRT-IS10 (irt ice2) 


IRT-IS10 (Irt-lce2) 


IRT-IS10 (irt-ice2) 


IRT-IS10 (irt-l ce2) 


IRT-IS10 (irt-lce2) 


Clean 


IRT-CS22 (irt-lce5) 


IRT-CS22 (Irt-iceb) 


IRT-CS22 (Irt-lce5) 


IRT-CS22 (Irt-ice5) 


IRT-CS22 (irt-ice5) 


IRT-CS22 (irt-ice5) 


IRT-CS22 (Irt-lce5) 



B-l 































































































































Date 


1 -May-02 


2-May-02 


3-May-02 


6-May-02 


Configuration Name 


IRT-CS22 (irt-ice5) 


IRT-CS22 (irt-ice5) 


IRT-CS22 (irt-ice5) 


IRT-CS22 (irt-iceS) 


Clean 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (Irt-lce3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (Irt-lce3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


Clean 


IRT-CS2 (irt tce4) 


IRT-CS2 (Irt-lce4) 


IRT-CS2 (irt-ice4} 


IRT-CS2 (Irt-lce4) 


IRT-CS2 (irt-ice4) 


IRT-CS2 (irt-ice4) 


IRT-CS2 (irt-ice4) 


IRT-CS2 (Irt-lce4) 


IRT-CS2 (irt-lce4) 


IRT-CS2 (irt-lce4) 


IRT-CS2 (Irt-lce4) 


Clean 


Clean 


IRT-IPSF22 (Irt-lce6) 


IRT-IPSF22 (irt-lce6< 


IRT-IPSF22 (Irt-lce6) 




IRT-IPSF22 (irt-ice6) 


IRT-IPSF22 (Irt-lceO) 


IRT-IPSF22 (Irt-ice6) 


IRT-IPSF22 (irt-lce6) 


IRT-IPSF22 (irt-ice6) 


IRT-IPSF22 iirt-lce6) 


Comments 


Baseline repeat 


Loose tape 


Loose tape 


Repeat R056 and R057 


Baseline repeat 


Baseline repeat 


Baseline repeat 


Repeat R08G 


Max out s 


I = 152°F 


B-2 










































































































































7-May-02 


8-May-02 


9- May-02 


10- May-02 


13-May-02 


14-May-02 


Run No. I 6 a (deg) Configuration Name 


15 IRT-IPSF22 (irt-ice6) 


IRT-IPSF22 (irt-ice6) 


IRT-IPSF22 (irt-iceO) 


Clean 


Clean 


Clean 


Clean 


IRT-IS10 (irt-ice2) 


IRT-IS10 (irt-ice2) 


IRT-IS10 (irt-ice2) 


IRT-IS10 (irt-lce2) 


IRT-IS10 (irt-ice2) 


IRT-IS10 (irt-ice2) 


IRT-IS10 (irt-lce2) 


IRT-IS10 (Irt-lce2) 


IRT-CS22 (Irt-lce5) 


IRT-CS22 (irt-lce5) 


IRT-CS22 (irt-ice5) 


IRT-CS22 (lrt-ice5) 


IRT-SC5 (irt-lce3) 


IRT-SC5 (Irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (irt-ice3) 


IRT-SC5 (Irt-ice3) 
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Clean 
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LS-IS10 (Is-ice2) 


LS-IS10 (Is-ice2) 


LS-IS10 (Is-ice2) 


LS-IS10 (Is-ice2) 


LS-IS10 (ls-ice2) 
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Loose tape. Terminate run 


Repeat R096 


Baseline repeat (Reunc = 1 -8x10®) 


ft 


IE1HES3E 


Repeal R032 


Repeat R033 


Repeat R034 


ReuAc = 0.5x10® 


Re*uc= 15x10® 


Reu*c = 1.0x10® 


R^uac = 1.8x10® 


Rei 4 Ac - 1 .8x10 (tape loose) 


Re UAC = 1.8x10®; Repeat R115 


Reutc - ' 5x10® 


RGuac = 1 > 


Baseline repeat 


Baseline repeat 


Baseline (using Model Constants Table 2). 


Baseline repeat 


Baseline repeat 


Clean 


LS-SC5 (Is-ice3) 


Baseline repeat 
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Date 

Run No. 

6 a (deg) 

Configuration Name 

20- May-02 
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LR-CS2 (Ir-ice4) 
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LR-CS2 (Ir-lce4) 
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LR-CS2 (Ir-lce4) 
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-2.5 

MES&mm 
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2.5 

LR-CS2 (Ir-lce4) 

21 -May-02 

232 

0 

LR-CS2 (Ir-lce4) 
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2.5 

LR-CS2 (Ir-lce4) 
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LR-CS2 (Ir-ice4) 
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1C 

LR-CS2 (Ir-ice4) 
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LR-CS2 (Ir-lce4) 
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LR-CS2 (Ir-ice4) 
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0 
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4 _ 4 _ i»*uy*02 
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0 

LR-1S10 (Ir-ice2) 
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-5 

LR-IS10 (Ir-ice2) 


242 

-2.5 

LR-IS10 (Ir-lce2) 
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2.5 

LR-IS10 (Ir-lce2) 
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5 

LR-IS10 (lr-ice2) 
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10 

LR-IS10 (Ir-ice2) 
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15 

LR-IS10 (Ir-lce2) 


247 

20 

■MUMIWW 

24-May-02 

264 

0 

Clean 

28-May-02 

276 

0 

Clean 

29-May-02 

291 

0 

Clean 

30- May-02 

299 

0 

Clean 

31 -May-02 

314 

0 

Clean 


315 

-15 

Clean (fixed transition) 


316 

-5 

Cle- -■ (fixed transition) 


317 1 

1 o 

Clean (fixed transition) 


318 1 

5 ; Clean (fixed transition) 


319 

* V 

Clean (fixed transition) 

3-Jun-02 

325 

9 ' ] 

Clean 


328 

O 

LS-CS10 (Is-icel) 


329 

-15 

LS-CS10 (Is-lcel) 


330 

-10 

LS-CS10 (Is-lcel) 


331 

-5 

LS-CS10 (Is-icel) 


332 

2 5 

LS-CS10 (Is-icel) 


333 

2.5 

LS-CS10 (Is-icel) 

4-Jun-02 
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0 

LS-CS10 (Is-icel) 
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5 

LS-CS10 (Is-lcel) 
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10 

LS-CS10 (Is-icel) 


337 

15 

LS-CS10 (Is-icel) 


338 ! 

20 ^ 

LS-CS10 (Is-icel) 


B-5 



Baseline repeat 

Use static tare R1113 
Repeated R1 99 



repeat 

B'sellne repeat 

Baseline repeat 

Baseline repeat 
Baseline repeat 

2 layers of aluminum strips at 2%c 
2 layers of aluminum strips at 2%c 
2 layers of aluminum strips at 2%c 
2 layers of akiminum strips at 2%c 
2 layers of aluminum strips at 2%c 



















































































































































Date 


4 June-02 


5-Jun-02 


6-Jun-02 


10-Jun-02 


11-Jun-02 


Run No. 

8a (deg) 

Configuration Name 

339 

0 

LS-IPSF22 (Is-ice6) 

340 

-15 

LS-IPSF22 (Is-ice6) 

341 

-10 

LS-IPSF22 (ls-ice6) 

342 

-5 

LS-IPSF22 (ls-ice6) 

343 

-2.5 

LS-IPSF22 (ls-lce6) 

344 

2.5 

LS-IPSF22 (ls-ice6) 

345 

5 

LS-IPSF22 (ls-ice6) 

346 

10 

LS-IPSF22 (ls-ice6) 
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15 

LS-IPSF22 (Is-ice6) 

348 

0 

LS-IPSF22 (ls-ice6) 

349 

15 

LS-IPSF22 (ls-ice6) 

350 

20 

LS-IPSF22 (ls-ice6) 

351 

10 

LS-IPSF22 (ls-ice6) 

352 

0 

LR-CS10 (Ir-icel) 

353 

-15 

LR-CS10 (Ir-icel) 

354 

-10 

LR-CS10 (Ir-icel) 

355 

-5 

LR-CS10 (Ir-icel) 

356 

-2.5 

LR-CS10 (Ir-icel) 

357 

2.5 

LR-CS10 (Ir-icel) 

358 

5 

LR-CS10 (Ir-icel) 

359 

10 

LR-CS10 (Ir-icel) 

360 

15 

LR-CS10 (Ir-icel) 

361 

20 

LR-CS10 (Ir-icel) 

362 

20 

LR-CS10 (Ir-icel) 

363 

0 

LR-CS10 (Ir-lcel) 

364 

0 

LR-IPSF22 (Ir-ice6) 

365 

20 

LR-IPSF22 (Ir-lce6) 

366 

15 

LR-IPSF22 (Ir-ice6) 

367 

10 

LR-IPSF22 (Ir-lce6) 

368 

5 

LR-IPSF22 (Ir-ice6) 

369 

2.5 

LR-IPSF22 (Ir-ice6) 

370 

-2 5 

LR-IPSF22 (Ir-6ce6) 

371 

-5 

LR-IPSF22 (Ir-ice6) 

372 

-10 

LR-IPSF22 (Ir-lce6) 

373 

-15 

LR-IPSF22 (Ir-lce6) 

374 

0 


376 

0 

Clean 

390 

0 

Clean 

396 

0 

Clean 

397 

0 

Clean 

399 

0 


400 

o ! 

| IRT-IS10 (irt-ice2) 


Comments 


Repeat R347 


Repeat R346 


Use Static Tare R1116 


Aileron binding. Terminate run 


Repeat R361 


Baseline repeat 


Baseline repeat 


Baseline repeat (automation wrong) 


Baseline repeat 


B-6 











































































































































Date 


1 2- Jun-02 


1 3- Jun-02 


17- Jun-02 


Run No. 

6a (deg) 

Configuration Name 

402 

0 

Clean 

403 

0 

LS-CS22N (Is-lce5n) 

404 

20 

LS-CS22N (Is-iceSn) 

405 

15 

LS-CS22N (Is-lce5n) 

406 

10 

L S-CS22N (ls-ice5n) 

407 

5 

LS-CS22N (ls-tce5n) 

406 

5 

LS-CS22N (ls-lce5n) 

409 

2.5 

LS-CS22N (ls-ice5n) 

410 

-2.5 

LS-CS22N (Is-lce5n) 

411 

-5 

LS-CS22N (ls-ice5n) 

412 

-10 

LS-CS22N (ls-ice5n) 

413 

-15 

LS-CS22N (Is-lce5n) 

414 

0 

LS-CS22N (ls-ice5n) 

417 

0 

LS-CS22N (Is-ice5n) 

418 

0 

LS-CS22N (ls-lce5n) 

419 

0 

LS-CS22N (Is-lce5n) 

420 

0 

LR-IPSF22 (lr-ice6) 

421 

0 

LR-IPSF22 (Ir-ice6) 

422 

0 


423 

0 

LR-IPSF22 (Ir-ice6) 

424 

0 

LR-IS10 (Ir-lce2) 

425 

0 

LR-SC5 (k-ice3) 

420 

0 

LR-SC5 (lr-lce3) 

427 

0 

LR-SC5 (Ir-ice3) 

428 

0 

LR-SC5 (Ir-ice3) 

429 

0 

LR-SC5 (Ir-lce3) 

430 

0 

LS-CS22S (ls-lce5s) 1 

431 

0 

LS-CS22S (Is-ice5s) 

432 

20 

LS-CS22S (ls-lce5s) 

433 

15 

LS-CS22S (Is Ice5s) 

434 

10 

LS-CS22S (Is-ice5s) 

435 ! 

5 

LS-CS22S (Is-lce5s) 

435 

2.5 

LS-CS22S (Is-lce5s) 

437 

-2.5 

LS-CS22S (Is-lce5s) 

4.38 

-5 

LS-CS22S (Is-lce5s) 

439 

-10 

LS-CS22S (Is-lce5s) 

440 

-10 

LS-CS22S (Is-ice5s) 

441 

-15 

L5> CS22S (Is-lce5s) 

442 

0 

LS-CS22S (Is-lce5s) 

444 

0 

Clean 

445 

0 

LR-CS22N (Ir-lce5n) 


Comments 


Baseline repeat 


Horn on pressure side split due to heat. 
Was glued back and reinforced with 
screws. 


Repeat of R407 


Repeat R364 and R374 


ReuAc= 1.0x10® 


Repeat R196 and R240 


Repeat R213 and R216 


Re-installed ice shape. Repeat R425 


Re^c = 1 0x10 s 


R&mac — 0.5x1 0 6 


Aileron binding. Terminate run 


Repeat R439 


Baseline repeat 
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Date 


1?-Jun-02 


18-Jun-02 


19-Jun-02 


20-Jun-02 


21 - Jun-02 


25-Jun-02 


Run No. 


446 


447 


6 a (deg) 


20 


Configuration Name 


LR-CS22N (Ir-ice5n) 


LR-CS22N (Ir-lce5n) 


LR-CS22N (Ir-ice5n) 


LR-CS22N (Ir-lce5n) 


Comments 


LR-CS22N 


LP-CS22N 


LR-CS22N 


LR-CS22N 


LR-CS22N 


LR-CS22N 


LR-Cw.22S 


LR-CS22S 


LR-CS22S 


LR-CS22S 


LR-CS22S 


LR-CS22S 


LR-CS22S 


LR-CS22S 


LR-CS22S 


LR-CS22S 


LR-CS22S 


(Ir-lce5n) 


(Ir-lce5n) 


(Ir-tce5n) 


(Ir-ice5n) 


(Ir-ice5n) 


(Ir-lce5n) 


(Ir-lce5s) 


(lr-lce5s) 


(Ir-ice5s) 


(Ir-ice5s) 


(Ir-ice5s) 


(Ir-ice5s) 


(Ir-ice5s) 


(Ir-lce5s) 


(Ir-ice5s) 


(Ir-ice5s) 


(Ir-lce5s) 


Loose tape. Terminate run 


Repeat R451 


LR-CS22S (lr-ice5s) 


LR-CS22S (Ir-ice5s) 


LR-CS22S (Ir-lce5s) 


LR-CS22S (Ir-lce5s) 


LR-CS22S (Ir-ice5s) 


LR-CS22S (Ir-ice5s) 


LR-CS22S (Ir-lce5s) 


LR-CS22S (Ir-ice5s) 


Clean 


LR-CS22N (Ir-ice5n) 


LR-CS22N (Ir-lce5n) 


LR-CS22N (Ir-ice5n) 


LR-CS22N (Ir-lce5n) 


LR-CS10 (Ir-tcel ) 


LR-CS10 (Ir-icel) 


LR-CS10 (Ir-lcel) 


LR-CS10 (Ir-icel) 


Clean 


Clean 


IRT-CS10 (Irt-icel) 


Loose tape. Terminate run 


Repeat R460, loose tape 


Repeat R460, loose tape 


Repeat R460, loose tape 


Repeat R460 


Repeat R464 


Sandpaper started to peel off at about 1 7°, 
but data looks good. 


Lost sandpaper 


Repeat R474 


iehhese 


ReuAc = 0.5x10 


Repeat R010 


Repeat R445 


Rbmac — 1 .0x10 


Re^c = 15x10' 


Repeat R352 


Re^c = 0.5x1 0 6 


Baseline repeat 


Flow-vlz, Re mc = 1.5x10® 


Flow-viz. Re^c = 1.5x10® 
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Date 


25-Jun-02 


Run No. 


8a (deg) 


0 



Configuration Name 


IRT-IPSF22 (lrt-lce6) 


IRT-IS10 (irt-lce2) 


IRT-SC5 (irt-ice3) 


IRT-CS2 (Irt-lce4) 


LR-IS10 (Ir-lce2) 


LR-SC5 (Ir-ice3) 


LR-CS? (Ir-ice4) 


LR-CS10 (Ir-icel) 


LR-IPSF22 (Ir-ice6) 


LR-CS22S (Ir-lceS?) 


LR CS22N (lr-ice5n) 


IRT-CS22 (Irt-lce5) 


Clear 


Clean 


Clean 


Comments 


Flow-viz, Reutc = 1.5x10 


Flow-viz, ReuAc = 1.5x10 


Flow-viz, Re^c = 1.5x10 


Flow-viz, Re^Ac = 1.5x10® 



LR-IPSF22 (Ir-ice6) 


LR-IPSF22 (Ir-ice6) 


IRT-IPSF22 (irt Ice6) 


IRT-IPSF22 (irt-lce6) 


IRT-IPSF22 (Irt-lce6) 


IRT-IPSF22 (irt-ice6) 


LS-IPSF22 (ls-lce6) 


LS-IPSF22 (Is-lce6) 


LS-IPSF22 (Is-iceS) 


LS-IPSF22 (Is-ice6) 


Clean 


LS-CS22S (Is-ice5s) 


LS-CS22S (Is-lce5s) 


LS-IPSF22 (Is-lce6) 


IRT-CS22 (irt-ice5) 


Clean 


Clean 


Clean 


Clean 


Clean 


IRT-CS22 (irt-lce5) 


Flow-viz, ReuAc= 1.5x10® 


Flow-viz, Re ht A C = 1.5x10® 




Flow-viz, Re^c = 1.5x10 


Repeat R520, Re^Ac = 1.0x10 


Repeat R521, fleuwc = 1.5x10® 


Flow-viz, ReuAc= 1.5x10® 


Flow-viz, ReuAc= 1.5x10® 


Take data with yam tufts on model, 
Re^Ac - 1.5x10® 


Basel!, te repeat 


Added 24-grlt to tip of horns of Ice shape 


Added 24-grit to tip of horns of ice shape 


Added 24-grtt to tip of horns of ice shape 


Rcuac = 1 .8x10® 


flew, v:« 1.5x10® 



R^uac = 


fle*Mc= 18x10® 


R&mac — 15x10® 


fle*4Ac= 1.0x10® 


fle«A C = 0.5x10® 


Baseline repeat 


Flow-viz, ReuAc = 1.5x10®, Repeat R582 


Flow-viz, ReuAc= 1-5x10® 


Flow-viz, fle*Mc= 1.5x10 


Checking flow angularity probe installation, 
q = 15 psf 


7 -hole flow angularity probes 


7 -hole flow angularity probes 


7-hoie flow angularity probes 


Gap opened 


Gap opened 


U-9/B-10 












































































































































APPENDIX C— PRESSURE DATA FROM ICING TESTS 


Surface pressures at 15%, 50%, and 85% semispan locations were obtained prior to the ice 
accretion tests. At velocity of approximately 150 mph, with the wing in normal position (refer to 
section 2.1.3), pressures were measured at angles of attack of -3°, -2.5°, -2°, -1.5°, -1°, -0.5°, 0°, 
0.5°, 1°, 1.5°, 2°, 4°, 6°. 8°, 10°, 11°, 12°, 13°, 14°, and 14.5°. Pressure coefficient <Cp) 
distributions of selected a are presented in figure C-l. 

Pressure measurements were performed with the icing research tunnel (IRT) electronically 
scanned pressure (ESP) system. Six 32-port (±5 psid) ESP modules were available in the IRT, 
providing a total of 192 pressure channels. One port in each module was used for a check 
pressure; thus 31 chant 'Is per module were available for test data, or a total of 186 ports. The 
E r P system applied a three-point pressure calibration to all port transducers. This on-line three- 
point calibration ensured that measurement errors were not greater than ±0.1% of full-scale. The 
standard calioration interval was every 400 cycles (approximately 15 minutes). 

The experimental presv i re data were used to validate Cp distributions from LEWICE analyses 
and wind btnnel tests at Wichita State University. 
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APPENDIX D— FLOW ANGULARITY STUDIES 


Two, seven-hole flow angularity probes were installed on the clean wing at two statio. 
corresponding to approximately 16% and 80% semispan. Seven-hole probes can measure t 
three components of velocity, the total pressure, and the static pressure at a point in the flo 
The flow probe provides results with high accuracy for angles as high as 75°. The data from ti 
seven-hole probes were used to estimate the local angle of attack of the wing due to geomeh 
twist of the wing, presence of the streamlined body, and vortices at the wing tip. Angie of ana 
sweeps were conducted for three aileron deflections (-15°, 0°, and 20°) and at Reynolds numb 
of 1.8 million based on wiu fe mean aerodynamic chord. Pressures measured from the seven-hc 
probes were recorded, installation of the seven-hole flow probes is shown in figures 2-52 
2-55. 


In order for the probes to provide accurate downwash angle measurements, they were attach 
onto the wing leading edge and aligned to the free -stream flow direction. Note that this line 
direction is not concurrent to the wing chcrdwise direction, and therefore, a set of brackets w 
specifically made to accommodate the installation. 

Effects of flow probes on the wing lift and drag coefficients are demonstrated in figures D-l aj 
D-2. It is observed that the flow probes improved lift performance of the wing; i.e., 1 
coefficient (Cl) was increased throughout the a-range and as a result Cl^uii was increased 
well. This is because the flow probes generated vortices that re-energized the boundary lay* 
which delayed transition and flow separation; hence, improved Cumii. From figure D-l, it 
observed that the offsets between the Cl values obtained with and without the flew protx 
amongs t all three aileron deflections of -15°, 0°, and 20 3 , were very similar, and that holier 1 
curves were produced with higher 5a- It is of interest to draw attention to the fact that the 1 
curve generated whh the flow probes at 5 a of -15° behaved very much like that from without t 
probes at 5 a of 0°. 

The same vortices that re-energized the b~undary layer also inadvertently increased dr 
coefficient (Cd), as shown in f ‘jure D-2. In general, the increment in drag contributed by t 
probes was only accountable when the wing was experiencing a positive a. Consider the case 
drag coefficients generated by the wing with the flow probes installed at 5 a of 20°, its Cd vaiu 
were the highest amongst all three aileron deflection configurations, with and without the profcx 
yet oniy at the range of positive a from 2° onwards. In contrary, its Cd values from a of -4° 
-8° were the lowest amongst all configurations tested. 

Downwash (c. up wash) angle as a function of geometric angle of attack is depicted 
figures D-3 and D-4. Results from figure D-3 were collected from the probe at 16% semisp 
location, whereas .hose from figure D-4 were collected from the one at 80% semispan. In be 
cases, the most negative aileron deflection (5a = -15°) configuration generated the larg< 
amo unt of downwash angles for the positive a-sweep. Under the same argument, the t< 
configuration of the most positive aileron deflection (5a - 20°) generated the largest upwa 
angles (negative downwash angles) during the negative a-sweep. Note that the downwash an{ 
curve in figure D-4 obtained near the tip of the semispan behaved less linearly than the o 
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observed towards the wing root, as shown in figure D-3. Also, the slope of the downwash angle 
with respect to the angle of attack observed at the near-tip section is smaller than that of the near- 
root section. 




FIGURE D-l . EFFECTS OF FLOW PROBES FIGURE D-2. EFFECTS OF FLOW PROBES 
ON CL; Re= \.8xl0 6 ON CD; Re= 1.8xl0 6 



-10 -5 0 5 10 la 

Angie of attack (dag) 



FIGURE D-3. RELATIONSHIP OF 
DOWNWASH ANGLE WITH GEOMETRIC 
ANGLE OF ATTACK; 16% SEMISPAN; 
Re= 1.8xl0 6 


FIGURE D-4. RELATIONSHIP OF 
DOWNWASH ANGLE WITH GEOMETRIC 
ANGLE OF ATTACK; 80% SEMISPAN; 
J2e=1.8xl0 6 
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APPENDIX E— SIMULATED FROST USING SANDPAPER 


As per Federal Aviation Administration request, a total of eight configurations of frost 
simulation using 2-grit sandpaper (times four coverage) were identified and investigated to study 
the effects of frost on aerodynamic performance of a swept wing. The use of sandpaper in 
simulating ice shape roughness aligns with current procedures used by airffamers during aircraft 
icing certification . 

Sandpaper tested were 40 grit and 150 grit Witn mean aerodynamic chord (MAC) of 1.25 ft, 
40-gnt provided normalized roughness (k/MAC) of 1x10 \ where the 150-grit sandpaper 
provided k/MAC of 2x1 O' 4 . Figure E-l shows the four different frost simulation coverage using 
40-grit sandpaper. The first configuration (figure E-l (a)) simulated frost coverage on the whole 
wing. From figures E-l(b) and E-l(c), the second and third configurations were to simulate frost 
coverage on the aft 87.5% and 35% of the wing, respectively. In addition, the fourth coverage 
(figure E-l(d)) was to simulate failed deicing fluid condition. All sandpaper coverage were 
applied on the wing upper surface only. 

Figure E-2 demonstrates the effects of simulated frost on Cl, Cd, Cm, and Ch of the wing using 
40-grit sandpaper, whereas the simulated frost effect using 150-grit roughness was shown in 
figure E-3. 


* Papadakis, M., Yeong, H.W., Chadrasekharan, V.R., Hinsloo, M., and Ratvasky, T.P., “Effects of Roughness on 
the Aerodynamic Performance of a Business Jet Tail,'' A1AA Paper 2002-0242, January 2002. 
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Pitching Moment Coefficient Coefficient 






Pitchinn Momorrt Coefficient ^ Coeffident 




(a) Lift coefficient (b) Drag coeffident 




(c) Pitching-moment coefficient (d) Hinge-moment coefficient 

FIGURE E-3. EFFECT OF FROST ON AERODYNAMIC PERFORMANCE; 150-GRIT 

SANDPAPER; Re = 1 .8xl0 6 ; 6 A = 0° 
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Lift coefficient 


APPENDIX F— EFFECT OF AILERON DEFLECTION (SUPPLEMENTARY) 



-10 -5 0 5 10 15 20 

Angie of attack (deg) 

(a) Lift coefficient (negative 6a) 


25 



(c) Hmge-moment coefficient (negative 6a) 




(b) Lift coefficient (positive 6a) 



-10 -5 0 5 10 I r 20 25 

Angie of attack (defc, 


(d) Hinge-moment coefficient (positive 6a) 



FIGURE F-l. EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; IRT-CSIO (ICING CONDITION 1); Re= 1.8xl0 6 
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Angi* of stuck (deg) 

(b) Lift coefficient (positive 5 a) 




(c) Hinge-moment coefficient (negative 6a) (d) Hinge-moment coefficient (positive 6a) 




FIGURE F-2. EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; IRT-IS10 (ICING CONDITION 2); Re= ’.SxlO 6 
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Lift coefficient Hinge moment c oefflclent Uft coeffident 


1.00 


Angie of stuck (deg) 

(a) Lift coefficient (negative 6 a) 


0*80 

0.60 

t 0 40 

| 0^0 

0 0.00 
£ 

3 -0^0 
-0.40 
-0.60 
-OJ80 


o(*j**j 

ItOkUM 


Angie of stuck (deg) 

(b) Lift coefficient (positive 6 a) 



(c) Hinge 

i.ao 
1.60 
1.40 
1.20 
1.00 4 


Angie of stuck (deg) 

moment coefficient (negative 6a) 
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(e) Lift coefficient vs 6 a 
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Aileron Deflection (deg) 

(f) Hinge-moment coefficient vs 6 a 


FIGURE F-3. EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; IRT-SC5 (ICING CONDITION 3); Re= 1.8xl0 6 
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-15 -10 -5 0 5 10 15 20 

Aileron Deflection (deg) 



(e) Lift coefficient vs 6 a 
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Aileron Reflection (dag) 

(0 Hinge-moment coefticient vs 8 a 


FIGURE F-4. EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; IRT-CS2 (ICING CONDITION 4); J& = 1.8xl0 6 
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FIGURE F-8. EFFECT OF AILERON DEFLECTION ON LIFT AND HL IGE-MOME1 
COEFFICIENTS; LS-CS2 (ICING CONDITION 4 ); Re= 1.8x10 s 
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FIGURE F-9. EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; LS-CS22N (ICING CONDITION 5);Re= 1.8xlU 6 
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MOURE F-10. EFFECT OF AILERON DEFLECTION ON LIF'" AND ?ilNGE-MOMLh 
COEFFICIENTS; LS-CS22S (ICING CONDITION 5); Re= 1.8x10° 
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Angta of attack (dag) 

(a) Lift coefficient (negative 5*) 


Angfcof Attack (dag) 

(b) Lift coefficient (positive Sa) 
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Ailaron Daflaction (dag) 
(e) Lift coefficient vs 5 a 
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(f) Hinge-moment coefficient vs 5 a 


FIGURE F-l 1 . EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; LS-IPSF22 (ICING CONDITION 6); Ke= 1.8xl0 6 
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Aiitron Ml«ction (dag) 

(e) Lift coefficient vs 5* 



AUorort Daftoctton (dag) 

(fj Hinge-moment coefficient vs S* 


FIGURE F-12. EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; LR-CS10 (ICING CONDITION I); Re= I.8xl0 6 
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Aileron Deflection (deg) 
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Aileron Deflection (dag) 


(e) Lift coefficient vs 6a 


1 »j Hinge-moment coefficient vs 6 a 


FIGURE F-13 EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; LR-IS10 (ICING CONDITION 2); Re= 1.8xl0 6 
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U , T OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
"• ~ - '-TS; LR-SC5 (ICING CONDITION 3); Re = 1.8xl0 6 
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(e) Lift coefficient vs 6* 


irj litnge-momem ccemctent vs o* 


FIGURE F-15. EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; LR-CS2 (ICING CONDITION 4); Re= 1.8xl0 6 
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FIGURE F-16. EFFECT OF AILERON DEFLECTION ON LIFT AND HINGE-MOMENT 
COEFFICIENTS; LR-CS22N (ICING CONDITION 5);Re= 1.8xl0 6 
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APPENDIX G— COMPUTATION OF LEWICE ICE SHAPE PROFILES 


The ice accretion code used in this study is LEWICE v2.0 [G-l], developed by NASA Glenn 
Research Center. LEWICE is commonly used in the industry to determine ice shape profiles, 
water droplet impingement patterns, water or icc mass flux, and chordwise extent of ice growth 
on the body of interest. Due to the physical modeling and computational methods employed in 
the software, LEWICE executes very rapidly and thus can be used to perform extensive analysis 
as long as the assumptions inherent in the software are not violated. This software uses the clean 
body geometry along with the aerodynamic and meteorological conditions to compute an ice 
shape. The output of the software is a 2-dimensional (2D) ice shape profile that car. be 
calculated for several locations along the span of the wing. The approach taken in this work was 
to calculate ice shapes at five locations along the span of the 3-dimensional (3D) wing and then 
construct a full 3D ice shape by lofting between sections. The five locations selected for the ice 
shape compulations were 0% (wing root), 15%, 50%, 85%, and 100% (wing tip) semispan. 

In order to use the 2D calculation for construction of the ice shape, adjustments were made to 
account for the 3D flow effects and how these effects impact the droplet trajectory, heat transfer, 
and ice growth calculations. From t' s work of Dorsch and Brun [G-2], the calculation of droplet 
impingement on a swept wing can be performed by taking a cross-section normal to the leading 
edge and adjusting the velocity by use of the following expression: 

V m = V^cosA (G-l) 

where V a is the velocity in the plane normal to wing leading edge. 

In addition to this velocity adjustment, a was modified to match the pressure profile at the 
leading-edge region as determined from a 3D Navier-Stokes calculation of pressures over the 
entire wing surface. 

The 3D Navier-Stokes computations were performed at Wichita State University with the clean 
wing. In he computations, the NASA Glenn Icing Research Tunnel (IRT) walls were included 
to simulate the tunnel wall effects on the wing flow field. The a used in the analysis were 4° and 
6°, to match the geometric angles of attack in the icing tests. Analysis pressures for stream wise 
sections at 15%, 50%, and 85% semispan were compared with experimental pressures 
distributions obtained in the IRT facility. Good correlation between experiment and analysis was 
demonstrated. From the computed flow fields, pressure distributions were obtained for the four 
wing sections normal to the wing leading edge (at 15%, 50%, and 85% semispan, and wing tip) 
and for the streamwise section at the wing root. 

The two adjustments described above allow the velocity and pressure profile o f the 2D 
calculation to be similar to those of the actual 3D flow. An alternative method to matching the 
pressure profile is to change the angle of attack in the 2D calculations so that the lift of the airfoil 
is equal to the lift of the wing at the cross section of interest or to simply adjust the angle of 
attack gee etrically. These approaches were not attempted for this investigation; rather, the 
former mch ad was employed throughout the study. 


G-l 



The exception to this approach was the calculation of the ice shape at the root of the wh'g. la 
that location, a section cut perpendicular to the leading edge would intersect the floor of the 
tunnel. In addition, the airflow in that region was affected by the presence of the tunnel fioor. 
Flow visualization suggests that the flow was parallel to the tunnel walls. As such, Me LIWICE 
calculation for the root section was performed using the geometry corresponding to a cross- 
sectional cut parallel to the tunnel walls. 

The temperature input for the LEWICE software was also adjusted to approximate the conditions 
present on the actual wing model. LEWICE requires the static temperature of the free stream as 
input. Since the input velocity was reduced, as indicated in equation G-l, the resulting total 
temperature profile for the airfoil would be reduced by an amount approximately equal to 



20 


(G-2) 


where C p is the specific heat of air. 

Thus, the input free-stream static temperature was modified such that the free-t-Tcam total 
temperature matched that of the 3D flow. The expression used to make this adjustment is 
provided in equation G-3. 


T SJ1 = T sJ + 


Vl 


20 


-sin 2 A 


(G-3) 


where T SrB and 2V are the static temperature in the plane normal to wing leading edge and free- 
stream plane, respectively. 

In summary, the aerodynamic and icing conditions input for LEWlCE analyses to generate the 
seven LEWICE ice shapes are documented in table G-l. 

REFERENCES 


G-l. Wright, W.B., < °r Manual for the NASA Glen Ice Accretion Code LEWICE (Version 
2.0),” NASA CR-1999-2094009, September 1999. 

G-2. Dorsch, R.G. and Brun, R.J., “A Method for Determining Cloud-Droplet Impingement on 
Swept Wings,” NACA TN 2931, April 1953. 
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TABLE G-l. FLOW AND ICING CONDITIONS INPUT FOR iJiWICE ICE SHAPES 





r 

(mm) 

RH 

(%) 

P italic 

(N/m 2 ) 

Iridic 

(K) 

Spanwi&e 

Section 

C 

(m) 

V 

(n>/sl 

a 

(deg) 

LS-CS10 

0.68 

20.0 

10.0 

100 

100000 

263 

Root 

0.640 

111.76 

4.0 

15% 

0.572 

98.678 

2.8 

50% 

0.440 

98.678 

22 

85% 

0.308 

98.678 

0.6 

Tip 

0.257 

98.678 

0.6 

LS-IS10 

0.65 

20.0 

10.0 

100 

100000 

263 

Root 

0.640 

67.055 

4.0 

15% 

0.572 

59.207 

2.8 

50% 

Eg® 

59.207 

22 

5% 


59.207 

0.6 

Tip 

0.257 

59.207 

0.6 

LS-SC5 

0.51 

14.5 

5.0 

100 

100000 

263 

Root 

0.640 

89.989 

4.0 

15% 

0.572 

79.456 

2.8 

50% 

0.440 

79.456 

12 

85% 

0.308 

79.456 

0.6 

Tip 

0257 

79.456 

0.6 

LS-CS2 

0.68 

20.0 

2.0 

100 

100000 

263 



111.76 

4.0 

15% 

0.572 

98.678 

2.8 

50% 

0.440 1 

98.678 

22 

85% 

0.308 

98.678 

0.6 

Tip 

C.257 

98.678 

0.6 

LS-CS22N 

0.68 

i 

20.0 

22.5 

100 

100000 

263 

Rcu 

0.640 

111.76 

wsm 

15% 

(■572 1 

98.6 7 8 

2.8 

50% ! A .440 

98.678 

22 

8-% 


98.678 

0.6 

Tip 

0.257 

98.678 

0.6 

LS-CS22S 

0.68 


22.5 

100 

100000 

263 

Root 

0.640 

111.76 

4.0 

15% 

0.580 

111.76 1 

2.8 

50% 

0.448 

111.76 

22 

85% 

0.314 
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APPENDIX H— COMPANION OF TWO- AND THREE-DIMENSIONAL WING ICE 

SHAPE TRACES AND LIFT DATA 

H.l COMPARISON OF TWO- AND THREE-DIMENSIONAL ICE SHAPES . 

The icing conditions for ice shape T RT-SC5, listed in tables 2-1 and 2-2 of this report were 
obtained by scaling the condidons used in an ice accretion test conducted with a two-aimerenanai 
(2D) 36-inch chord GLC-305 airfoil [H-l]. 

Tue scale condidons were determined udng the Ruff meth >d [H-2] for sea level testing with 
constant velocity. This method was develc*' id for scaling 2D models to permit simulating a 
larger model by testing a smallt model, o; tc permit tesing at one set of test conditions to 
simulate conditions that were not attainable in the facility being used. For either simulation, the 
objective is to produce a scaled ice shape that simulates the reference shape; in fact, when 
normalized by the model chord, the two shapes should match in overall size and in locutions and 
size of individual features. The Rut! method maicnes scale and reference values of tho*e 
nondhnensional similarity parameters that have been demonstrated to have the strongest 
influence on the ice accretion quantity and shape. The dimensionless simib*"‘\ parameters 
include the modified inertia parame f - r , Ko, which relates to the droplet trajectories; the 
accumulation parameter. A., which is a measure of the quantity of ice that can potentially accrete 
the model; and the freezing fraction, no, which is the proportion of water impinging at 
stagnation that actually freezes. Another energy parameter frat can be used is the water energy 
transfer parameter, <p. This j •ameter has units of temperature and collects all the *erms in the 
ene.gy balance that relate u> energy carried to the surface by water drops. For rime ice, it is 
convenient and acceptable to set the scale velocity equal to the reference, although oilier methods 
of finding scale velocity have been shown to give better results for glaze ice. The five equations 
formed by matching the scale and reference values of the three noudimensiooal similarity 
parameters, plus the water energy parameter and the velocity can be solved to fitu the scale test 

r-nndiiirtnK 


The method outlined in the previous paragraph has only been validated to scale 2D (straight 
wing) models. There have been no studies of 2D to three-dimensional (3D) (*wept wing) 
scaling, and a ctucial question is how valid are the 2D parameters, some of which appiy only to 
the stagnation line, might be for swept wings. However scaling betwe" 1 two geometrically 
similar swept wing models of different size was studied by ONERa [H-3 and H-4] by applying 
most of the same similarity parameters that Ruff later used. Aitnocgh direct comparisons of 
scale and reference ice shapes were not possible, photographs suggested that quantities and 
features of the reference accretions were simulated well by the scale tests. These results 
provided some er ouragoineut tiiat the Ruff method might be valid for 2D to "D scalirg. Tc test 
thus possibility in a very limited way, the present study included rime tests with a 2D model as a 
reference and compared those icc shapes with ones from a 3D model at scaled conditions. 

The strategy was to anply the Ruff scaling calculation* outlined above using rime tests w ith a 30- 
in chord, 2D GLC-305 model as the reference. The scale model wa* to have the same rarfou 
form as the reference but w ith a chord of 17.4 inches. Rime ice, even for swept wings, has a 
fairly simple shape , lacking the 3D scallops of swept wing glaze. Therefore, the assumption was 
that rims ice accretions would be no different whether on a 1 7.4-in chord 1 D airfoil or on a 3D 
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airfoil with a mean aerodynamic chord (MAC) of l" 1 .- itches, at least at the MAC. The lea 
edge radius for the GLC-305 airfoils was taken to be 0.0134c, where c is the chord, 
reference conditions and resulting scale conditions -re sho wn in the following table: 



c 

(in) 

V 

(mph, 

(deg) 

Ttoul 

(°F) 

MVD 

(pm) 

LWC 

(gm 3 ) 

Time 

(minutes) 

Po 

(%) 

Pc 

Reference 

36.0 

201.3 

6 

11.7 

:o.o 

0.40 j 

16.'!' 

66.3 

l.< 

Scaled 

17.4 

201.3 

6 

11.7 

12.8 

0.58 1 
1 

5.6 

i 

66.3 

l.< 


The relative thickn ess of ice accreted at the leading edge is proportional to the product of 1 
parameters: the collection efficiency at stagnation, po, the accumulation parameter. A*, ant 
stagnation freezing fraction, no. Thus, for time ice (no of 1), the product |3oA« gives a measu 
the relative size of the ice accretions. Note that for proper scaling, this product must nr 
between scale and reference. Other suniimty parameters are not important for rime scaling. 


Some adjustments to the calculated scale conditions were necessary. First, because the ! 
drop size determined by the Ruff method w»s smaller than the minimum value for whicl 
icing research tunnel (1RT) cloud was calibrated (14.5 pm), the minimum median volum 
diameter (MVD) for the IRT had to be substituted. For a straight wing, this increase in drop 
would result in a stagnation collection efficiency about 6% higher than the reference. Fa 
swept wing, the assumption was made that the change would have the same proportional e 
as that for the straight wing. Thus, the second adjustment w to reduce the spray time free 
to 5 minutes to compensate for the increased collection efficiency. Then the revised i 
conditions were: 



c 

(in) 

V 

(mph) 

(deg) 

Ttoul 

(°F) 

MVD 

(pm) 

LWC 

(g/m 3 ) 

Time 

(minutes) 

Po 

(%) 

Po- 

Reference 

3o.O 

201.3 

Si 

11.7 

20.0 

0.40 

16.7 

66.3 


Scaled 

17.4 

201.3 

S3 

11.7 

14.5 

0.58 

! 

70.2 



Although the scaling calculations were made using a chord of 17.4 inches, the icing tests 
actually performed with a swept GLC-305 airfoil with a MAC of 18.72 inches. Botl 
collection efficiency and accumulation D&rameter are dependent on chord. The local cho 
which icing traces were made also varied with spanwise location; therefore, the scale 
conditions shewn above were valid only for one location. In addition to the change in MAC 
liquid water content (LWC) actually tested was not the value desired due to an error in the 
calibration discovered after the completion of tests. Thus, while die spray bar conditions 
set to give an LWC of 0.58 gin 3 , as recommended from the scale calculations, the actual vai 
believed to have been 0.51 g/m\ 

Figure H-l compares a trace of the ice shape obtained from the icing tests with the 2D GLC 
with three traces of the ice shape obtained from the 3D sw«pt wing icing tests at the N. 
Gienn IRT facility. The traces of the 3D ice shape were obtained at wing stations A, B, and < 
shown in figure 2-5 of this report. The ice traces at stations A and B were taken normal t< 
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swept wing leaping edge, while the trace at station C was taken in the stream wise direction. A 
discussed in the main body of this report, the swept wing had a GLC-305 airfoil section that wa 
in the stream wise direction. Thus, sections taken normal to the leading edge of the swept win 
were slightly thicker than the streamwise GLC-305 section anu n^d a chord length that wa 
approximately 99% of the streamwise airfoil chord length. The airfoil and ice traces presented i 
figure H-l have been normalized with the normal or streamwise chord of each airfoil sectioi 
These nondimcnsional tracings show good agreement between the 2D trace and that of cut A c 
the swept wing, but the relative quantity of ice for cuts B and C are far less than that for the 21 
trace. The fot’owing table shows the relevant similarity parameters for the conditions actual! 
tested at the * ! *vC ana at each of the tracing locations: 



x/c 

FIGURE H-l. COMPARISON OF THE ICE SHAPES FOR THE 2D AND 3D MODELS 



c 

(in) 

V 

(mph) 

a 

(deg) 

Ttoui 

(°F) 

MVD 

(pm) 

LWC 

(&'m 3 ) 

Time 

(minutes) 

ESqI 

PoAc 

Reference 

36.00 

201.3 

6 

11.7 

20.0 

0.40 

16.7 

66.3 

1.065 

Scaled (MAC) 

18.72 

201.3 

6 

11.7 

14.5 

0.51 

5.0 

68.9 

0.8 1: 

Scaled (Cut A) 

13.80 


6 

11.7 

14.5 

0.5 1 

5.0 

74.1 

1.1K 

Scaled (Cut B) 

19.30 

201.3 

6 

11.7 

14.5 

0.51 

5.0 

68.3 

0.781 

Scaled (Cui C) 

E&3 

201.3 

6 

11.7 

14.5 

0.51 

5.0 

63.2 

0.55i 


Note that the scale value of fioAc is only to the reference at cut A, where the scale ar 
reference are within 1 1%. The parameter Ac is diiectly proportional to LWC, and the uncertain) 
in the LWC calibration for the IRT is generally quoted as ±10%, so these numbers aj 
reasonably close. The good agreement of the ice tracings for the 2D and 3D at cut A i 
therefore, consistent with the PoAc value there, and this result shows that for this rime case, . 
least, it is possible to scale effectively from 2D to 3D. The much smaller scale shapes at cuts 
and C are predicted by the significantly lower (3oAc values at those locations. 
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Although the results for nearly matched PoAc were encouraging, the tests were very limited in 
scope, including only rime ice at one set of reference conditions. For glaze ice, the scallop 
formation on swept wings of large sweep angles makes a direct comparison with 2D ice shapes 
impossible. However, additional tests for glaze ice at several freezing fractions are encouraged 
to determine if useful similarities in 2D and 3D shapes can be identified. For future tests, it 
should also be recognized that the matciting of reference and scale similarity parameters applies 
only at one swept wing spanwise location. 

H.2 COMPARISON OF 2D AND 3D CLEAN AMD ICED WING AERODYNAMIC 
PERFORMANCE. 


In this section, the aerodynamic lift performance of the clean and iced GLC-305 2D wing model 
presented in reference H-l is compared with the lift performance of the clean and iced GLC-305 
3D wing presented in the main body of this report. The 2D and 3D iced wing configurations 
considered below are the ones presented in figure H-l. The aerodynamic performance 
comparisons will be limited to linear lift slope and maximum lift coefficient. 

In comparing the aerodynamic performance of the two wing models, the following differences in 
model geometry' and flow-field behavior should be considered: 

• Model size and geometric features. The 2D model was a straight infinite span planar 
wing with a chord of 36 inches. The 3D model was a finite span swept wing with 
geometric twist and taper, and with a MAC of 18.72 inches. 

• The lift data for the 2D model were obtained at Reynolds numbers in the range of 3 to 7.5 
milli on compared to the 1.8 million Reynolds number used in the 3D wing tests. 

• The flow field about a finite swept vting is three-dimensional and is, therefore, inherently 
different from 2D flow about an infinite straight wing. Three-dimensional effects include 
spanwise flow, leading edge and wing tip vortices, down wash effects, and complex 
separation patterns at high angle” of attack as the wing approaches stall. As a result, the 
stall behavior of the 2D infinite span wing will in general be different from that of the 
finite span wing with the same airfoil section. 

• Aerodynamic per ormance of finite swept wings depends on aspect ratio (AR), taper ratio 
(X), sweep angle (A), Reynolds number (Re), wing twist, and potential fuselage/wing 
interaction effects. Note that the 3D wing was tested with a fuselage like body, as shown 
in figure 2-18 of this report. 

From experimental work and classical aerodynamic theory , a number of simple equations have 
been developed for relating the lift slope of a straight finite wing to that of a 2D wing with the 
same airfoil section. For planar untapered swept wings, basic swept wing theory can be applied 
to correct the lift slone of straight wings for the effects of wing sweep as discussed in 
reference H-5. Sweep corrections account for the fact that, for swept wings only the velocity 
component normal to the wing leading edge is responsible for wing lift and surface pressures. 
The velocity component tangential to the wing leading edge is important only for the 
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determination of the frictional stresses on the surface. Corrections for sweep effects depend to 
some extent on the method used to convert a straight wing into a swept wing. In general, sweep 
is introduced by rotating or by shearing a straight w ing. If rotation is applied to sweep the wing 
(bent-back or yawed wing), then the airfoil section of the unswept wing is the same as the airfoil 
section normal to the leading edge of the swept wing. A sheared wing is obtained by shearing 
backward (or forw ard) every section of the unswept airfoil, leaving its shape and lateral position 
unchanged The swept finite GLC-305 wing was obtained by applying the shear method. Thus, 
the airfoil section of the straight 2D wing and the stream wise section of the swept wing were the 
same. However, the airfoil section, normal to the swept wing leading edge was slightly thicker 
than that of the standard GLC 305 profile. Note that the simple sweep theory requires the use of 
the lift data for the airfoil section normal to the wing leading edge. A 2D viscous flow analysis 
was conducted with the XFOIL code [H-6] to assess the difference in lift performance between 
the GLC-305 airfoil and the thicker airfoil section normal to the wing leading edge. The results 
showed that the airfoil normal to the wing leading edge had a slightly higher maximum lift 
coefficient but in general the lift performance of both airfoils was very similar over the linear lift 
range. Thus, it was decided to use the baseline GLC-305 lift data in all calculations presented 
below. 

Experimental lift curves for the 2D and 3D clean wing are compared in figure H-2. The main 
differences between the 2D and 3D wing lift curves include the following: 

• Considerably Higher lift slope for the 2D wing (i.e., higher lift at a given angle of attack). 

• Angle of zero lift is negative for the 2D wing but slightly positive for the 3D wing. 

• Maximum lift is considerably greater for the 2D wing. 

• Stall behavior for the 2D and 3D wings is considerably different A gradual stall is 
observed to take place for the 3D wing compared to a sharp stall for the 2D wing. 



FIGURE H-2. COMPARISON OF LIFT CURVES FOR THE CLEAN 2D AND 3D WINGS 
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In the following, the lift characteristics of the clean and iced GLC-305 finite wing will be 
estimated from the lift properties of the clean and iced 2D wing using simple equations 
developed for swept wings. It should be emphasized that stall behavior cannot be predicted by 


The method used to estimate the finite swept wing lift data from that of the 2D wing was as 
follows: 

• Apply standard sweep, aspect ratio, and compressibility corrections to estimate the 3D 
wing Lift slope. 

• Estimate effective twist angle for the 3D wing and use this value to correct the 3D wing 
geometric angle of attack. 

• Estimate the effect of Reynolds number on lift slope using experimental data and a 2D 
viscous analysis panel code. 

• Apply a simple formula to estimate maximum lift for the 3D wing from that of the 2D 
wing. 

Due to downwash and geometric twist effects, the lift curve slope of the 3D swept wing was 
reduced compared to that of the 2D wing model. Knowing the 2D lift slope (ao) and the angle of 
attack for zero lift (a.L-o), one can predict the lift curve of a finite swept wing from equation H-l 
[H-7] which takes sweep angle and aspect ratio into consideration. Note that all lift slopes in 
equation H-l are per radian. From reference 24, the 2D slopes (ao) of the clean and iced wing 
were 0.0973 and 0.0977 per degree, respectively. Corresponding lift slope values per radian 
were 5.574 and 5.598. From equation H-l, the lift slopes (CLa)of the clean and iced finite swept 
wings were calculated to be 4.114 and 4.127 per radian or 0.0718 and 0.072 per degree, 
respectively. 



where AR= 6.8, A^ = 22° and A4L = 0.185 for the swept finite wing tested. 

Typically, planar wings with infinite and finite span and identical airfoil sections have the same 
angle of zero lift (ai-«). The GLC-305 swept wing model, however, had a washout of 4° at the 
wing tip, which reduced the effective angle of attack of the 3D wing, thus ( 0 : 1 - 0)30 and ( < 21 - 0)20 
were not the same. For a finite wing with washout, (< xl = o )3 d should be the difference of (<*1=0)20 
and Gen, where 0eu is the effective angle of twist of the whole wing. Note that at zero lift, the 
downwash angle generated by a finite wing is practically zero, thus only wing geometric twist 
causes the 2D and 3D geometric angles of attack to differ. To estimate the effective twist angle 
of the finite swept wing, the span wise lift distribution of the 3D wing must be considered. Using 
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the span wise lift distribution, the spanwise location, where the wing local lift times the span of 
the wing is approximately equal to the total lift of the wing (integral under the curve), is 
determined. From the limited surface pressure data ob ‘’lined with the 3D wing, an approximate 
spanwise lift distribution was computed and the data was used to determine the spanwise 
location corresponding to the effective twist angle. This location was close to the 50% semispan 
for the angles of attack presented in figure H-3. The geometric twist at the 50% semispan station 
(05o) was -1.14°. Thus, the effective twist of the wing at low angles of attack was -1.14°. 
Therefore, the difference between the 2D and 3D geometric angles of attack at zero lift was 1.14° 
(the 3D angle of attack was 1.14° higher than the 2D geometric angle of attack). 



O 6 10 16 20 25 30 35 40 46 60 65 60 
SpHMist station (in) 

FIGURE H-3. SPANWISE LIFT DISTRIBUTIONS OF THE CLEAN SWEPT WING 

The effect of Reynolds number on the lift behavior is demonstrated in figure H-4 using the 
experimental data of reference H-l and computational data obtained with the XF01L code [H-6]. 
Both the experimental and analysis data indicate a significant increase in maxim um lift and angle 
of stall as expected. In addition, a small increase in lift slope is observed as Reynolds number is 
increased. Since the airfoil section for the 3D was constant from root to tip, the Re effects for the 
3D wing should be similar to that observed with the 2D airfoil. 

The experimental lift curves of the clean and iced wings are compared in figures H-2 and H-5 
respectively. Each figure provides the 2D and 3D experimental lift curves and the linear part of 
the lift curve few the 3D wing (line with diamond symbols) that was computed from simple wing 
sweep theory'. The results indicate that the estimated 3D wing slope is in good agreement with 
the experimental data of the 3D wing. 

Next, an estimate of the maximum lift coefficient for the 3D wing is computed using simple 
wing sweep theory. According to Hoemer [H-5], for a swept wing can be predicted from 
the maximum lift coefficient (C/mu, m) of a 2D wing using equation H-2. Again, the formula 
provided is strictly applicable to planar wings and assumes that the airfoil section, in this case 
GLC-305, and the flow arc normal to the wing leading edge. 

CL^ = Cl mxlD • cos 2 A c , 4 (H-2) 
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where A ^4 is the quarter-chord sweep angle, which was 25° for the 3D GLC-305 swept wing. 
Experimental data presented in reference H-5 indicate that equation H-2 underestimates CL.* for 
sweep angles greater than about 20°. 


§ 




FIGURE H-4. EFFECT OF REYNOLDS FIGURE H-5. COMPARISON OF LIFT 
NUMBER ON GLC-305 LIFT SLOPE CURVES FOR THE ICED 2D AND 

3D WINGS 


Lift slopes, maximum lift coefficients, and angles of stall for the 2D and 3D clean and iced wing 
models arc summarized in table H-l. The stall angles p~ovided in table H-l were obtained from 
the experimental data. Due to span wise flow, leading edge and tip vortex structures, washout, 
and downwash effects, the stall mechanism of finite swept wings is considerably different from 
that of straight 2D wings and in general the angle of stall for 3D wings is higher. 

TABLE H-l. SUMMARY OF LIFT CURVE SLOPES AND CL ^ 



Clean Wing 

Iced Wing 

Lift curve slope 
(per degree) 

2D Exp 

0.0973 

0.0977 

3D 

Exp. 

0.0765 

0.0775 

Calc* 

0.0718 

0.0720 

CLniix 

2D Exp 

1.0850 

0.9702 

3D 

Exp. 

0.8738 

0.8955 

Calc.* 

0.8912 

0.7969 

<w(deg) 

2D Exp 

10.54 

13.23 

3D Exp 

13.78 

15.80 


* Estimates for 3D wing obtained from the 2D experimental data of reference H-l and simple 
sweep wing theory. 
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H.3 RELATIONSHIP BETWEEN 2D AMD 3D WING GEOMETRIC ANGLE OF ATTACK. 


To estimate the effective angle of attack of the 3D wing, the geometric twist and downwash 
angle effects need to be estimated and the geometric angle of attack must be adjusted for these 
effects. These corrections will provide an equivalent geometric angle of attack for the 2D wing. 
Note that for the same lift, the geometric angle of attack of a finite wing is greater than that of an 
infinite span wing. According to Hoemer [H-5], the characteristic of pressure distributions in the 
vicinity of 50% semi-span of a swept wing is the same as that found in conventional irfinite span 
straight wing sections. Therefore, by matching pressure distributions at 50% semispan of the 
clean 3D wing model to the coefficient pressure distributions of the clean 2D wing model, the 
relationship between the geometric angles of attack of the 2D and 3D wing models can be 
estimated. Note that this method assumes the Reynolds numbers of the 3D and 2D flows are the 
same, which is not the case here. However, as it was shown in figure H-4, the effect of Re over 
the linear poition of the lift curve was small. The results obtained by applying the pressure 
matr.hing method are presented in figure H-6. By fitting a least-squares straight line to the data 
presented in figure H-2, the relationship between the 2D and 3D geometric angles of attack can 
be obtained. Note that the difference (a 20 - am) provides an estimate of the average geometric 
twist and downwash angle at 50% semispan. 



FIGURE H-6. RELATIONSHIP OF GEOMETRIC ANGLES OF ATTACK OF THE TWO- 
DIMENSIONAL AND THREE-DIMENSIONAL CLEAN WING MODELS 


H-9 



H.4 REFERENCES. 


H-l. Addy, H.E., Broeren, A.P., Zoeckler, J.G., and Lee, S., “A Wind Tunnel Study of Icing 
Effects on a Business Jet Airfoil,” ALAA Paper 2003-727, January 2003. 

H-2. Ruff, GA, “Analysis and Verification of the Icing Scaling Equations,” AEDC-TR-85- 
30, Vol. 1 (Rev), March 1986. 

H-3. Armand, C., Charpin, F., Fasso, G., and LeClerc, G., “Techniques and Facilities Used at 
the ONERA Modane Centre for Icing Tests,” AGARD-AR- 127, Appendix A6, 
November 1978. 

H-4. Charpin, F. and Fasso, G., “Icing Testing in the Large Modane Wind-Tunnel on Full- 
Scale and Reduced Scale Models,” NASA TM-75373, March 1979. 

H-5. Hoemer, I.S.F. and Borst, H.V., “Fluid Dynamic - Lift,” Hoemer Fluid Dynamics, P. O. 
Box 342, Brick Town, NJ, 1975. 

H-6. Drela, M., “XFOIL 6.6 User Primer,” MIT Aero and Astro Engineering, March 19%. 

H-7. McCormick, B. W., “Aerodynamics Aeronautics and Flight Mechanics," 2 nd edition, John 
Wiley and Sons, 1995. 


H-10 



